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PROPELLANT HEATING AND DOSE TO CREW DURING ENGINE OPERATION 

by Gerald P. Lahti 
Lewis Research Center 

SUMMARY 

A FORTRAN IV computer program, QADHD, based on the LASL code QAD4, was 
written to evaluate the propellant heating in a nuclear rocket stage and to calculate the 
gamma and fast neutron doses received by a crew from the radiation emitted by the nu- 
clear engine during the operating period. The QADHD is a combination of three codes, 
QAD4, QADH, and QADD. The LASL QAD4 is a general purpose program written to 
calculate gamma and fast neutron doses due to a volume -distributed source in a complex 
source-shield geometry that is generally describable by quadratic surfaces. QAD4 is a 
point-kernel code that employs infinite -medium buildup factors to calculate the gamma 
dose and employs the Albert -Welton kernel to calculate the fast neutron dose. The QAD4 
included here was revised slightly from the LASL version to facilitate its use in the 
QADHD, as well as to simplify data input. 

QADH utilizes the QAD4 routines to evaluate local heating throughout the propellant 
(and tank wall) and then integrates these heating rates over the entire volume of the pro- 
pellant. Assuming that this heat is distributed uniformly throughout the propellant by 
natural convection (the complete-mix model) permits the rate of temperature rise of the 
propellant to be determined. QADH then determines the total temperature rise and accu- 
mulative heat input to the propellant as a function of reactor operating time. In some 
nuclear rocket configurations, the propellant may serve to shield the crew from radia- 
tions emitted by the engine. In these cases QADD utilizes the QAD4 routines to calculate 
biological dose rates to the crew during engine operation while the propellant depth in the 
tank (or amount of shielding the propellant provides) changes with time. QADD then inte- 
grates these dose rates over engine operating time to obtain a total dose-time history. 


INTRODUCTION 


In nuclear rocket vehicles, some special problems that require evaluation are the 


propellant heating and the dose to the crew from radiation emitted during the operation of 
the reactor. Excessive heating of the propellant may produce pumping problems so that 
either shielding, additional pressurization, or venting of the propellant may be required 
to avoid the problem. In a manned stage, the dose to the crew can control the shielding 
requirements. In many configurations, the propellant can serve as a shield between th$ 
crew and the engine. During engine operation, however, the consumption of propellant 
results in a continuous reduction of this shielding. 

The computer program QADHD was written to evaluate the heat deposition and tern - 
perature rise of the propellant and the dose to a crew during reactor operation. The 
code uses the Los Alamos Scientific Laboratory (LASL) point-kernel line-of- sight code 
QAD4 (ref. 1) as a basis for its calculations. QAD4 evaluates dose or heating occurring 
at a point due to a distributed volume source of radiation in a complex source -shield 
geometry that is generally describable by quadratic surfaces. QAD4 is used to deter- 
mine the local heating values throughout the tank, and the QADH portion of QADH then 
integrates by Gauss quadrature to obtain the total heating rate in the propellant. By 
assuming that convective currents distribute the heat uniformly through the propellant 
(the complete -mix model), the rate of temperature rise of the propellant in the tank 
can be determined. By relating the propellant volume and heating rates and temper- 
ature rise rates to the engine operating time, the code determines a heat-absorption and 
temperature-rise-time history of the propellant. 

QAD4 routines are used to determine the crew dose rate for a number of discrete 
propellant depths, and the QADD portion of QADHD determines the time for the propel- 
lant depth to change from one level to another from the engine run time, and then QADD 
integrates these dose rates over the total operating time to obtain the total dose in- 
curred. In a line-of-sight code such as this, the dose incurred from scattered gammas 
is estimated by using infinite -medium buildup factors. Implicit in these buildup factors 
is the assumption that all material in the neighborhood of the line-of-sight path experi- 
ences a similar flux. In a nuclear rocket stage, there may be substantial amounts of 
propellant and parts of the propellant tank and other structures outside the direct line of 
sight from source to detector, which are in very much higher flux fields than those in' the 
direct line of sight. From these regions, the scattered radiation received by the crew 
could be very much greater than that estimated by QADD; hence it is important to apply 
QADD (and QAD4) only to configurations where the calculations will be valid. Because 
of many similarities, the programs QAD4 (a slightly revised version of the LASL QAD4), 
QADH, and QADD were combined into a single code QADHD. The QAD4 included here 
has been revised slightly from the LASL version to facilitate its use in QADHD as well 
as to simplify data input. QADHD then may be used to perform three types of calcula- 
tions: dose or heating rates at detectors, propellant heating history, or crew dose his- 
tory. 
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To reduce data input normally required for the LASL QAD4, the version of QAD4 in- 
cluded here utilizes a library of gamma total cross sections, neutron-removal cross 
sections, parameters for infinite -medium buildup factors, and several gamma-response 
functions included internally in the code. 

v In general, in QADHD the geometry is that describable by quadratic surfaces identi- 
cal to the LASL QAD4 with a limitation that the emptying propellant tank (used in the 
QADH and QADD calculations) must be parallel to the z-axis (the axis of source sym- 
metry). Most nuclear rocket geometries fit this description. 


DESCRIPTION OF CODE 
QAD4 Method of Solution 


Computer output . - The QAD4 method of solution, which is also that used to calcu- 
late local heating and dose rates in QADH and QADD, is a line -of -sight technique that 
proceeds as follows. QAD4 divides the source region into a number of small source 
volumes and concentrates the total source strength of this volume at its geometric cen- 
ter. Knowing the complete geometry and material arrangement, QAD4 traces a ray 

from the center of each small source volume (the source point) to the receiver point and 

o 

keeps a tally of the geometric distance traveled and the mass thickness (g/cm ) of each 
element encountered on the traverse. 

The dose (or heating) rate from gamma rays is calculated as follows: 
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where 


* y < r o’ V Vj 
\K’ z o’Vo> 


D (r , z , cp ) 
y v 0 O *0' 
r un 


Vvv^ 


uncollided gamma dose rate (or heating rate) at point r , z Q , cp Q from 
source gammas of group I 

total (uncollided plus scattered) gamma dose rate (or heating rate) at 
point r , z q , cp Q from source of gammas of energy group I 

uncollided gamma dose rate (or heating rate) at point r , z Q , cp Q from 
all source gamma energy groups 

total (uncollided plus scattered) gamma dose rate (or heating rate) at 
point r Q , z Q , q> Q from all gamma energy groups 


(Units of D are rad/hr or other units depending on units of Kj. ) 



number of gamma energy group 
number of source point 

distance from J^ 1 source to receiver point, cm 

11. iL 

source strength of J in source point in I in energy group, MeV/sec 

output conversion factor for I th group for dose, (rad/hr)/ 
(MeV)/cm^-sec); for heating, (cal/g-sec)/(MeV/cm^-sec) 



where 


PI 


— ] (P^mj dimensionless e-folding length encountered between 


1L J.t_ 

J U1 source point and receiver point for I in energy group gammas 


M 

M 

\ p )m 

B i(bu) 


number of elements used 

1 L 1 L 

total gamma mass absorption coefficient for M l element, I U1 
energy group, cm 2 /g 

i.U J-U 

total mass thickness of M n element encountered between J Ln 

2 

source point and receiver point, g/cm 

infinite -medium buildup factor corresponding to energy group I and 
e -folding length bjj for appropriate representative material and 
receiver response 
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The neutron dose rate, estimated from the Albert -Welton kernel, is as follows: 


*' 


D (r , z , 
n v o’ o’ 



exp(-c)hj 



where 


D n ( r 0 ’ z 0 > dose rate from neutrons at point r Q , z q , (p Q , rem/hr 


J 

S J 


c 


where 



number of source point 

th 

source strength of J in source point, fissions/sec 



tween J 
hydrogen 


(pt) M j dimensionless fast neutron e-folding lengths be- 


source point and receiver point from all elements except 


2 th 

fast neutron removal cross section, cm /g for M element 

3 

equivalent thickness of hydrogen with density 0. Ill g/cm encountered 
between J** 1 source point and receiver point, cm 


The empirical constants are 


a l 

2. 187x10”^ cm^ (rem/hr)/(fission/sec) 

a 2 

0.29 (dimensionless) 

a 3 

0.83 (dimensionless) 

“4 

0.58 (dimensionless) 


Note that with this set of constants the Albert-Welton attenuation kernel will give 

meaningful results only if hj ^ 0. 5 centimeter. 

Output from the QAD4 section of QADHD is a table of r , z q , <p Q , D , D , D neut , 

7 un * 

IERR3, D , and D , as defined previously. If IERR3 = 1 on the output, no hydrogen 

y I y I 

un 
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was encountered on the traverse and the fi neu j- is in error. H IERR3 = 0, some hydro- 
gen was encountered; however, if hj is less than 0. 5 centimeter the D neut may still 
be in error. 

Library data . - To reduce data input normally required for the LASL QAD4, a li- 
brary of basic data used by QAD4 (and subsequently by all other segments of this code) f 
has been included in the program. This library includes for all elements from Z = 1 to 
Z = 92, the total mass absorption coefficients for gammas p/p. Tables of p/p against 
E are stored as FORTRAN BLOCK DATA for values of energy E of 0. 05, 0. 06, 0. 08, 
0.10, 0.15, 0.20, 0.30, 0.40, 0.50, 0.60, 0.80, 1.00, 1.50, 2.00, 3.00, 4.00, 5.00, 
6.00, 8.00, and 10.00 MeV. Data are from NBS circular 583 (ref. 2) and from LA-2237 
(ref. 3) without coherent scattering for the elements not in NBS-583. The gamma mass 
absorption coefficients for gamma energies listed as input for a given problem are then 
obtained from this master table by linear interpolation. 

A number of infinite -medium buildup factors are fitted to polynomials of the form 

B(b, E) = El A^E* 

i j 

or to other similar polynomials in APEX 510 (ref. 4). The bivariant coefficients A^. 
were included in the QADHD library for dose and energy absorption buildup in water, 
aluminum, iron, and lead. The selection of one of these buildup factors is input to the 
program. QAD type buildup -factor coefficients are calculated from these A^ for the 
given input gamma energies. 

Neutron removal cross sections S^/p (cm /g) are included for all elements from 
atomic number Z = 1 to Z = 100. These data are from reference 5. Cross sections 
for elements not in this reference were calculated from the fitted curves 


S R _ 0. 6 
p A 


A < 15 


= 0. 162 A" 0, 516 A > 15 


where A is the atomic number. 

The output conversion factor Kj, listed in the previous section, is assembled from 
two pieces of library data as follows. QAD calculates the energy flux E<t> at a receiver 

O 

point in MeV/(cm -sec). This flux is then multiplied by the mass energy absorption co- 
efficient P a /p (cm /g) and a numerical conversion factor to give the final desired output 
form: 
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E$ 


MeV 




a /cm 


* conversion 


Vcm -sec/ 


' rad/hr \ 
(MeV/g-secJ 



The library includes the mass energy absorption coefficients p /p for muscle (table 

r 

IA1 of ref. 6) that normally are used in QAD4 and QADD calculations, and it includes the 
mass energy absorption coefficients p a /p for hydrogen (ref. 6) that normally are used 
in QADH calculations. The tables of p„/p against E are stored for the same values 
of E as indicated for the p/p tables. Again, values for input gamma energies are in- 
terpolated from this master table. 
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For dose rate, the numerical conversion factors in the library are 5. 76756x10 
(rad/hr)/(MeV/g-sec) and for heating in liquid hydrogen with a density of 0. 07 gram per 
cubic centimeter 2.67557X10" 5 (cal/cm** sec)/(MeV/g-sec). The QAD conversion factor 
Kj is then the appropriate product, Kj = conversion x (p a /p) j. Although these P a /p 
and conversions are internal to the program, the user may load alternate ones if desired. 
The constants in the semiempirical Albert -Welton fast -neutron-dose kernel are 


ag = 0. 29 
a g = 0. 83 
a ^ = 0. 58 

The leading coefficient a j actually consists of the product of two terms, a "flux -to - 
dose" conversion a Q and a part to normalize the attenuation curve a*, such that the 
extrapolation of the hydrogen attenuation kernel 

a 2 , a 
a*h expl-a^h 

approaches 1.0 as h approaches 0. 

-3 

For neutron dose, values of a ^ stored as library data are a^ = 2. 187X10 
cm^ (rem/hr)/(fission/sec) and for neutron heating of liquid hydrogen are 
c*! = 9. 14x10"*^ cm^ (cal/cm^-sec)/(fission/sec). 

The value of 9. 14x10" 14 for calculating neutron heating in liquid hydrogen must be 
regarded as a rough estimate only. This value was obtained in the following manner 
from the used for dose calculations: 

Assume: 



1 tissue rad of fast neutrons = 10 rem 
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1 water rad = 1 tissue rad 


All energy deposited by fast neutrons in water is due to recoils of hydrogen in water * 

Therefore, i 

1 rem = 0. 1 tissue rad 
= 0. 1 water rad 

= 10 erg/g water (since 1 rad = 100 erg/g) 

Q 

= 10 erg/cm water 

O 

= 90 erg/g hydrogen (in l^O) (since density of hydrogen is 0. Ill g/cm in HgO) 
Therefore, 

a -2 19 10~ 3 cm2 ( rem / hr ) v 90 erg / g H x hr y 0. 239 cal y g H 
1 (fission/sec) rem 3600 sec 1Q 7 erg cm 3 

= 9. 14 10- 14 cm 2 (<=al/cm 3 -sec) 

(fission/sec) 


QADH Method of Solution 

Computer output . - QADH calculates the temperature history of the propellant during 
reactor operation. The QAD4 routines compute local heating rates, and the QADH inte- 
grates these over the propellant volume to obtain the total heating rate as a function of 
the quantity of propellant in the tank. If it is desired, the contribution to propellant heat- 
ing by energy deposited by gammas in the tank wall and tank bottom may also be evalu- 
ated. Assuming that convective currents distribute this heat uniformly throughout the 
propellant permits the temperature rise rates to be determined. When propellant volume 
is related to engine operating time, the code determines a heat absorption and tempera- 
ture history of the propellant. The computation proceeds as follows. The propellant 
tank is divided into a number of axial slices, and the Z -coordinates of the slices are 
loaded as input, as indicated in sketch (a). A number of receiver points are located in 
each slice, the number of receivers being determined by the order of Gaussian integra- 
tion (quadrature) specified for each of the three coordinate directions. This number may 
be 2, 4, 8, or 16 in the R-direction, 2, 4, 8, or 16 in the Z -direction (for each slice), 
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( 3 ) 

and 1, 2, 4, 8, or 16 in the -direction. (Only one point in the ^-direction is required 
when the tank and source are symmetrical about the same axis. ) These receivers are 
located at abscissas corresponding to the order of Gaussian quadrature specified. Thus, 
as few as 4 or as many as 4096 receivers are located in each slice. At each receiver 
point, the local heating rate is calculated as indicated in the section QAD4 Method of 
Solution. The total heating rate in the slice is obtained by completing the Gaussian quad- 
rature. 

The heating rate for the i slice is obtained from 


q. = f Zi+l f f w(R, Z, <p ) R dR dZ dcp 
1 J Z i J R *V 

by Gaussian numerical integration where W(R, Z, q>) is the local heating rate from neu- 
trons or gammas in cal/(cm^)(sec). The first page of QADH output is a table of 

1 z i < z i + i- z i> VOL i DTME i \ \ itj 


where 


VOL. 

DTIMEj 


1U o 

volume of i slice, cm° 

time to expel volume of i^ slice from tank, sec 

f Vi 

cal/sec deposited by gammas in i slice of propellant 
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q„ cal/ sec deposited by neutrons in i** 1 slice of propellant 

II* 



If wall heating is calculated, a table similar to that just given is printed out, and is fol- t 
lowed by a single line of output calculated for bottom heating. The terms q. in this 
table are the heating rates in the wall slices (which correspond to the propellant slices). 
These values are then summed by code with the q^ obtained for the propellant to obtain 
the total heating in each slice of propellant. 

The total heating rate into the propellant when the propellant level is at any slice co- 
ordinate is obtained by summing up the heating rates of all the slices from the tank bottom 
up to that level. When it is assumed that this heat is then redistributed by convective 
currents uniformly throughout the propellant by turbulence within the propellant (the 
complete-mix model), then the rate of temperature rise is determined by using the den- 
sity and specific heat and volume of propellant up to that level. Similarly, rates of heat 
input and temperature rise are computed for each position considered. 

These data are presented in the next pages of QADH output as a table of 


1 \ Sd Sii Si ^yi *ni Si 

where 

cal/sec being added by gammas to entire tank when liquid level is 
cal/sec being added by neutrons to entire tank when liquid level is Z^ 

Si Si + Si 

T^ rate of temperature increase of propellant, °C/sec from gamma heating when 
liquid level is Z^ 

T n i rate of temperature increase of propellant, °C/sec from neutron heating when 
liquid level is Z ^ 

T t l 


^yi ^ 

1=1 
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1=1 


T . = 

y i 


Qyj 


P c p £ VOL, 


1=1 



Qni 

p Cp £ VOL, 


where p is the density of liquid hydrogen, which is assumed constant (0. 07 g/cm 3 ), and 
Cp is the specific heat of liquid hydrogen, assumed constant (2.25 cal/(g °C). 

From the time required for the propellant level to drop from one position to the next 
and the arithmetic average of the rate of temperature rise and rate of heat input over this 
interval, an incremental temperature increase and heat input is evaluated for this time 
period (DTIMEj). 

These data are presented in the next page of QADH output as a table of 
i ACVi AQn! AQa AT rl AT ni AT,, 

where 

AQ^ amount of heat from gammas added to propellant in time DTEME^, cal 

AQni amount of heat from neutrons added to propellant in time DTEME^, cal 

AQ ti AQ yi + AQ ni 

A T^i temperature rise of propellant from gammas in time DTJME-, °C 

AT n i temperature rise of propellant from neutrons in time DTIME^, °C 

AT ti AT yi + AT ni 

and where 


aQ ri = 5 % + <W * DTME i 
a <5ni - \ «5ni + <W * MIME, 
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aT ri - \ (t yi + W * DTIME i 
AT ni - \ <T„i + W . DTIMEj 

V 

Finally, to accumulate these values of incremental temperature rise and heat input and to 
present them as a function of reactor operating time, a new index j is defined as 
j = NSLICE + 1 - i. Thus, at any time TIMEj, where TIMEj_j = DTIME.^g^jcj; and 

TIMEj = TIME. j + DT I ME i = NSLICE+l-j’ the fluid level is at Z i=NSLICE+l-j' TIME j 
is the time from start of engine operation. 

The final page of output is then 

J TIMEj Q,j Qhj Qq T vj T nj Tq 


where 

heat absorbed in propellant from gammas from start of reactor operation to time 
TIMEj, cal 

heat absorbed in propellant from neutron from start of reactor operation to time 
TIMEj, cal 

S'] + Sij 

temperature rise in propellant due to gammas from start of reactor operation to 
time TIMEj, °C 

temperature rise in propellant due to neutrons from start of reactor operation to 
time TIMEj, °C 

T • + T . 

Yi n] 

initialization, it was assumed that 

Si = AQ y, i=NSLICE 

Sil = AQ n, i=NSLlCE 

T yl = AT y, i=NSLICE 
T nl = AT n, i=NSLICE 


Q • 
^y] 


Sij 

Si 

T yj 


t . 
n] 


T ti 

For 
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Then the following computations were made: 

* Si = Si-1 + A S i=NSLICE+l-j 

( Sj = Sj-1 + AQ n, i=NSLICE+l-j 

T ri = T yj-i + AT r,i=NSLiCE+i-j 
T nj = T nj-1 + AT n, i=NSLICE+l-j 

These accumulated increments then give the temperature rise history and total heat- 
added history. 

It was assumed that all heat added to the propellant goes into raising the temperature 
of the propellant. If, however, one has a tank of propellant at a temperature near satu- 
ration at the start of operation, the temperature will increase only to saturation and then 
boiling will begin. For such a case, the total temperature rise calculated by QADH will 
of course have no real meaning. However, from the total heat -added -time history output 
of QADH and from the heat of vaporization (108 cal/g of LHg), sufficient information is 
available to calculate the boiloff after saturation temperature has been attained. 

Intermediate output . - If it is desired, intermediate output (local heating rates and 
results of various integrations) may be obtained, which may be of interest if regions of 
high local flux are anticipated (e. g. , away from the shadow of a partial shadow shield). 
Control for this output is by the parameter NSKIP that is loaded as data on the tank 
emptying card QH1. The intermediate output is as follows. If NSKIP £ 3, write 

H R,Z ,cp W y /Buildup W y W n 


where 


n 

R, Z ,<p 
W y /Buildup 



receiver point number 
receiver coordinates, R and Z, 
local heating rate at R, Z, cp due 

local heating rate at R, Z, cp due 
buildup factor) 

local heating rate at R, Z, cp due 


cm, and cp, rad 
to uncollided gammas 

to all gammas (uncollided rate times 
to neutrons 


If NSKIP s 2, write 
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and 


f w v d r, 

^(p Y 


f W d cp; f (W + W)d<p; 

J(p •'cp r 



dcpR dR; 



dcpR dR; 





ff NSKIP s 1, write 


f Z i+l [ fy/ d< p R dR dZ; /* Z i+l /* f w n d<pR dR dZ; /' Zl+1 f W +W n )d(pRdRdZ 
^z. •'R'V r •'Z. •'R'V n y Z. ^R'V r 


All values are printed as soon as they are calculated. The integrals are of interest only 
in the checking of calculations. 


QADD Calculations 

Method of solution . - QADD calculates the time -integrated dose at a detector during 
the operating time of the reactor when the propellant level in the tank between the source 
and detector is dropping. The tank is divided into a number of axial slices and the Z- 
(axial) coordinates of these slices Z (I) are loaded as input. Gauss quadrature is used to 
evaluate the volumes of these slices, and the total tank volume is also determined. From 
the total tank volume and the propellant flow rate, the time for the propellant level to 
drop to any of the slice coordinates is determined. In the code, a boundary describing 
the propellant level is designated MOVER, and this boundary takes on successive values 
of the slice coordinates. The dose rate is calculated, as indicated in the preceding sec- 
tion, now with the MOVER positioned at each axial coordinate. For each discrete propel- 
lant level I, one then has total neutron and .gamma dose rates £> n (I) and 6^(1). The dose 
obtained while the propellant level drops from one position to the next lower position is 
equal to the product of the time elapsed during this level change and the average dose 
rate during the level change, or 

ADOSE (I) = I |D (I) + D (I + 1)] x ATIME(I) 

ADOSE n (I) = I [D n (I) + D n (I + 1)] x ATIME(I) 

2 
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where 


•ADOSE y (I) gamma dose obtained while level drops from level (I + 1) to I, rad 
D (I) gamma dose rate with propellant at level (I), rad/hr 

i 

i> y (i+i) gamma dose rate with propellant at level (I + 1), rad/hr 
ATIME(I) time for propellant level to drop from level (I + 1) to I, sec 
ADOSE n (I) neutron dose obtained while level drops from level (I + 1) to I, rem 
D n (I) neutron dose rate with propellant at level (I), rem/hr 

D n (I + 1) neutron dose rate with propellant at level (I + 1), rem/hr 

QADD also evaluates the accumulative dose as a function of engine operating time: 


DOSE y (I) = 


NSLICE 

E 


J=I 


ADOSE y (J) 


NSLICE 

DOSE n (I) = £ ADOSE n (J) 

J=I 


NSLICE 

TIME (I) = J ATIME(J) 

J=I 

where 

TIME (I) time elapsed while propellant moves from initial level to level I, sec 

DOSE y (I) total gamma dose obtained from time = 0 to TIME (I), rad 
DOSE n (I) total neutron dose obtained from time = 0 to TIME(I), rem 
NSLICE total number of tank slices considered 

Computer output. - The computer output from the QADD section QADHD is a table 
of 

I Z(I) D y (I) ADOSE y (I) DOSE y (I) 

D n (I) ADOSE n (I) DOSE n (I) TOTDOS(I) 
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where TOTDOS(I) is DOSE^(I) + DOSE n (I). It should be reiterated that the dose calcu- 
lated here is good only if the large propellant tank is fully shielded from the source which 
results in slowly varying fluxes in the propellant. If high flux regions are present in 
some parts of the tank, unnoticed by the line of sight, then the dose may be underesti- 
mated by orders of magnitude. , 


DATA INPUT DESCRIPTION 
Source Specification 

The source volume is subdivided into a desired number of smaller volume elements. 
The source strength of each of these is then represented by a point source at the center 
of the volume element. The source is specified by the total source strength So, 
the location of subdividing surfaces, and a relative source strength distribution. In 
QADHD the source must have cylindrical symmetry. The subdividing surfaces of the 
source are specified by inputting a set of R, Z, and <p; that is (R^, Rg, • • RjJ, 

(Zj, Z2, . . ., Zjyj), and (<pj, The source will then be contained in the 

volume Rj to R^, Z^ to Z M , and cp x to cp^. Point sources will be located at the 
center of each of the orthogonal volumes described by these surfaces. 

The relative spatial source strength in the R- and Z -directions is then specified in 
one of two ways (no cp variation is allowed). Values can be given for and 

T) g hi the cosine distribution 

P(R, Z) = A cos |j(R - £ 2 ) cos tJj( Z - tj 2 ) 

where this expression is valid in the range Z^ ^ Z ^ Zj^ and Rj < R < R^. The con- 
stant A is determined by the program by a normalizing definition, 

S = / Rl / Zm /*^ N P(R, Z)R dR dZ d <p 
° J z x 'cp 1 

An alternate means of specifying relative source strength is to give a relative power on 
the R- and Z -coordinates previously given. The relative power at each point need only 
be consistent for its particular coordinate. Again, normalization similar to that just in- 
dicated will be performed by the program. The source strength for each equivalent point 
source will then in essence be calculated from 
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P(R, Z)R dR dZ d(p 



1 where Sj mn is the source strength at point ton and P(R, Z) is the normalized relative 
source distribution. The source as specified does not have to fill a geometric QAD re- 
gion and it is not necessary for the specified source to be in only one region; it may ex- 
tend over several. 

As used in QADHD, the units of the total source S Q are in fissions per second or 
captures per second to correspond to a convenient gamma spectrum given in MeV per 
fission or MeV per capture (the final gamma source strength is MeV/sec) and an a j in 
the neutron dose kernel having units in dose rate per fission per second. However, other 
units may be used provided that dimensional consistency is maintained through the calcu- 
lation. 

Geometry 

The geometric description of problems is the same as that used in the LASL QAD4. 
The geometry description is illustrated by the reactor-shield-propellant tank configura- 
tion shown in figure 1. This configuration could represent a stage of a space vehicle that 
uses a clustered set of nuclear rocket engines, where, because of the symmetry of the 
arrangement, only one engine, its associated propellant tank, and an adjacent propellant 
tank are sufficient for calculational purposes. It is assumed that the crew compartment 
is not in this stage and that propellant heating is an important problem in the stage, hi 
figure 1, the entire reactor assembly is represented as a carbon core, the shadow shield 
is composed of zirconium hydride, and the tanks contain liquid hydrogen propellant. 

Boundaries . - All space is divided into regions bounded by quadratic surfaces that 
are described by equations of the following types: 

Type 1: 

AX 2 + XXo + BY 2 + YYo + CZ 2 + ZZo - K = 0 

Type 2: 

A(X - Xo) 2 + B(Y - Yo) 2 + C(Z - Zo) 2 - K = 0 

Type 3: 

(X - Xo) 2 + (Y - Yo) 2 - K = 0 
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X - K = 0 


Type 5: 

I 

Y - K = 0 


Type 6: 


Z - K = 0 


A, B, C, Xo, Yo, Zo, and K are constants for each surface specified in the input. Types 
4, 5, and 6 are X-, Y-, and Z -planes, respectively. Type 3 is a circular cylinder 
parallel to the Z-axis. Types 1 and 2 may be used to obtain other geometric shapes. 

For example, a cone with its axis parallel to the Z-axis (sketch (b)) is described by 

(X - Xo) 2 + (Y - Yo) 2 - (Z - Zo) 2 tan 2 <* = 0 



An oblate (or prolate) spheroid with its axis parallel to the Z-axis (skecth (c)) is de- 
scribed by 


(X 


Xo) 2 + (Y - Yo) 2 + R 2J 


(Z - Zo) 


( Z l - Z J 




The configuration shown in figure 1 is arranged into a set of QAD regions by the use 
of bounding surfaces, as shown in figure 2. These bounding surfaces are presented in 
the following table. 


Boundary, 

number 

Equation 

QAD 

equation 

type 

Parameters 
for input 

1 

tsi 

ii 

o 

6 

o 

n 

2 

Z * 130 

6 

K = 130 

3 

Z = 160 

6 

K = 160 

4 

Z = 500 

6 

K = 500 

5 

Z = 750 

6 

K = 750 

6 

Z = 1600 

6 

K = 1600 

7 

X 2 + Y 2 = (50) 2 

3 

Xo = 0, 
Yo = 0, 

K = 2500 

8 

X 2 + Y 2 = 600 2 

3 

Xo = 0, 

Yo = 0, 

K = 360000 

9 

(X - 1250) 2 + Y 2 = 600 2 

3 

Xo = 1250, 
Yo = 0, 

K = 360000 

10 

X 2 + Y 2 = 2000 2 

3 

Xo = 0, 

Yo = 0, 

K = 4xl0 6 

11 

X 2 + Y 2 = (Z - 150) 2 tan 2 45° 

2 

A = 1 
B = 1 
C = -1 
Xo = 0 
Yo = 0 
Zo = 150 
K = 0 

12 

(X - 1250) 2 + Y 2 = (Z - 150) 2 tan 2 45° 

2 

A = 1 
B = 1 
C = -1 
Xo = 1250 
Yo = 0 
Zo = 150 
K = 0 


Regions . - Regions are defined as either positive or negative with respect to any of 
their boundaries. K f (X, Y, Z) = K is a boundary of region R, then let a quantity called 
the residual be defined as r(X, Y, Z) = f (X, Y, Z) - K. If a point X Q , Y Q , Z Q is in 
region R, then the geometry routine requires that the product j * r (X Q , Y Q , Z Q ) must be 
negative. The factor j is equal to +1 or -1 depending on whether r(X Q , Y Q , Z Q ) is 
negative or positive. If j = +1, then the region is considered positive with respect to the 
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(d) 

any point in region Rg, such as Z Q = 15. For 

r < X C V Z o> - '«o- Y 0' Z 


boundary, or if j = -1, the region is 
negative with respect to the boundary. 
The sign of j (+ or -) is affixed to the* 
QAD boundary number in the input QAD 
region description. For example, con r 
sider a region Rg bounded by planes 
Z = 0 (boundary Bl) and Z = 20 (bound- 
ary B2), shown in sketch (d). Consider 
boundary Bl (type 6), K = 0; hence, 

J - K = Z - 0 = 15 
o o 


The residual is positive; thus j must be -1 and Rg is negative with respect to Bl. 
For boundary B2, K = 20, and 


r(X rt , Y n , Z J = f(X„, Y . Z ) - K = Z - 20 = 15 - 20 = -5 
O O 0 o’ o’ o o 

Since the residual is negative, j = +1 and Rg is positive with respect to B2. Similarly, 
it is seen that the region within a cylinder is positive with respect to the cylindrical sur- 
face. 

Another method of assigning signs to boundaries, which yields the same correct end 
result, is to determine, if a particle in region R crosses a Z -boundary, the sense of its 
motion in Z (+ or -). The boundary is then given that sign. Similarly for cylinders and 
cones, if the particle is going in the positive direction of r, then the radial boundary 
is +. 

Regions are classified as either inside or outside regions. Once a ray from a source 
point passes into an outside region, it is assumed that it will encounter no other regions 
on its way to the receiver. Regions 11, 12, and 13 in figure 2 are outside regions. The 
input required for each region includes (1) the number of boundaries associated with the 
region (NBNDZN). NBNDZN is given a (+) sign if the region is am inside one and a (-) 
sign if the region is an outside one; (2) the composition number for the region (NCMPZN); 
(3) the boundary number of each boundary associated with the region, to which has been 
affixed the appropriate sign (+ or -) that defines whether the region is positive or nega- 
tive with respect to that boundary; and the probable region entered after crossing a given 
boundary. On the region card (read in with a Format (7(15, 14)), columns 1 to 5 contain 
±NBNDZN, 6 to 9 contain NCMPZN, 10 to 14 contain ± boundary number, 15 to 18 contain 
the probable region entered when crossing this boundary, 19 to 23 contain ± another 
boundary number, 24 to 27 the probable region entered when crossing this boundary, and 
so forth for up to six boundaries per region. The region numbers are assigned in the 
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order that the region cards are read in. For example, the region card for region R2 in 
sketch (c) is 


Columns 

1-5 

6-9 

10-14 

15-18 

19-23 

24-27 

+2 

Composition 

number 

-B1 

R i 

+B2 

R 3 


The material composition of each region is specified by assigning each 
composition a number and constructing a matrix with the appropriate densities of all ele- 
ments used in the problem, as shown in the following sample problem. If an element is 
not used in a composition, its density in that particular composition would be zero. The 
composition number is then part of the region specification. For the sample problem, then: 


Regions in which composition 
is used 


Composition 


3 

Density, g/cm 

Hydrogen, 

Carbon, 

Zirconium, 

i H 

6 C 

40 Zr 

0 

1.5 

0 

.11 

0 

5.49 

.07 

0 

0 

0 

0 

0 


(1) Core 

(2) Shield 

(5, 7, 8, 10) Liquid hydrogen 

All others Void 


The complete set of region cards for the sample problem (fig. 2) is then 


Column 


1-5 

6-9 

10-14 

15-18 

19-23 

24-27 

— 

28-32 

33-36 

37-41 

42-45 

46-50 

51-54 

73-80 

Number of 

Composition 




• Boundary, Probable region 




Region 

boundaries 

number 












+3 

1 

-1 

13 

+7 

3 

+2 

2 





1 

+3 

2 

-2 

1 

+7 

3 

+3 

4 





2 

+4 

4 

-1 

13 

+10 

12 

+3 

4 

-7 

1 



3 

+3 

4 

-3 

2 

+10 

12 

+4 

5 





4 

+3 

3 

-4 

4 

+11 

6 

+5 

8 





5 

+5 

4 

-4 

4 

+10 

12 

+5 

9 

-11 

5 

-12 

7 

6 

+3 

3 

-4 

4 

+12 

6 

+5 

10 





7 

+3 

3 

-5 

5 

+8 

9 

+6 

11 





8 

+5 

4 

-5 

6 

+10 

12 

+6 

11 

-8 

8 

-9 

10 

9 

+3 

3 

-5 

7 

+9 

9 

+6 

11 





10 

-1 

4 

-6 

8 









11 

-3 

4 

-10 

3 

-1 

13 

+6 

11 





12 

-1 

4 

+1 

1 









13 
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QADH and QADD Tank-Emptying Geometry 


General . - In addition to the general geometry described in the previous section, 
additional parameters are required to define the mode of tank emptying. QADD and 
QADH propellant tanks are constrained to be circular in cross section (in the X-Y plane), < 
and the tank centerline must lie parallel to the Z-axis and in the <p> - 0 plane. The tank 
may consist of only one part or two parts, a tank head and a main tank body. The sur- 
faces describing the tank must be fully rotated type 2 or type 3 QAD boundaries. Several 
possible tank shapes are illustrated in figure 3. The QAD boundary numbers denoting the 
tank head and main tank body are specified in input as JHEAD and JTANK, respectively. 
Figure 3(c) shows that there are two ways to specify a one-piece tank; in this case the 
specification of JTANK = 0 or JHEAD = 0 is a flag to indicate a one -part tank and is not 
a QAD boundary number. 

To simulate the propellant tank emptying, the tank is divided into a number of axial 
slices (NSLICE). The Z -coordinates of these are loaded as input. The i of that Z 
which divides the tank head from the body is designated KNUCKL. The volume of each 
slice is determined by Gauss quadrature and the total tank volume determined. From 
this and the reactor operating time (BURN), the flow rate and level of propellant in the 
tank as a function of time are determined. Up to this point, QADD and QADH function in 
the same manner. 

QADH calculations . - The information given in the previous section is sufficient for 
QADH to calculate propellant heating, as indicated in the section QADH Method of Solu- 
tion. However, if gamma heating in the tank wall is to be considered as part of the pro- 
pellant heating, then the QAD boundaries defining the outside wall of the tank head and 
the main tank body must be specified. These are designated JWHED and JWTNK, re- 
spectively. (The inside walls are assumed to be the previously inputted JHEAD and 
JTANK, respectively. ) Also, the Z-coordinate of the bottom of the tank wall must be 
specified to calculate the contribution from the tank bottom. 

In the sample problem illustrated in figure 2, for heating calculations of tank I, 
JHEAD = 11 and JTANK = 8. For heating calculations of tank II, JHEAD = 12 and 
JTANK = 9. For both tanks, select NSLICE = 8, Z = 500, 550, 600, 650, 750, 900, 

1100, 1300, and 1600 centimeters; therefore, KNUCKL is 5. For tank I, concentric 
with the core, no variation in local heating in the <p -direction is expected; hence, set 
NPHI equal to 0 because only one point in this direction {cp = 0) is necessary for heating 
calculations. 

QADD calculations . - For the QADD calculation of dose received while the tank is 
emptying, a QAD boundary, defining the propellant level in the tank, is designated as 
MOVER. To simulate emptying, this boundary takes on the successive values of the Z - 
coordinates of the slice divisions. The dose rate is determined with the MOVER at its 
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initial position and each of the NSLICE values of which are input. (Note that only 
NSLICE values of Z ^ are required for QADD while values of NSLICE + 1 are required 
‘for QADH. ) The MOVER is the boundary that separates the liquid in the tank from the 
emptied portion of the tank. Thus, it is important to distinguish between the MOVER 
. boundary and the upper boundary of the tank and to insert between them a region of void 
or pressurizing gas. As the position of MOVER changes, this region increases in size. 
For example, consider the configuration in figure 4, which represents a stage with a 
single nuclear engine and which contains the crew compartment (a wall of which is shown 
as region 10 in the figure). The dose to the crew is an important problem in the opera- 
tion of this stage. The reactor, shield, tank sizes, composition, and arrangement are 
similar to those illustrated in figure 2. Note that the tank has been extended a small 
amount to distinguish between the upper tank boundary 9 and the MOVER boundary 6 and 
that region 8 inserted between these boundaries is considered to be a void region. The 
dose to the crew located at a position Z = 2100 centimeters, and R = 0 centimeters is 
desired. 

For this QADD calculation, JHEAD = 11, JTANK = 8, and MOVER = 6. Select a set 
of values of Zj = 500, 550, 600, 650, 750, 900, 1100, and 1300 centimeters (NSLICE = 8) 
for fluid levels. Also KNUCKL= 5. (Z 5 = 750 is the plane which divides the tank head 
from the main tank body. ) 

One additional piece of information is necessary for proper execution of the calcula- 
tion. The MOVER boundary 6 ultimately will be the upper boundary of region 5 rather 
than boundary 5. To take care of this situation, boundary 6 is also specified as a bound- 
ary of region 5. For this case then, the QAD region card for region 5 must be 


Card column 

1-5 

6-9 

10-14 

15-18 

19-23 

24-27 

28-32 

33-36 

37-41 

42-45 

+4 

3 

-4 

4 

+11 

6 

+5 

7 

+6 

8 


The boundary condition (+6 8) is redundant until the MOVER enters the tank head. Then 
(+5 7) becomes redundant. 


QADHD DATA INPUT 
Program Notes 

Data input to QADHD consists of (1) a control -card selecting mode of calculation for 
QAD4, QADD, or QADH; (2) a set of data, called general QAD data, that describes 
source, geometric regions and boundaries, composition and source spectrum; and (3) 
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receiver cards for QAD4, tank geometry and receiver cards for QADD, or tank geometry 
cards for QADH. The detailed version of the general QAD data, although virtually the 
same as that for the LASL QAD, is presented here for the sake of completeness. 

In its present form, the QADHD program is designed to handle 50 regions, 50 bound- 
aries, 40 compositions, 20 elements, 30 energy groups, and 20 source points in each of < 
the three coordinates, R, Z, and <p (a total of 8000 source points). With these limits, 
the program does not occupy the maximum 32K storage available in the core. The user, 
by increasing the dimensions of any or all of these arrays, may alter any of these limits. 

O 

A subroutine WEIGHT was added to calculate the volume (cm ) and weight (g) of all 
inside QAD regions defined by planes QAD type-3 cylinders and QAD type-2 quadratic 
surfaces fully rotated about an axis parallel to the Z-axis. In a general purpose program 
such as this, most of the computation time is spent in tracing rays through the various 
regions. There is one ray traced from each source point to each receiver point, and the 
total running time was observed to be proportional to the number of these rays. QADHD 
calculates 3500 to 5000 source-to-receiver rays per minute on the IBM-7094 n, depend- 
ing on the complexity of the geometry. 
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QADHD Input Flow 


Read MASTER card 

Card 

column 


FORMAT (15, 14) 


1-5 


6-9 


MASTER = 1 QAD4 problem 
= 2 QADD problem 


NERGY 


Is NERGY > 0? 


> 


Yes 


1A 


No 


IB 


Read general QAD 
data. 



= 3 QADH problem 

< 0 Use output conversion factors internal to code 
(listed on p. 7) 

> 0 Read new output conversion factors 

o 

Read gamma energy absorption coefficients (cm /g); 
20 values corresponding to energy grid for gamma 
mass attenuation coefficients listed on page 6. 
FORMAT (8E9. 4) 

Read DOSE CONVERSION FORMAT (2 E 9. 4) 

Card 

column 

1-9 Dose conversion factor for gammas 
10-18 for Albert -Welton kernel 

Comment 
Control 
Source cards 
Region cards 
Boundary cards 
Element cards 
Composition matrix 
E for gamma energy groups 
Source spectrum for gammas 


fSee pages 26-29 


QAD4, page 30 


Read receiver cards. 


-? ...i 

[ 

=3 

=2 

1 

QADD, page 31 

'"Read QADD cards. 
Then go to 1 . 


QADH, page 32 


''Read QADH cards. 
Then go to 1 . 




Symbol means read a card or set of cards as described. 
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General QAD Data For QADHD 


The following block of data is required for all problem types (QAD4, D, H) and de- 
scribes the source and geometry of the problem. 


✓ 

G1 



G2 


Card column 


Read any alphabetic, numeric information FORMAT (12A6) 
Card column 1 *~72 

Read control card FORMAT (6(15, 14)) 

Variable 


1-5 

6-9 

10-14 

15-18 

19-23 

24-27 

28-32 

33-36 

37-41 

42-45 

46-50 


51-54 


LSO 

Number of source points on R coordinate axis 

MSO 

Number of source points on Z coordinate axis 

NSO 

Number of source points on <p coordinate axis 

MAT 

Number of elements used 

NCOMP 

Number of compositions 

NREG 

Number of regions 

NRGY 

Number of gamma energy groups 

NBOUND 

Number of boundaries 


Not used 

NZSO 

Most probable region for source 

ISRC 

1 power defined by cosine distribution 


P(R, Z) = cos £j(R - £ 2 )cos t^(Z - r^) 

2 power defined by point -by-point tabulation 
0 use source defined by previous case 
INEUT >0 make neutron calculations 

<0 bypass neutron calculations 



ISRC=0, go to G10. 


j 


Otherwise 



Read source card. 


Card column Variable 


FORMAT (5E9. 4) 


1-9 

10-18 

19-27 

28-36 

37-45 



^2 


Total source strength, fissions/sec or captures/ sec 

Parameters of cosine power distribution for ISRC=1 
(not used if ISRC=2) 
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G4 



G5 



G6 


R- coordinates of source volumes 

(LSO+1) values 

Z -coordinates of source volumes 

(MSO+1) values 

<p- coordinates of source volumes 

(NSO+1) values 


(0 — for QADH concentric tank) 
(0 — ^271 for QADH offset tank) 


c 


If ISRC * 2, go to G10. 


D 


Otherwise read 



G7 


s' 

G8 



G9 


Relative value of source at R-coordinates 
of card G4 

Relative value of source at Z -coordinates 
of card G5 

Relative value of source at ^-coordinates 
of card G6 

(Note: Regardless of values on this card, 
QADHD will assume no variation of power 
in ^-direction; the contents of this card 
are ignored. ) 


s' 

G10 


Region cards (NREG cards required) 

See pages 19-21 for more detail. The first 
card read is region 1, the second card is 
region 2, and so forth. 


FORMAT (8E9. 4) 
FORMAT (8E9.4) 
FORMAT (8E9.4) 


FORMAT (8E9. 4) 
FORMAT (8E9.4) 
FORMAT (8E9.4) 


FORMAT (7(15, 14)) 


Card column 


1-5 


6-9 

10-14 

15-18 

19-23 

24-27 

28-32 

33-36 

37-41 

42-45 

46-50 

51-54 

55-59 

60-63 


± number of boundaries for this region; if region is an outside region, the 
sign is -. 

composition number for this region 

± boundary number; region is + or - with respect to boundary 
probable region entered if this boundary is crossed 
± boundary 
probable region 
± boundary 
probable region 
± boundary 
probable region 
± boundary 
probable region 
± boundary 
probable region 


5“ Up to six boundary surfaces may be used* to define a region. 
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Boundary cards (NBOUND cards required) FORMAT (15, 14, 7E9. 4) 

Card column 

Boundary number 

Equation type (1) AX 2 +XXo+BY 2 +YYo+CZ 2 +ZZo-K = 0 
(general quadratic) 

(2) A(X -Xo) 2 +B(Y - Yo) 2 +C (Z - Zo) 2 -K = 0 
(spheroid, conic, or cylinder) 

(3) (X-Xo) 2 +(Y - Yo ) 2 -K = 0 
(circular cylinder parallel to Z -axis) 

(4) X-K = 0 (X -plane) 

(5) Y-K = 0 (Y-plane) 

(6) Z-K = 0 (Z -plane) 

A 
B 
C 
Xo 
Yo 
Zo 
K 

NBLD, Element card FORMAT (8(15, 14)) 

Card column 
1-5 


6-9 

10-14 

15-18 


MATZ(MAT) Atomic number of MAT th element used 


NBLD= 1 Water dose buildup factor used 

2 Aluminum dose buildup factor used 

3 Iron dose buildup factor used 

4 Lead dose buildup factor used 

5 Water energy absorption buildup factor used 

6 Aluminum energy absorption buildup factor used 

7 Iron energy absorption buildup factor used 

8 Lead energy absorption buildup factor used 
MATZ (1) Atomic number of first element used 
MATZ (2) Atomic number of second element used 
MATZ (3) Atomic number of third element used 


1-5 

6-9 


10-18 

19-27 

28-36 

37-45 

46-54 

55-63 

64-72 
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Composition matrix cards 


FORMAT (8E9. 4) 


G13 


o 

Enter the table of density (g/cm ) of elements. 
At least one card per composition. The first 
card (cards) read are for composition 
number 1, and so forth. 


Card column 


1-9 Density of element 1 in this composition 

10-18 Density of element 2 in this composition 

19-27 Density of element 3 in this composition 



G14 

Gamma energy card 

FORMAT (8E9. 4) 



A list of the mean gamma energies (MeV) 
for each gamma energy group. 

NRGY values required. 



G15 

Source spectrum, gammas 

FORMAT (8E9.4) 


A list of the MeV/fission or MeV/capture 
in each gamma energy group. 

NRGY values required. Units to be 
consistent with S Q . 

(S Q ) * (Source spectrum) must be MeV/sec 






QAD4 Receiver Cards 


Read only if MASTER = 1. 



FORMAT (3E9.4) 


(Note: The R-coordinate also functions as a signal 
for the ending of QAD4 case for proper return of 
the program for a new case. After the calculation 
for the final receiver is made, another receiver 
card with R < 0 is required for proper return to 
the main program. ) 





QADD Cards 


Read only if MASTER = 2. 

Tank emptying card 


QD1 


FORMAT (3(15, 14), E9. 4) 


Card column 


1-5 


6-9 

10-14 

15-18 

19-23 

24-27 

28-36 


QD2 


QD RECVR 


G 


Card column 

1-9 

10-18 

19-27 


Variable 

MOVER 

NZ 

NSLICE 

KNUCKL 

JHEAD 

JTANK 

BURN 


QAD boundary number of top of propellant 

at beginning of engine operation 

1, 2, 3, or 4; 2 points will be used to 

get the volume of each tank slice 

Number of tank slices used to simulate emptying 

The KNUCKL*'* 1 Z -slice divides tank head from body 

QAD boundary number of tank head 

QAD boundary number of tank body 

Engine run time, seconds 


Z- slice coordinate card 


FORMAT (8E9.4) 


List of Z -coordinates of fluid level in tank used to simulate tank 
emptying. Begin with smallest value of Z (coordinate of tank 
bottom generally). NSLICE values required. The top of the 
propellant Z (NSLICE + 1) is taken as the Z -coordinate of the 
QAD boundary MOVER. Z (KNUCKL) divides the tank head 
from the tank body. 


Is R coordinate > 0 ? 


Receiver card 


R- coordinate of receiver 
Z -coordinate of receiver 
(^-coordinate of receiver 
No r 


FORMAT (3E9.4) 


Yes 


Make QADD 
calculation. 


Finished with 
QADD case. 
Return to read 
new master. 


(Note: Make sure that the QAD region specification (card G10) 
for the tank head includes the redundant specification for the 
MOVER boundary as indicated on p. 23. ) 







QADH Cards 


Read only if MASTER = 3. 


Read tank emptying card. 


FORMAT (4(15, 14), E9. 4, 15) 


Card column 


NPHI = 0, 1, 2, 3, or 4. The number of Gaussian abscissas in 

the fluid in the ^-direction is Z NPHI . 

NR = 1, 2, 3, or 4. The number of Gaussian abscissas in the 

fluid in the R-direction is 2^. 

NZ = 1, 2, 3, or 4. The number of Gaussian abscissas in the 

fluid in the Z -direction per slice is 2** Z . 

NSLICE The number of tank slices used to simulate tank emptying 
NSKIP = 0. Print final values only (normal option). 

>1. Also print out integrals over the slice as soon as 
they are calculated. 

>2. Also print out integrals over <p and R as soon as 
they are calculated 

>3. Also print out local heating rates at all points. 
KNUCKL The KNUCKL th ZSLICE coordinate divides the tank head 
and body. 

JHEAD The QAD boundary number of the tank head 

JTANK The QAD boundary number of the tank body 

BURN The engine burn time, seconds 
JWALL = 1. Make wall heating calculation. 

* 1. Bypass wall heating calculation. 

(Note: The total number of receivers in the propellant is 
NSLICE *2< NPHI+NR+NZ >.) 


Read tank ZSLICE coordinates. 


FORMAT (8E9. 4) 


A list of Z -coordinates of propellant level is used to simulate tank 
emptying. Begin with the smallest value of Z (coordinate of tank 
bottom). NSLICE + 1 values are required. The last coordinate is 
the top of the tank. Z(KNUCKL) divides the tank head from the tank body. 


Make QADH 
calculation of 
heating in fluid. 


Is JWALL = 1 ? 


Obtain time- 
temperature 
history. 


Return to main 
program for new 
MASTER. 


Go to next page 
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QW1 


Head wall heating control card 


FORMAT (15, 14, 15, 219, 14, E9. 4) 


Card column 

1-5 NPHI =0, 1, 2, 3, or 4. The number of Gaussian abscissas in the 

wall in the ^-direction is 2 NPHI . 

6-9 NR = 1, 2, 3, or 4. The number of Gaussian abscissas in the 

wall in the R-direction is 2 NR . 

10-14 NZ = 1, 2, 3, or 4. The number of Gaussian abscissas in the 

wall in the Z -direction is 2 NZ . 

15-23 NSKEP = 0. Print final values only (normal option). 

>1. Also print out integrals over the slice as soon as they 
are calculated. 

>2. Also print out integrals over R and cp as soon as they 
are calculated. 

>3. Also print out local heating rates at all points. 

24-32 JWHED The outer QAD boundary number of wall of tank head (inner 

boundary is assumed to be JHEAD) 

33-36 JWTNK The outer QAD boundary number of wall of tank body (inner 

boundary is assumed to be JTANK) 

37-45 ZPRIME Z-coordinate of bottom of tank wall 



(Note: Sample problem n calculates wall heating. See p. 34 for details. ) 
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EXAMPLES OF DATA INPUT 


A complete set of data input for the calculation of propellant heating in tank I of fig- * 
ure 2 is presented in table I. Part of the input was specified in previous sections of this 
report. To complete the required data input, it is assumed that the total operating power, 
is 3. 3x10*9 fissions per second (1000 MW) and the relative power distribution is constant 
in the R- and Z -directions. (The flat power distribution is conveniently given by specify- 
ing a cosine power distribution with coefficients and ? 7 j both equal zero. ) The 
source is subdivided into volumes bounded by the surfaces R = 0, 15, 30, 40, 50 centi- 
meters, Z = 0, 30, 60, 90, 105, 120, 130 centimeters, and <p = 0, 1. 57, 3. 14 radians. 
Thus, there are (4 x 6 x 2) 48 radiation source points defined for this case. 

Thirteen gamma energy groups were selected for this sample problem; the mean 
gamma energy for each group is listed on cards G14 of table I; the gamma source 
strength (MeV/sec)/(fission/sec) corresponding to these energy groups is listed on 
cards G15. 

Two-point Gauss quadrature has been selected for both R- and Z -integrations of 
heating in a tank slice. Because of symmetry, only one point is used in the (^-direction. 
Thus, there are (2x2x1) four receivers located in each tank slice (more are generally 
recommended in the R-direction). Eight tank slices are used; therefore, a total of (8x4) 
32 receivers are spotted throughout the tank. 

The total number of geometry rays traced by QAD is equal to the product of the num- 
ber of source points and the number of receiver points. In this case, (48 x 32) 1536 such 
rays are drawn. Since the IBM 7094 n can handle 3500 to 5000 of these per minute, an 
execution time of less than 1/2 minute is expected. 

A second sample problem is illustrated in figure 5. This calculation includes tank 
wall heating. Complete data input are presented in table II. 
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TABLE I. - FORMAT FOR CALCULATION OF PROPELLANT HEATING IN AXIAL TANK I OF FIGURE 2 



50 51 52 55 54 55 56 57 58 59 60 61 62 65 61. 65 *6 6 7 68 6 9 70 71 72|75 74 75 76 77 7B 79 80 





























































TABLE H. - FORMAT FOR CALCULATION OF CORE GAMMA HEATING IN TANK WITHOUT SHIELD 



5 5 56 5 7 58 5 9 6 0 61 62 6 3 6 <. 65 66 67 68 69 70 71 72 73 7 k 75 76 77 78 79 80 


































TABLE H. - Continued. FORMAT FOR CALCULATION OF CORE GAMMA HEATING IN TANK WITHOUT SHIELD 
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TABLE n. - Concluded. FORMAT FOR CALCULATION OF CORE GAMMA HEATING IN TANK WITHOUT SHIELD 























FORTRAN LISTING AND SAMPLE PROBLEM OUTPUT 


The complete IBM 7094-11 FORTRAN IV listing of the QADHD program is now pre- * 
sented along with the computer output of QADHD for the propellant heating calculation in 
tank I of figure 2 for which input data were previously described. 


SIEFTC QAOMST DECK 

COMMON / DOSE/ XSECE 0 ( 20 , 2 ) , ALF(2), DC0N(2) 

COMMON /STDATA/ 

1 NAME (100) , ENRG (20) , XSECO (20 , 100 ) , XS ECE A ( 20 ) , 

2 C (4,6,9) , EMI N33 ( 8) , NA MB ( 8 ) , XSECNT(IOO) ,DCONST 

COMMON /Q ADATA / 

1 NSET (16) , R SO (21) , ZSO (21) , PHISO (21) , 

2 NBNDZN(50) , NCMPZN(50) , LBD (6,50) , NTRYZN(6,50) , 

3 NEQBO (50) , ABO (50) , BBD (50) , CBO (50) , 

4 XOBD (50) , YOBO (50) , ZOBD (50) , DBD (50) , 

5 HY0RAK30) , XSECN (30) , XSECG (20,30), BILDO (30) , 

6 BILD1 (30) , 8 IL02 (30) , BIL03 (30) , GAMEN (30) , 

7 CONV (30) , COMP (20,40), EBAR (30) , MATZ (20) , 

8 FL (21) , FM (21) , FN (21) , 

9 ASO , XI1S0 , XI2S0 , ETA1S0, ETA2S0, A LF1 , ALF2 , ALF3 , 

X ALF4 , EPSLN , FUDGE , NBLO 

10 FORMAT ( 15, 14) 

11 FORMAT ( 8E9.4 ) 

101 READ( 5,10) MASTER , NE RG Y 
GO TO ( 1,2,3) , MASTER 

1 J =2 

IF(NERGY) 20,20,30 

2 J =2 

IF(NERGY) 20,20,30 

3 J = 1 

IF(NER GY) 20,20,30 
20 DO 25 1=1,20 
25 XS ECEA ( I)=XSECEB(I,J) 

AL F 1=ALF( J ) 

DCONST=OCON( J) 

GO TO 6 

30 READ( 5,11 ) ( XSECEA ( I ) , 1=1,20) 

READ( 5,11 ) DCONST, ALF 1 

6 GO TO (7, 8,9), MASTER 

7 CALL QAD4 
GO TO 101 

8 CALL Q ADO 
GO TO 101 

9 CALL QADH 
GO TO 101 
END 
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$ I BFTC Q4QAD DECK 


SUBROUTINE QAD4 
COMMON /3 ADA TA / 


1 

NSET 

( 16) , 

R SO 

(21) 


zso 

(21 ) 


t 

PHIS 0 

(21 ) 

2 

N BNDZ N (50) , 

NGMPZN 

(50) 


LBD 

(6,5 

0) 

* 

NTRYZN ( 6 » 50 ) 

3 

NEQBD 

<50 , 

ABD 

( 50) 


BSD 

(50) 


t 

CBD 

(50) 

4 

XOBD 

150) , 

Y03D 

( 50) 


ZO0D 

(50) 


t 

DBD 

(50) 

5 

HYDRATC30) , 

XSECN 

(30) 


XSECG 

(20, 

30) 

t 

BILDO 

(30) 

6 

BILD1 

(30) , 

3IL32 

(30) 


BILD3 

(30) 


♦ 

GAMEN 

(30) 

7 

CON V 

( 30) , 

COMP 

(20 >40) 


EBAR 

(30) 


* 

MATZ 

(20) 

8 

FL 

(21) 

FM 

(21) 


FN (21) 


t 



9 

A SO 

, XI ISO 

, XI2S0 

, ETA1 

SO, ETA2S0, 

ALF 1 

, ALF2 , A 

X 

ALF4 

, EPSLN 

, FUDGE 

, NBLD 








COMMON 

/GAUSSQ/ 

V( 4 f 1 6 

)» W(4, 

16) 






COMMON 

/3 ADMAN / 










1 

SREG 

( 20 ) , 

AGAM 

(30) 

» 

DGAM 

(30) 


t 

ID 

(12) 

2 

UDGAM 

( 30) , 

UAGAM 

(30) 









EQUIVAL 

EN CE ( L SO 

,NSET( 1) ) , (MSO 

* 

NSET ( 2 ) ) , ( NSO, 

NSET (3) 

). 


ALF3 , 


I(NRGY,NSET( T) >,< ISR:,NSET( II) ) * ( NRC VR ,NSET ( 9 ) ) 

4 FORMAT ( 8E 9 .4 ) 

5 FO RM A T ( 1 1 9HO RECEIVER R(RCVR) Z(RCVR) PHKRCVR) 

10T G MAX TOT G NEUTR OOSE ER GROUP MIN G DOSE MAX G 

6 FORMAT ( IH0I6,5X1P6E 12.4, I 2 , 1 5 , 5 X2E 1 2 . * / ( I 9 1 , 5 X2 E 12 .4 ) ) 

7 FORMAT ( 7HO RCVRI3,4H R=1PE12.4,4H Z=E12.4,5H PH I = F 1 2 . 4 , 1 1H 

1URCE PT3I3 ) 

10 FORMAT ( 1H 1 , 5X , 12A6) 

11 FORMAT ( 7H0 RCVRI3,4H R=1PE12.4,4H 
1URCE P T3 I 3 t 8H TYPEI2 t 29H GEOMETRY 

20 L IN 1 = 5 5 
11 = 0 
N ERR = 0 

21 CALL INPUT! ID( I) ) 

23 I F ( ISRC)25,2 5, 24 

24 CALL SOURCE 

25 READ (5 , 4 ) RRC , ZRC , PH I R: 

I F ( RRC )999,27, 27 

27 11=11+1 


Z=E1 2 • 4 #6 H PHI= El 2* 
ERROR AT THIS POINT) 


QADC0 170 
Q A DC 0 1 8 0 
Q ADC0330 
MIN TQA0C0340 
DOSE) 0 A DC 0 3 5 0 
QADC0360 
S0QADC0370 
0ADC0380 
0ADC0390 
4, 1 1H S OQ ADC0400 
0 A DC 0 4 1 0 
0 A DC 04 20 
QADC0430 
0 A0C0440 
0 A DC 0 4 5 0 
Q A DC 04 60 
Q A DC 04 70 


QADC0500 


XRC=RRC <t CO S( PH IRC I 
YRC=RRC*SIN(PHIRC) 
NSFTI 1 5 ) = 0 
I ERR 3 =0 
ON EUT = 0 .0 
DGAM T = 0.0 
UOGAM T =0 • 0 
DO 30 J = l , NR GY 
OGAM ( J )=0.0 
30 UDGAM ( J )=0 .0 
DO 60 N = 1 1 N SO 
SINS=S IN ( PHI SO (N ) ) 
CO S S = CO S( PHI SO (N) ) 
00 60 L = 1 » L SO 


0ADC0510 
0ADC0520 
0 ADC 0 530 
0 A DC 05 40 
Q ADC0550 
Q ADCC560 
0 A DC 05 70 
0A0C058D 
0 ADC0590 
0 ADC0600 
0 A DC 061 0 
0ADCC620 
0 ADC0630 
Q A0C0640 
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X $Q=R S0( L )*COS S 
YSO=RSO!L )*SINS 
FLN=FL !L ) *FN < N ) 

DO 60 M = 1 , M SO 

NS ET( 1 5) =N SET! 15)+1 

I F( NSET( 15 )- NSE T { 13) ) 36,34,33 

33 IFCNSFTI 15J-NSET! 14) >35,35,36 

34 WRITE (6, 10) ( IDII ) , I = 1 , 12 ) 

L IN 1-5 5 

35 WRITE I6,7)II,RRC,ZR:,PHIRC,L,M,N 
CALL SLITE ( 1> 

36 FLMN-FLN*FM( M ) 

39 CALL LENGTH! SR EG ( l ) , S SO , XSO , YSO , ZSO I 
IF! I ERR 1 ) 40, 40 , 73 

40 CALL KERNEL! SREG! 1 ) , SSO , A NE UT , AG AMI 
I ERR 3=MAX 0! I ER R 3, 1 ERR 2 ) 

A=FLMN/$S0**2 
DN EUT = DNE UT+A* ANEUT 
00 45 J = 1 , NR GY 
OGAM! J )=DGAM ! J ) +A* AG A M ( J ) 

45 UDGAM! J )= UDGAM (J)+A*UAGAM(J) 

60 CONTINUE 

ONEUT=ONEUT/ 12 ,5663706 
DO 65 J = 1 , NR G Y 

A = GAM EN( J )*CONV( JI/12. 5663 706 
DGAM ( J )=A*DGAM ( J ) 

UDGAM ( J ) = A* UOGAM! J ) 

UOGAM T=UDGAM T+UDGA M( J ) 

65 DGAMT = DGAMT+OGAM( J ) 

67 IF! 55-LIND68, 68,70 

68 WRITE (6, 10MIDU), 1=1,121 
WRITE 16,5) 

L IN l^NRGY 

70 WRITE !6, 6)11, RRCtZRS, PHI RC,UDGAMT 
1GAMI J ) , DGAM ( J ) , J-l , NRS Y) 

IF( MSETI 9 1 ) 71,71,707 

707 PUNCH 708, R RC , ZRS , ONE UT 

708 FORMAT! 2F 8 . 2 , 1 PE 10. 3 > 

PUNCH 709, ! DGAM ( J ) , J- 1 ,NRG Y) 

709 FORMAT I IP 8E9 *2 ) 

71 L IN l-L IN 1 +NRGY+ 1 
GO T3 25 

73 WRITE !6, 11)11 ,RRC , ZRC , PH I RC , L , M, N 
N ERR=N ERR 4- 1 
L IN 1-5 5 

IF! 2-N ERR ) 74, 74, 25 

74 READ i 5 ,4)XX 
IF! XX ) 999 , 74, 74 

999 RETURN 
END 


Q A DC 0 6 5 0 
0 ADC0660 
0A0C0670 
Q A DC 0 6 8 0 
0ADC0690 
0 A DC 07 00 
Q A0C0710 t 

0A0C0730 

0 ADC0750 
QADC0760 

) ,XRC,YRC ,ZRC, IERR1) 0ADC0770 

QADC0780 

fUAGAMd) , I ERR2 ) 0ADC0790 

0ADC0800 
0 A DC 08 10 
QADC0820 
0A0C0830 
0ADC0840 
0 A DC 0 8 5 0 
0 A DC 0 8 6 0 
0ADC0870 
0 ADC0880 
0 A DC 0 8 9 0 
0 ADC0900 
QADC0910 
Q ADC0920 
Q A DC 0 9 3 0 
0 ADC0940 


0 ADC0970 

DGA MT ,DNEUT , IERR3, ! J,UD 

QADC0990 


Q ADC 1000 
0ADC1010 

IERR1 

0ADC1030 
GADC1040 
0 ADC 1050 


0 ADC 1080 
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$ I BFTC QDQAD DECK 


SUBROUTINE Q ADO 
COMMON /0 ADA TA / 

1 NSET 1161 , R SO (21) , ZSO (21) , PHISO (211 , 

2 NBNDZN(50> , NCMPZN{50) , LBD (6,53) , NTRYZN<6,50) , 

3 NEQBD (50) , A3D (50) , BBO (50) , CBD (50) , 

4 XOBO (50) , Y03D (50) , ZOBD (50) , DBD (50) , 

5 HYDRATI 30 ) , XSECN (30) , XSECG (20,30), BILDO (30) , 

6 BILD1 (30) , BIL02 (30) , BIL03 (30) , GAMEN (30) , 

7 CONV (30) , COMP (20,40), EBAR (30) , MATZ (20) , 

8 EL (21) , FM (21) , FN (21) , 

9 ASO , XIISO , XI2S0 , ETA1S0, ETA2S0, ALF1 , ALF2 , ALF3 , 

X ALF4 , EPSLN , FUDGE , NBLD 

COMMON /GAUSSQ/ V(4,16), W(4,16) 

COMMON /O ADHD/ 


1 

Z 

( ICO) 

t 

vz 

(16) 

9 

wz 

(16) 

9 

VR 

( 16) 

9 

2 

WR 

< 16) 

» 

VP 

(16) 

9 

WP 

(16) 

9 

VOL 

(LOO) 

9 

3 

DT IMF 

( 100) 

9 

RANGE 

(100) 

9 

vv 

(100) 

9 




4 

CTIME 

(100) 

9 

DEPTH 

(100) 

9 

DGOT 

(ICO) 

9 

DNOT 

(100) 

9 

5 

DAVE 

(100) 

9 

ONAVE 

(100) 

9 

OELG 

(100) 

t 

DELN 

(100) 

» 

6 

CUMG 

( 100) 

9 

CUMN 

UOO) 

9 

TOTAL 

(100) 

9 

DUMMY 

(1600) 


C0MM3N 

/Q ADMAN / 










1 

SR EG 

(20) 

9 

AGAM 

(30) 

9 

DGAM 

(30) 

9 

ID 

(12) 

9 

2 

UDGAM 

(30) 

t 

UAGAM 

(30) 









EOUIVALENCE(LSO,NSET( 1)) , ( MSO , NSET ( 2 ) ) , I NSO, NS ET { 3 ) ) , QA0C0170 

1(NRGY,NSFT( 7) ) , ( ISR: ,NSET( 11) ) , ( NRC VR , NSET ( 9 ) ) QADC0180 

2 FORMAT ( 14H0 FLOW RATE = 1PE12.4.9H CM3/SEC/16H TANK VOLUME = 

1E12.4, 5H CM3) 

3 FORMAT ( 20HL QAOO C ONTROL , 61 6 , 1 PE12 . 4 ) 

4 FORMAT ( 8E9.4 ) 

5 FORMAT ( 1H I 6, 3X1PE 12. 4 ,36X ,El 2 . 4 ,1 2 X ,E12 . 4 ) 

6 FORMAT! IH 16 , 3 X 1P9E 1 2. 4) 

7 FORMAT ( 7H0 RCVRI3,4H R=1P612.4,4H Z=E12.4,6H PHI= E12. 4, 1 1H SO 

1URCE PT3I3 ) 

8 FORMAT ( 3( 15, I4),E9.4) 

9 FORMAT ( 27H0 RECEIVER R, Z, PHI * 1 P3 El 2. 4/1HL62 X, 7 HGAMMA — 15X, 

X9HN EUTRQN — /118H I DEPTH ( I ) VOL(I) DTIME(I) 

Y T IME( I) DOSE RATE CUM. DOSE DOSE RATE CUM. DOSE TOTAL DO 

ZSE ) 

10 FORMAT ( 1H 1 , 5X , 12A6) 

11 FORMAT ( 7H0 RCVRI3,4H R=IPE12.4,4H Z=E12.4,6H PHI=E12.4, 11H S0QADC0400 

1URCE P T3 1 3, 8H TYP5I2,29H GEOMETRY ERROR AT THIS POINT) 0ADC0410 

PI = 3. 141593 
20 11*0 
NERR=0 

CALL INPUT! 1 0 1 1 ) ) 

I F( ISR C ) 25, 2 5, 24 

24 CALL SOURCE 

25 R E AD( 5 , 8 )MOVER ,NZ,NSLICE, KNUCKL , JHEAO , JT ANK, BURN 

WR ITE ( 6,3 ) MOVER ,NZ,NSLIC£, KNUCKL , JHE AD, JT ANK, BURN 

NGAUSZ =2**NZ 
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DO 252 I = 1 »NGA USZ 

vzu)=vhz,i ) 

252 WZ( I ) = W(NZ,I ) 

J J =NSL ICE + 1 

READ ( 5,4)<Z(I ) , I = l ,NSLI CE > 

Z< JJ )=DBD( MOVER) 

IF(JTANK) 27,27,26 

26 IFINEQBDI JTANK).EQ.3) A3D ( JTANK) =1 . 

27 IF(JHEAD) 29,29,28 

28 IF(NEQBD< JHEAD) .EQ.3) A3D ( JHEAD) =1 . 

CALCULATE VOLUME 

29 SUMV=0. 

00 295 NSL=1»NSLICE 
SUMVZ =0. 

RANGE( NSL )= Z ( NSL-t-1 )-Z (NSL) 

J B=JHEAD 

IF( NSL .GE .KNUCKL ) JB = JTANK 

00 293 KZ =1, NGAUSZ 

ZZ = RAN GE( NSL )* VZ ( KZ ) +Z < NSL) 

RR=S3RT( (DBD( JB)-CBO{ J3)*(ZZ -ZOBO(JB) ) ** 2 ) / ABO ( JB ) ) 
293 SUMVZ=WZ(KZ)*PI*RR*RR«- SUMVZ 
VOL(NSL)= RANGE (NSL )* SUMVZ 
295 SUMV=SUMV + VOL(NSL) 

FLOW=SUMV/BURN 
WR I TE ( 6, 2 ) FLOW, SUMV 
FF= 1 , / FL 0 W 
SUMT=0. 

DO 297 1 = 1, N SL ICE 

IJ=JJ- I 

VV( I ) = VOL 1 IJ ) 

OT IME f I ) = VV< I )*FF 
SUMT=SUMT+OT IME ( I ) 

297 CT IME I I)=SUMT 

298 READ (5,4) RRC , ZRC , PH I RC 
I F( RRC ) 999,299,299 

299 XRC=RR C*CO S( PH IRC ) 

YRC=RRC*SIN(PHIRC ) 

SUMG=0. 

SUMN=0. 

DO 99 NS=1,J J 
IJ=NSL ICE +2-NS 
OFPTH( NS ) = Z(IJ)-Z(1> 

OBDIMOVER )=Z( I J ) 

1 1 = 1 1 + l 
NSETI 1 5) = 0 
IERR3 =0 
ON EUT = 0.0 
DGAMT = 0.0 
UDGAM T=0. 0 
DO 30 J=l ,NRGY 
DGAMt J )=0 .0 

30 UDGAM ( J )=0.0 
DO 60 N=l ,NSO 
SINS=SIN(PHISO(N)) 

COSS=COS( PHI SO ( N ) ) 

00 60 L=1 ,L SO 


QA0C0610 
0 A0C0620 
0 A0C0630 
0A0C0640 


XSO=RSO(L )*COSS 
YSO=R SO! L )* S IN S 
FLN=FL(L)*FN(N> 

. DO 60 M=1 , M SO 

NSET( 1 5) =N SE T( 15)*1 
t F( 64 S E T t 1 5 )-NSET ( 13))36,34,33 

33 I FC N SFTt 15)-NSET( 14) > 35,35,36 

34 WRITE (6,10)(ID(I),I=1,12) 

35 WRITE (6,7)II,RRC,ZR:,PHIRC,L,M,N 
CALL SLITE ( I) 

36 FLMN=FLN*FM(M) 

CALL L BIGTH! SREG! 1) ,SS0,XS0,YS0,ZS0IM) ,XRC , YRC ,ZRC , IERR1 ) 

I F ( IERR1 ) 40, 40, 37 

37 WRITE (6, 11)11, RRC,ZRC,P-IIRC,L,M,N,IERR1 
N ERR =N ERR *1 

IFI 2-NERR >298,40,40 

40 CALL K ERNEL ( SR EG! I) *SSO,ANEUT , AG AMI 1) , UAGAMI I ) , I ERR 2) 

I ERR 3=MAX0! I ERR 3, IERR2) 

A=FLWN/SS0**2 

DNEUT =DNE UT+A*ANEUT 

DO 45 J=1 ,NRGY 

DGAM! J )=OGAM(J)+A*AGAM(J) 

45 UDGAM J)=UDGAM( J)+A*UAGAM( J) 

60 CONTINUE 

DNEUT=DNEUT/ 12. 5663706 
DO 65 J = 1 »NR GY 

A = GAM ENI J )*C0NV!J)/12. 5663706 
DGAMI J )=A*DGAM ( J ) 

UDGAM! J)=A*UDGAM! J ) 

UDGAM T =UD GAHT+UDGAM! J ) 

65 DGAMT=DGAMT+DGAM ( J ) 

DGDTI NS) = DGAMT 
ONDTI N S)=DNEUT 
I FI NS- 1 ) 70,70,69 

69 DAVE!NS) = (DGTZ +DGDTI NS ) ) *. 5 
DNAVE!NS)=!DNTZ+DNDT(NS) )*. 5 
T = DT I MF ( NS- 1) /3600. 

DELG(NS)=DAVE!NS)*T 
DELN(NS)=DNAVE(NS)*T 
SUWG=SUMG+DELG (NS ) 

SUMN = SUMN +DELN ! NS ) 

CUH G1 N S ) = SUMG 
CUMN(NS) = SUMN 
TOTAL ( NS ) = SUNG +SUMN 

70 DGTZ=DGDT (NS ) 

99 DNTZ = DNDT I NS ) 

WR ITE I 6, 1 0 ) ( ID I I ) , I = 1 , 12 ) 

WR ITEI 6, 2 ) FL OW, SUM V 
WRITE16.9) RRC ,ZRC,PHIR; 

1 = 1 

WRITE! 6,5) I,DEPTH(1),DG0T!1),DNDT!1> 

WR ITE! 6,6) I I , DEPTH! I > , VV! I -1 ) ,DTI ME ( 1-1 ) , CT I ME ( I- 1 ) , DGDT ( I), 
XCUWG1I) , DNDT ! I ) ,CUNN ! I ) ,TOTAL(I) ,I=2,JJ) 

GO TO 298 
999 RETURN 
END 


QA0C0650 
QA0C0660 
0A0C0670 
Q A0C0680 
QA0C0690 
QA0C0700 
0A0C0710 


QADC0750 

QADC0760 


0A0C0790 
QADC0800 
0ADC0810 
QA0C0820 
QADC0830 
QADC0840 
0ADC0850 
QADC0860 
QADC0870 
QADC0880 
0 A0C0890 
QADC0900 
QADC0910 
QADC0920 
0ADC0930 


45 


SIBFTC QHQAD DECK 


SUBROUTINE QAOH 
COMMON /Q ADA TA / 


1 

NSET 

(16) 

R SO 

(21) 

t 

ZSO 

(21) 


PHISO 

(21 ) 

2 

NBNDZN( 50 ) , 

NGMPZN 

(50) 

* 

LBD 

(6,50) 


NTRYZN(6 ,50 

3 

NEQBD 

(50) , 

A3 3 

( 50) 

t 

BBD 

(50) 


CBD 

(50) 

4 

XOBD 

(50) , 

Y08D 

(50) 

t 

ZOBD 

(50) 


D8D 

(50 ) 

5 

HYOR A T( 30 ) , 

XSECN 

(30) 

» 

XSECG 

(20,30) 


BILDO 

(30 ) 

6 

BILOl 

(30) , 

3ILD2 

(30) 

? 

BILD3 

(30) 


GAMEN 

( 30 ) 

7 

CON V 

(30) 

COMP 

(20 ,40) 

» 

EBAR 

(30) 


MATZ 

(20 ) 

8 

FL 

(21) , 

FM 

(21) 


FN (21) 




9 

A SO 

, XII SO 

» XI2SQ 

, ETA1S0 

1, ETA2S0, ALF1 

, ALF2 . 

X 

ALF4 

, EPSLN 

, FUDGE 

, N8LD 







COMMON 

/G A US SO/ 

V( 4» 1 6) f W ( 4 * 

16 

) 





COMMON 

/O ADHD/ 









1 

Z 

(100) , 

vz 

(16) 

t 

WZ 

(16) 

♦ 

VR 

(16) 

2 

WR 

(16) , 

VP 

< 16) 

» 

WP 

(16) 


VOL 

(100) 

3 

OTIME 

(100) , 

RANGE 

(100) 

f 

vv 

(100) 

t 



4 

SRG 

(16) , 

SRN 

(16) 

* 

SPG 

(16) 

♦ 

SPN 

( 16) 

5 

SG 

(100) , 

SN 

(100) 

? 

ST 

(100) 

* 

QNDOT 

(100) 

6 

QGOOT 

(ICO) , 

QTDOT 

(100) 

f 

TGDOT 

(100) 

t 

TN90T 

(100) 

7 

TTOOT 

(100) , 

DQG 

(100) 

t 

DON 

1100) 

» 

DQT 

( 100) 

8 

DTG 

(100) , 

DTN 

(100) 

t 

DTT 

(100) 

t 

OG 

(100) 

9 

ON 

(100) , 

QT 

(100) 

t 

TG 

(100) 

f 

TN 

(100) 

X 

TT 

(100) , 

T 

(100) 

t 

JTANK 

, JHEAD, 


KNUCKL 

, NS L I 

COMMON 

/Q ADMAN / 









l 

SREG 

(20) f 

AGAM 

< 30) 

t 

DG AM 

(30) 

* 

ID 

(12) 

2 

UDGAM 

(30) 

UAGAM 

(30) 








ALF3 


EQUI VALENCE ( L SO, NSET ( m t (MSO,NSET(2) ) » ( NSO , NS ET ( 3 ) ) , 
1( NRGY , NSE T ( 7 ) ) , ( ISRCtNSET(ll)),( NRC VR , NSET (9 ) ) 

2 FORMAT 1 14HO FLOW RATE = 1PE12,4,9H CM3 / SEC/ 16 H 
1E12.4,5H CM3) 

3 FORMA T ( 20H1 HEATING C ONTROL , 81 6 f l PE12 . W , I 5 ) 

4 FORMAT ( 8E 9 .4 > 

5 FORMAT (IH0,68X, 27HFLUI0 HEATING RATES CAL/SEC/7H 

l) Z(l + l)-zm VOLU) OTtMF(I) GAMMA 

2 TOTAL) 

6 FORMAT 4 1H I 6 f 5 X IP 7E 1 2. 4 ) 

7 FOR MAT < 7HO RCVRI3,4H R=1PE12.4,4H 
1URCE PT3I3 ) 

8 F0RMAT(4( 15,14) ,E9.4, 15) 

9 FORMAT (IP3EI2.4) 

10 FORMAT { lHlfSXf 12A6) 

11 FORMAT ( 7H0 RCVRI3.4H R=1PE12.4,4H 
1URCE PT3I3,8H TYPE12,29H GEOMETRY 


TANK VOLUME = 


I »13X, 76 HZ 
NEUT 


Z = E1 2 • ^ , 5 H PHI*EI2.4 f 11H 


Z = E1 2 • ^ , 6 H PHI=E12.4, I1H 
ERROR AT THIS POINT) 


12 

FORMAT ( 7HL 

Z 

= /(1P10E12.4) 


P I = 3. 

141593 



20 

11=0 





NERR = 

0 




CALL 

INPUT ( 

ID ( 

U) 


IF( ISRCI25, 

25, 

24 


QADC0170 

QADC0I80 


SOQAOCO 370 
Q A DC 0 3 8 0 


QADC0390 
S 0QADC0400 
0ADC04 10 
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24 CALL SOURCE 

25 REAO( 5,8) NR HI , NR , NZ , NSL I C E , N SKI P , KNUC KL , JHF AD , JT ANK , 8U RN, JW AL L 
WRITE! 6, 3) NPHI , NR , NZ , NSL l CE , NSKI P , KNUCKL , JHEAD, JT A NK , RURN, JWALL 
NGAUSR =2*4NR 

DO 251 1 = 1, NGAUSR 
VR! I) = V!NR,I ) 

251 WR( I ) = W< NR , I ) 

NGAUSZ =2< C *NZ 

DO 252 1=1, NGAUSZ 
VZ( I) = V! N Z , I ) 

252 WZ( I ) =W{ N Z » I ) 

IF! NPHI ) 253,253,254 

253 NPHI=0 
NGAUSP =1 
VP! 1) = 0. 

WP( 11=1. 

GO TO 256 

254 NGAUSP = 2**NPHI 

DO 255 1=1, NGAUSP 

VP! I ) = V! NP HI , I ) 

255 WP< I) = W«NPHI , IJ 

256 J J =NSL ICE + 1 

READ ! 5, 4 M Z ( I ) ,1=1, JJ) 

WR ITE ( 6, 12 )( Z ( I),I=l,JJ) 

IF! JTANK+JHEAD ) 999,999,259 
259 IF(JTANK) 27,27,26 

26 IFtNEOBD! JTANK ) .EQ.3) A3D ( JTA NK) =1 . 

27 IF! JHEAD) 29,29,28 

28 IF!NEQBD! JHEAD). EQ.3) A8D « JHEAD) =1 . 

29 SUWV=0. 

DO 99 NSL =1, NSL ICE 
SUMZG=0. 

SUMZN = 0. 

SUMVZ=0. 

RANGE! NSL ) = Z ( NSL+ 1 > - Z ( NSL) 

J B = JHFAD 

IF! NSL .GE .KNUCKL ) JB = JTANK 

DO 85 KZ = 1, NGAUSZ 

ZRC = RANGE! NSL )* VI l K Z ) + Z l NSL) 

RR=S3RT! ( DBO! JB)-CB3! J3 )*( ZRC-ZOBD! JB) ) **2 ) / ABD! JB) ) 

SUMRG=0. 

SUM RN = 0. 

DO 75 KR=1, NGAUSR 
RPR IM E = RR*VR(KR) 

SUMPG = 0. 

SUMPN = 0. 

DO 70 KP = 1, NGAUSP 
PHRIME = P I* VP ( KP ) 

11=11+1 

XRC = RPR IME * COS(PHRIME) +XOBD!JB) 

YRC = RPRIME * S IN ( PHR I ME ) 

RRC = SQR T ( XRC * XRC + YRC ♦YRC ) 

PHIRC = A TAN 2 (YRC, XRC) 

C COMPUTE DOSFS AT ( RRC , ZRC , PHI RC ) REGULAR Q AD CALC. HERE TO 65 

NSET! 1 5 ) =0 
I ERR 3 =0 


QADC0530 

QA0C0540 
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DNEUT=0.0 
DGAMT=0.0 
UDGAMT =0.0 
00 30 J=1,NRGY 
0GAM( J )=0.0 
30 UDGAM(J)=0.0 
00 60 N=1,NS0 
SINS = SIN(PHISO(N)> 

00SS=C0S( PHI SO (N ) ) 

00 60 L = 1 , L SO 
XS0=RS0(L)*C0SS 
YS0=RS0(L)*SINS 
FLN = FL (L ) * FN { N ) 

00 60 M=1 , M SO 
NSETI 1 5 ) =N SE T ( 15H-1 

1 F( NSETI 15 )-NSET( 13) ) 36,34,33 
33 IFC NSET( 15)-NSET( 14) ) 35,35,36 
3 A WRITE (6,101(10(11,1=1,12) 

35 WRITE (6,7)II,RRC,ZR:,P-IIRC,L,M,N 
CALL SLITE (1) 

36 FLMN=FLN*FM(M> 

CALL L ENGTH( SREG(l) ,SSO,XSO, YSO,ZSO( M) ,XRC ,YRC »ZRC , I ERR 1 ) 

I F ( IERR1IA0, AO, 37 

37 WRITE (6, 11)II,RRC,ZRC,PHIRC,L,M,N,IERR1 
N ERR=N ERR 1 

I Ft 2-N ERR ) 20, AO, AO 

AO CALL KERNEL(SREG(1),SS0,ANEUT,AGAM(1) ,UAGAM<1) ,IERR2) 

I ERR3=MAX0( I ERR 3, I ERR 2) 

A=FL3N/SSO**2 

DNEUT = ONEUT + A*ANEUT 

DO A 5 J = 1 , NR GY 

OGAM( J )=DGAM (J ) +A*AGAM( J ) 

A 5 UDGAH( J ) = UDGAM (J)+A*UAGAM(J) 

60 CONTINUE 

ON EUT=DNEUT/12. 5663706 
00 65 J = 1 ,NR GY 

A=GAMEN( J )*CONV( J) /12. 5663706 
DGAMI J ) = A*DGAM(J) 

UDGAM ( J ) =A*UOGAM( J ) 

UDGAMT =UD GAM T + UDGA M ( J ) 

65 DGAMT = OGAMT«-DGAM( J ) 

IF( NSK IP .GE. 3) WRITE (6,6) 1 1 ,RRC , ZRC , PHI RC ,UDG AMT , DGAMT , DNEJT 
C INTEGRATE OVER PHI 

SUMPG = SUMPG + WP ( KP ) *DSA MT 
70 SUMPN = SUMPN +WPIKP) * ONEUT 
C END 0 H I LOOP. GET FINAL INTEGRALS OVER PHI 

SPGKR ) = 2. *P I* SUMPG 
SPN ( KR ) = 2 . *P I* SUMPN 

IFINSK IP .GE.2) WR I TE ( 6 , 9) SPG ( KR ) , SPN( KR) 

C INTEGRATE OVER R 

F = RPR IM E*WR ( KR ) 

SUMRG = SUMRG + F*SP5(KR) 

75 SUMRN = SUMRN ♦ F*SPN(KR) 

C END R LOOP. GET FINAL VALUES OF INTEGRALS OVER R 

SR G( < Z )= SUMRG*RR 
SRN ( < Z )= SUMRN*RR 


Q A DC 0 5 5 0 
QADC0560 
0 A OC 0 5 7 0 
QA0C0580 
Q A0C0590 
0 ADC0600 
0A0CC610 
0 A0C0620 
0 A0C0630 
0 A0C06A0 
0A0C0650 
0 ADC0660 
0 A0C0670 
0A0C0680 
0A0C0690 
OAOCO 700 
Q A DC 0 7 1 0 


0A0C0750 

QA0C0760 


0ADC0790 
0A0C0800 
0 A OC 0 8 1 0 
0 AOCO820 
0ADC0830 
0A0C08A0 
0ADC0850 
QA0C0860 
OAOCO 8 70 
QADC0880 
Q A0C0890 
0ADC0900 
0A0C0910 
QA0C0920 
0A0C0930 
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IF< NSK IP.GE.2) WR I TE ( 6» 9) SRG (KZ) tSRNIKZ) 

C INTEGRATE OVER Z 

SUMVZ=WZf KZ)*PI*RR*RR«- SUMVZ 
SUWZG = WZ(KZ)*SRG(KZ)«-SUMZS 
85 SUMZN = WZ(KZ)*SRN(KZ)+SUMZN 
C ENO OF Z LOOP. GET FINAL VALUES OF INTEGRALS OVER Z. 

SG(NSL> = RANGE I NSL ) *SUMZS 
SN { NSL )= RANGE(NSL)*SUMZN 
ST ( N SL ) = SG(NSL)*SN(NSL) 

VOL ( N SL I - RANGE (NSL )* SUMVZ 
SUM V=S UMV +VO L( NSL ) 

IF(NSKIP.GE.l) WRITE(6,9) SG ( NSL> , SN( NSL) , ST ( NSL) 

99 CONTINUE 

FLOW=SUMV/BURN 
WR I TE ( 6# 2 ) FLOW, SUMV 
FF=l./a3W 
DO 101 1 = 1, NSL ICE 
101 DTIMEl I) =VOL ( I )*FF 

WRITE (6, 10) I ID II), I *1,12) 

WR I TE ( 6,5) 

WRITE( 6,6) ( l,Z( I), RANGE (I) ,VOL(I) ,DTIME< I) ,SG(I) ,SN( I ) ,ST (I) , 1 = 1, 
XNSL ICE ) 

WRITE( 6,6) JJ,Z(JJ) 

IF( JWALL.EQ.l) CALL QADW 
CALL FINAL( 10, NSL ICE) 

999 RETURN 
END 
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$ I BFTC OWQAO DECK 


SUBROUTINE QADW 
C CALCULATE WALL HEATING 

DIMENS ION WGIIOOI, WN(IOO), WT(IOO), WOLUOOI, WCNVI30) 
COMMON /3ADATA/ 


1 

NSFT 

! 16) , 

R SO 

(21) 


ZSO 

(21) 

t 

PH IS 0 

(21 ) 


2 

NBNDZN( 50 ) , 

n:mpzn(50) 


LBO 

(6,53) 

f 

NTRYZ N ( 6 , 50 ) 


3 

NEQBD 

(501 , 

A83 

(50) 


BBD 

(50) 

9 

CBO 

(50) 


4 

XO0O 

(50) , 

YOBD 

(50) 


ZOBD 

(50) 

9 

DBD 

(50 ) 


5 

HYQRAT(30) , 

XSECN 

(30) 


XSECG 

(20,30) 

t 

3110 0 

(30 ) 


6 

BIL01 

(30) , 

BILD2 

(30) 


BILD3 

(30) 

t 

GAMEN 

(30) 


7 

CONV 

(30) , 

COMP 

(20, AO) 


FBAR 

(30) 

♦ 

MATZ 

(20) 


8 

FL 

(21) f 

FM 

(21) 


FN (21) 

, 




9 

ASO 

, X I i SO 

t XI2S0 

, ETA1S0 

, ETA2S0, ALF1 

, ALF2 , ALF3 

X 

AL FA 

, EPSLN 

, FUDGE 

, NBLD 








COMMON 

/GAUS SO/ 

VI 4, 1 6) , W ! 4 , 

16) 






COMMON 

/QADHO/ 










1 

Z 

(100) , 

vz 

(16) 

t 

WZ 

(16) 

t 

VR 

(16) 


2 

WR 

(16) , 

VP 

(16) 

t 

WP 

(16) 

♦ 

VOL 

( 100 ) 


3 

DT IM E 

( 100 ) , 

RANGE 

( 100 ) 

t 

VV 

( 100 ) 

* 




A 

SR G 

(16) , 

SRN 

(161 

♦ 

SPG 

(16) 

9 

SPN 

(16) 


5 

SG 

( 100 ) , 

SN 

( 100 ) 

» 

ST 

( 100 ) 

t 

QNDOT 

( 100 ) 


6 

QG03T 

( 100 ) f 

QTOOT 

( 100 ) 

t 

TGO0T 

( 100 ) 

, 

TNDOT 

( 100 ) 


7 

TTODT 

( 100 ) , 

OQG 

( 100 ) 

» 

DQN 

( 100 ) 


DOT 

t 100 ) 


8 

OT G 

( 100 ) , 

OTN 

( 100 ) 

? 

DTT 

( 100 ) 

9 

QG 

( 100 ) 


9 

QN 

( 100 ) , 

OT 

( 100 ) 

» 

TG 

( 100 ) 

9 

TN 

( 100 ) 


X 

TT 

( 100 ) , 

T 

( 100 ) 


JTANK, 

JHEAD, 


KNUCKL 

, NS LI CE 

COMMON 

/O ADM AN / 










1 

SR EG 

( 20 ) , 

AGAM 

(30) 

, 

OGAM 

(30) 

9 

ID 

( 12 ) 

9 

2 

UOGAM 

(30) , 

UAGAM 

(30) 









EQUIVALENCE! LSO,NSETI 1) I , <MSO,NS£T!2>) , ( NSO, NS ET ( 3 ) ) , 

1! NRGY.NSET! 7) ) , ( I SRC ,NSET< 11) } , ! NRCVR , NSET (9 >> , C INFJT , NSET ( 12 I ) 

2 P0RMATQ6H0 WALL VOLUME = 1PE12.4,5H CM3) 

3 FORMAT ( 20H1 W-HEATING CONTROL , 6 16 , l PE12 . 4 > 

4 FORMAT ( 8E9 .4 ) 

5 FORMAT! IH0,68X,27H WALL HEATING RATES CAL/SEC/7 H I,I3X,76HZ!I 

1) ZH + D-ZU) VOL ! I ) OTIME(I) GAMMA NEUT 

2 TOTAL) 

6 FORMAT ( 1H 1 6 , 5 X IP 7E 12. 4) 

7 FORMAT ! 7H0 RCVRI3,4H R = 1PE12.4,4H Z=E12.4,6H PH 1= E 1 2 .4 , 1 1H S0QADC0370 

1URCE PT3I3 ) QADC0380 

8 FORMAT ( 4! 15, 14 >,E9.4) 

9 FORMAT! IP 3E12.4) 

10 FORMAT ( 1H1, 5X, 12A6) QADC0390 

11 FORMAT < 7H0 RCVR13,43 R=1PE12.4,4H Z=E12.4,5H PHI= E12.4, 1 IH S00ADC0400 

1URCE PT3I3.8H TYPEI2.29H GEOMETRY ERROR AT THIS POINT) QA0C0410 

PI=3. 141593 
INEUT = P 
2 C 11=0 
N ERR = 0 

READ! 5,8) NP HI , NR , NZ , NOUMM Y ,N SKI P , NOUM, J WHED , JWTNK , ZPR I ME 
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WRITE(6,3) NPHI ,NR,NZ,NSKI P, JWHED, JWTNK , ZPRIME 
READ(5,4) (WCNV(I ),I = 1,NRGY) 

WR ITE C 6, 12 )( WCNV( I ) ,I=1,NRGY> 

12 FORMAT ( 14HL WALL CNV = / ( IP 1 OE 12 . 4) ) 

NGAUSR =2**NR 
DO 251 1=1, NGAUSR 
VR( I)=V(NR,I ) 

251 WR ( I ) = W< NR , I ) 

NGAUSZ =2**NZ 

DO 252 1 = 1, NGAUSZ 

vzc d=v(nz,i ) 

252 WZ( I ) =W(N Z , I ) 

IFCNPHI) 253,253,254 

253 NPHI=0 
NGAUSP =1 
VP( 1)=0. 

WP( 11*1. 

GO TO 256 

254 NGAUSP = 2**NPHI 

DO 255 1=1, NGAUSP 

VP C I > = Vf NPHI, II 

255 WP( I) = W ( NP HI , I ) 

256 I Ft JWHED+JWTNK ) 999,999,259 
259 I FI JWTNK) 27,27,26 

26 IF(NEQBDt JWTNK). EQ. 3) ABD( JWTNK) =1 . 

27 IF(JWHED) 29,29,28 

28 IFCNEQBDI JWHED). EQ. 3) ABD I JWHED) =1 . 

29 SUM V=0 • 

DO 99 NSL =1, NSL ICE 
SUMZG=0. 

SUMZN = 0. 

SUM VZ =0. 

RANGE ( NSL )= Z t N SL + 1 ) - Z ( NSL ) 

IF! NSL-KNUCKL) 295,296,296 

295 J B=J HEAD 
JC =J W HED 
GO TO 297 

296 J B=JTANK 
JC=JWTNK 

297 DO 85 KZ=1, NGAUSZ 

ZRC = RANGE(NSL)*VZ(KZ)+Z (NSL) 

RR=SORT( ( DBD( JB )-C BD ( JB ) * ( ZRC-Z OBD < JB) ) **2 ) /ABD( JB) ) 
RW=SQRT( ( DBD( JC )-CBO ( JC ) #( ZRC-Z OBD ( JC ) ) **2 > / ABD ( JC ) ) 
SUMRG=0. 

SUMRN = 0. 

DO 75 KR=1, NGAUSR 

RPR IM E = RR +VR (KR)*(RW-RR) 

SUMPG = 0. 

SUMPN = 0. 

DO 70 KP = 1, NGAUSP 
PHRIME = PI*VP(KP) 

11 = 11+1 

XRC = RPR IME * COS(PHRIME) «-XOBD(JB) 

YRC = RPR IME * SIN ( PHR I ME ) 

RRC = SQRT( XRC*XRC + YRC*YRC) 

PHIRC = A TAN 2 (YRC, XRC) 



C COMPUTE OOSFS AT ( RRC , ZRC , PH I RC ) REGULAR QAO CALC* HERE TO 65 

NSFT( 1 5 ) = 0 
l ERR 3 =0 
ON EUT = 0*0 
DGAMT=0.0 
UDGAM T =0* 0 
DO 30 J = 1 » NR GY 
DGAM{ J )=0 .0 
30 UDGAM { J 1 = 0*0 
DO 60 N=1 1 N SO 
S IN S = S IN( PHI SO (N) ) 

COSS=COS( PHI SO (N) > 

00 60 L = 1 ? L SO 
X SO=R SO ( L )*COSS 
YSO*RSO(L )*SINS 
FLN=FL(L)*FN(N) 

00 60 M = i , M SO 

NS ET( 1 5 ) =N SE T ( 15) + 1 

IF( NiS ET| 1 5 ) - N S E T ( 13) ) 36, 34 , 33 

33 IF(NSET<15)-NSET<14))35 # 35,36 

34 WRITE (6, 10)( ID< I ) ,1=1,12) 

35 WRITE (6,7)11, R RC »ZRC, PH IRC, L,M,N 
CALL SL ITE (1) 

36 FLMN=FLN*FM( M ) 

CALL L ENGTH ( SR EG ( 1 ) , S S 0, X SO , Y SO , ZSO ( M) , XRC , YRC , ZRC , I ER R 1 ) 

IF< IERR1)40, 40,37 

37 WRITE 16, 11)11, RRC , ZRC, PH I RC , L , M, N, I FRR1 
N ERR=N ERR +1 

I F ( 2-NERR ) 20,40,40 

40 CALL KERNEL! SREGI 1),SS0,ANEUT ,AGAM(l) ,UAGAM<1) , IERR2) 

1 ERR 3=MAX 0 ( I ERR 3, IERR2) 

A=FLM N /SSO** 2 

DNEUT = DNEUT+A*ANEUT 

DO 45 J =1 , NRGY 

DGAM{ J )-DGAM ( J )+A*AGAM( J) 

45 UDGAM { J ) = UDGAM ( J ) +A* UA GA M ( J) 

60 CONTINUE 

DN EUT = DNEUT/ 12. 5663706 
DO 65 J =1 ,NR GY 

A=GAM EN( J )*WCN V( J ) / 12 . 5663 706 
OGAM! J )=A*DGAM( J ) 

UDGAM ( J ) =A*UDGAM< J ) 

UDGAM T =UD GAM T + UDGA M( J ) 

65 DGAMT = DGAMT +OGAM I J ) 

IF( NSK IP.GE.3 ) WRITE (6,6) I I ,RRC , ZRC, PH I RC , UDG AMT , DGA MT , DNEUT 
C INTEGRATE OVER PHI 

SUMPG = SUMPG + WP( KP) *05AMT 
70 SUMPN = SUMPN +WPIKP) * DNEUT 
C END PHI LOOP. GET FINAL INTEGRALS OVER PHI 

SPGUR ) = 2**PI*SUMP3 
SPN ( K R ) = 2 **P I*$UMPN 

IFINSK IP.GE.2) WR I TE ( 6 , 9) SPG ( KR ) , S PN ( KR) 

C INTEGRATE OVER R 

F = RPRIME*WR(KR) 

SUMRG = SUMRG ♦ F*SP3(KR) 

75 SUMRN = SUMRN + F*SPN(KR) 


Q ADC 0 530 
0ADC0540 
0 A DC 0 5 5 0 
QA0CC560 
Q A DC 0 5 7 0 
0 A DC 0 5 8 0 
0 A DC 05 90 
Q A0C0600 
Q A DC 0 6 1 0 
0 A DC 0 6 2 0 
0 ADCC630 
Q A DC 0640 
QADC0650 
OA DC 06 60 
0ADC0670 
QADC0680 
OA DC 06 90 
0ADC0700 
0 A DC 0 7 1 0 


Q ADC0750 
OADC0760 


0 ADC0790 
QADCC800 
Q A DC 0 8 1 0 
0 A DC 0 8 2 0 
Q A DC 08 30 
QA0C0840 
0 A DC 0 8 5 0 
0 A DC 0 8 6 0 
Q ADC 08 70 
0 A DC 0 8 8 0 

OADC 0900 
0 A DC 0 9 1 0 
OADC 09 20 
0ADC0930 
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C END R LOOP. GET FINAL VALUES OF INTEGRALS OVER R 

SRGIKZ ) = SUMRG* (RW-RR ) 

SRN I K Z ) = SUMRN* ( RW-RR ) 

IF! NSK IP.GE.2) WR I TE ( 6 » 9) SRG ( KZ ) , SRN( KZ ) 

C INTEGRATE OVER Z 

SUM VZ = WZ ( KZ ) *P I*(RW+RR)*!RW-RR) + SUMVZ 
SUMZG = WZIKZ)*SRGIKZ)*SUMZ3 
85 SUMZN = WZIKZ)* SR NIKZ)-* SUMZN 
C END DF Z LOOP. GET FINAL VALUES OF INTEGRALS OVER Z. 

WG( NSL ) = RANGE INSL)*SUMZS 
WNINSL) = RANGE INSL )*SUMZN 
WT( NSL )- WGINSD+WNINSL) 

WOL ( N SL ) = RANGE INSL )* SUM VZ 
SUMV=SUMV+WOLI NSL ) 

IFINSKIP.GE.il WRITE!6,9) WG ( NSL I , WNI NSL) , WT ( NSL I 
99 CONTINUE 
J J=NSL ICE ♦ 1 
WRITE! 6» 2 ) SUM V 
WRITE(6,10)I IDU 1,1=1,12) 

WRITE! 6,5) 

WRITE! 6,6) ( I,Z! I ), RANGE! I ) ,WOL!I) , D TI ME ! I ) , WG ! I ) ,WN II) , WT ! I ) , 1= 1, 
XNSL ICE ) 

WR ITE ! 6, 6 ) J J, Z ( J J ) 

DO 101 1 = 1, NSL ICE 
SGI 11 = SGI I ) + WG(I ) 

SNI I ) = SN( I ) * WN l I ) 

101 ST! I ) = ST! I )+WT ! I ) 

IF! ZPR IME.NE. 0.0) CALL QAOB 
1 |NPHI,NR,NZ,NSKIP,JWHED,JWTNK,WCNV,ZPRIME) 

999 RETURN 
END 
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$ I BFTC QBQAD DECK 


SUBROUTIN E QADB 

1 (NPHI.NT ,NZ,N SKIP, JWHED.JWTNK, WCNV.Z PRIME) 
C CALCULATE HEATING FROM TANK BOTTOM 


DIMENSION WG ( 100)* WN(IOO) 

, WT(IOO) , 

WOL <100 ) 

, WCNV (30 ) 

COMMDN 

/3 ADA TA / 









1 

NSET 

( 16) 

» 

R SO 

(21) 

9 

ZSO 

(21 ) 

♦ 

PH ISO 

(21 ) 

2 

N BNDZN ( SO) 

t 

NO MPZ N ( 50 ) 

» 

L80 

(6,50 ) 


NTRYZN<6 ,50 ) 

3 

NEQBD 

( 50 ) 

9 

A3D 

(50) 

r 

BBD 

(50) 


CBO 

(50 > 

4 

X OBO 

(50) 

t 

Y0B3 

( 50) 

» 

ZOBD 

(50) 

9 

OBO 

(50) 

5 

H Y DR A T ( 30) 

t 

XSECN 

(30) 

» 

XSECG 

(20,30) 

* 

BILOO 

(30 ) 

6 

BILD1 

( 30) 

t 

B I LD2 

( 30) 

9 

BILD3 

(30) 

♦ 

GAMEN 

(30 ) 

7 

C3NV 

(30) 

9 

COMP 

(20, 

40) , 

FBAR 

(30) 

» 

MATZ 

(20 ) 

8 

FL 

(2 1) 

9 

FM 

(21) 

t 

FN (21) 

t 




9 

AS3 

, XI ISO 

, XI2S0 

, ETA1S0 

, ETA2S0, ALF1 

, ALF2 , 

ALF3 

X 

AL F4 

, EPSLN 

, FUDGE 

, NBLD 








C0MM3N 

/G A US SO/ 

V ( 4, 1 6 ) , W ( 4 , 

16) 






COMMDN 

/3 ADHD/ 










1 

Z 

(100) , 

VZ 

(16) 

9 

WZ 

(16) 

? 

VR 

( 16 ) 

t 

2 

WR 

(16) , 

VP 

(16) 

9 

WP 

(16) 


VOL 

(100) 

t 

3 

DTIMF 

(100) , 

RANGE 

(100) 

t 

vv 

(100) 

9 




4 

SR G 

(16) , 

SRN 

(16) 

t 

SPG 

( 16) 

f 

SPN 

( 16 ) 

t 

5 

SO 

(100) , 

SN 

( 100) 

t 

ST 

(100) 

9 

QNDOT 

(100 ) 

* 

6 

QGDOT 

(100) , 

QTDQT 

(100) 

» 

TGDOT 

(100) 

9 

TNDQT 

(100) 

9 

7 

TTD3T 

(100) , 

DOG 

(100) 


DQN 

(100) 

t 

DOT 

( 100) 

9 

8 

DTG 

(100) , 

OTN 

(100) 

t 

DTT 

(100) 

» 

OG 

( 100) 

9 

9 

ON 

(100) , 

OT 

(100) 

♦ 

TG 

(100) 

» 

TN 

(100) 

9 

X 

TT 

(100) , 

T 

(100) 

t 

JTANK, 

JHEAD, 


KNUCKL 

, NSLICE 

COMMON 

/Q ADMAN/ 










1 

SR EG 

(20) , 

AGAM 

(30) 

9 

DGAM 

(30) 

t 

ID 

(12) 

9 

2 

UDGAM 

(30) , 

UAGAM 

(30) 









EQUIVALENCE! L SO, NSET(l) ) , IMSO,NSET<2) ) , ( NSO, NS ET ( 3 ) ) , QADC0170 

II NRGY,NSET(7) ) , < I SRC, NSET! Ill) , ( NRC VR » NSFT ( 9 ) t QADC0I30 

2 FORMAT 1 16HL BOTT VOLUME = 1PE12.A.5H CM3) 

5 FORMAT ( 1H0,68X,27H30TTM HEATING RATES CAL/SEC/7 H I,13X,76HZ(I 

I) Zll + n-ZII) VOL ! I ) DTI ME ( I ) GAMMA NEUT 

2 TOTAL) 

6 FORMAT ( 1H 1 6 , 5 X IP 7E 1 2 . A) 

7 FORMAT I 7H0 RC VR I 3 * AH R=1PE12.A,AH Z=E12.A,SH PHI= E12.A, 1 IH S0QA0C0370 

1URCE PT3I3 ) QAOC03BO 

9 FORMAT! IP 3EI2. A) 

10 FORMAT! IHL ,5X, I2A6) 

11 FORMAT ( 7H0 RCVRI3.AH R=1PE12.A,AH Z=E12.A,6H PH I = E 1 2 . A, 1 1H SOQAOCOAOO 

UJPCE PT3I3.8H TYPEI2,29H GEOMETRY ERROR AT THIS POINT) QAOCOAIO 

Pl=3.lA1593 
20 11=0 
NERR=0 
NGAUSR =2**NR 
00 251 1= 1 »NGA USR 
VR! I ) = V I N R * I ) 

251 WR( I ) =W( NR » I ) 
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NGAUSZ=2**NZ 
DO 252 I = l,NGAUSZ 
VZ( I) = V( M Z » I ) 

252 WZ{ I ) = W( N Z , I ) 

IFINPHI ) 253,2 53,254 

253 NP H I =0 
NGAUSP = 1 
VP( 1 ) =0. 

WP( 1)=1. 

GO TO 29 

254 NGAUSP = 2**NPHI 

DO 255 1=1, NGAUSP 

VP ( I ) = V ( NP H I ,1 J 

255 WP( I) = WCNPHI , I) 

29 SUMV=0. 

NSL =1 
SUMZG=C. 

SUMZN=0. 

SUM VZ = 0. 

RANGE! NSL )= Z( 1I-ZPRIME 
I FI NSL-KNUCKL ) 295,296,2 96 

295 J C =J W HED 
GO TO 297 

296 J C =J WTNK 

297 DO 85 KZ = 1 , NGA USZ 

ZRC = RANGE(NSL)*VZ(KZ)+ZPRIME 

RW= SQR T { ( DBD ( JC )-CBD( JC )* ( ZRC-ZOBD( JO )**2)/ABD(JC) ) 

SUM R G= 0 • 

SUMRN=0. 

DO 75 KR = 1 , N GA USR 
RPR IM E = RW*VR(KR) 

SUM PG = 0. 

SUM PM = .0 • 

DO 70 KP = 1, NGAUSP 
PHR IMF = Pl*VP(KP) 

11=11+1 

XRC = RPR I ME * COS I PHR I ME ) +XOBDUB) 

YRC = RPR IME * S IN ( P HR I ME ) 

RRC = SQRT( XRC* XRC + YRC* YRC ) 

PH IRC = A TAN 2 I YRC,XRC ) 

C COMPUTE DOSES AT I RRC , ZRC , PH I RC ) REGULAR QAD CALC. HERE TO 65 

NISETI 1 5 ) =0 
I ERR 3 =0 
DN EUT = 0.0 
DGAMT = 0.0 
UDGAM T =0 . 0 
DO 30 J = 1 , NR GY 
OGAM( J ) = 0 . 0 

30 UDGAM ( J ) = 0 .0 
DO 60 N = 1 , N SO 

S IN S= S INI PHISO(N) ) 

COSS=COSI PHISO(N) ) 

DO 60 L = 1 , L SO 
XSO=RSO( L )*COSS 
YSO=R SOIL )*SINS 
FL N =FL (L )*FN(N ) 


QADC0530 
0 A DC 0 5 4 0 
QADCC550 
GADC0560 
Q A DC 0 5 7 0 
0ADC0580 
Q ADC0590 
0 ADC0600 
0 A DC 0 6 1 0 
0 A DC 0 6 2 0 
GADC0630 
0 A DC 0640 
0 ADC0650 
0 ADC0660 
0 A DC 06 70 
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** 


00 60 M=1 , M SO 
NSETI 15)=NSETI 15) + i 

I FI NSETI 15)- NSETI 13) ) 36, 34,33 
33 IF( NSETI 15)-NSET( 14) ) 35,35,36 
3 A WRITE <6, 10) ( 101 1 ),I=1 ,12) 

35 WRITE <6,7)II,RRC,ZRC,PHIRC,L,M,N 
CALL SLITF (1) 

36 FLMN=FLN*FM(M) 

CALL LENGTH! SR EG(l) , SSO, XSO, Y SO, ZSO I M) , XRC ,YRC ,ZRC , IERRI ) 
IFI IERRDAO, AO, 37 

37 WRITE 16, 11)11, RRC,ZRC,PHIRC,L,M,N, IERRI 
N ERR=N ERR +1 

1 F( 2-NERR ) 20, AO, AO 

AO CALL KERNEL! SREG(l), SSO, ANEUT.AGAMU) ,UAGAM(1) ,IERR2) 
IERR3=MAX0< IERR3, IERR2) 

A=FLMN/SS0**2 

DNEUT=DNEUT+A*ANEUT 

DO A5 J=1 ,NRGY 

OGAM( J )=0GAM( J)+A*AGAMIJ) 

A 5 UDGAMI J)=UDGAMI J)*A*UAGAM( J) 

60 CONTINUE 

0NEUT = 0NEUT/12. 5663706 
DO 65 J =1 ,NR GY 

A = GANEN«J )*WCNV(J)/12. 5663706 
DGAMI J )=A*DGAM< J) 

UDGAM ( J ) =A*UDGAM ( J ) 

UDGANT =UDGAMT +UOGAMI J ) 

65 DGAMT =DGAMT+OGAM( J ) 

IFINSK IP. GE. 3) WRITE I 6, 6) II ,RRC,ZRC,PHI RC , UOG AMT , DGAMT , DNEUT 
C INTEGRATE OVER PHI 

SUMPG = SUMPG ♦ WP(KP)*DGAMT 
70 SUMPN = SUMPN ♦ WP I KP ) * ONEUT 
C END PHI LOOP. GET FINAL INTEGRALS OVER PHI 

SPGKR) = 2.*P I*SUMPG 
SPNKR ) = 2.*P I*SUMPN 

IFINSK IP. GE. 2) WRITEI6.9) SPG I KR) , SPNI KR) 

C INTEGRATE OVER R 

F = RPR IM E*WR I KR ) 

SUMRG = SUMRG + F* SPG ( KR ) 

75 SUMRN = SUMRN ♦ F*SPNIKR) 

C END R L03P. GET FINAL VALUES OF INTEGRALS OVER R 

SR Gl < Z )= SUMRG*RW 
SRNKZ >= SUMRN*RW 

IFINSK IP. GE. 2) WRITEI6.9) SRG I KZ) , SRNI KZ) 

C INTEGRATE OVER Z 

SUMVZ=WZ(KZ)*PI*RW*RW*SUMVZ 
SUMZG = WZCKZ)*SRG(KZ ) +SUMZS 
85 SUMZN = WZ(KZ)*SRN(KZ)+SUMZN 
C END DF Z L03P. GET FINAL VALUES OF INTEGRALS OVER Z. 

WGI NSL )= RANGE (NSL )*SUMZ3 
WN ( NSL )= RANGE I NSL )* SUMZN 
WTI NSL )= WGINSL M-WNINSL) 

WOL l N SL ) = RANGEINSL )* SUMVZ 
SUM V = SUM V + WO LI NSL ) 

IFINSK IP. GE. 1) WR I TE C 6 ,9) WG ( NSL ) , WN I NSL ) , WT ( NS L ) 

99 CONTINUE 


0ADC0680 

0A0C0690 

0A0CC700 

0A0C0710 


0 A DC 0 7 5 0 
0A0C0760 


Q A DC 07 90 
0ADC0800 
Q A DC 0 8 1 0 
0ADC0820 
Q A DC 08 30 
0A0C08A0 
QADCC850 
0 ADC0860 
QADC0870 
0 A DC 0 8 0 0 

0ADC0900 
QADC09 10 
QADC0920 
0ADC0930 
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1 = 1 

WR l TE ( 6» 2 ) SUM V 

WR ITE( 6, 10)1 ID< l), 1=1, 12) 

WR I TE I 6» 5 ) 

WRITE (6, 6) I » Z ( I ) » RANGE 41) » WOL (I )*0TIME(II«W6(I)*WN(I)fWT(I) 
SGI 1 ) = SGI lUWG(l) 

SN ( 1 ) = SN < 1)+WN(1> 

ST C 1 > = ST ( 1 ) ♦WT ( 1 ) 

<599 RETURN 
END 



SIBFTC INQAD 


DECK 


SUBROUTINE INPUT(ID) 

COMMON / 3 A DA T A / 

1 NSET (16) , R SO (21) , ZSO (21) , PHISO (21) , 

2 NBNDZN( 50) , NCMPZN(50) , LBD (6,53) , NTRYZN(6,50) , 

3 NEQBD (50) , ABD (50) , BBD (50) , CBD (50) , 

4 XOBD (50) , YOBD (50) , ZOBD (50) , DBD (50) , 

5 HYDRAK30) , XSECN (30) , XSECG (20,33), BILOO (30) , 

6 BILD1 (30) , BILD2 (30) , BILD3 (30) , GAMEN (30) , 

7 CONV (30) , COMP (20,40), EBAR (30) , MATZ (20) , 

8 FL (21) , FM (21) , FN (21) , 

9 ASO , XliSO , XI2S0 , ETA1S0, ETA2S0, ALF1 , ALF2 , ALF3 , 

X ALF4 , EPSLN , FUDGE , NBLD 

COMMON /STDATA/ 

1 NAME (100) , ENR3 (20) , XSECO ( 20 , 100 ) , XS ECE A( 20 ) , 

2 C (4,6,9) , EMI N3B ( 8) , NA MB ( 8 ) , XSECNT(IOO) ,DCONST 

DIMENSION ID ( 12) , NAM(20) 

EOU I VALENCE ( L SO, NSET ( 1) ) , ( MSO , NSET ( 2 ) ) , ( NS0,NSET(3) ) , 

1( MAT, N SET ( 4) ) , ( NCOMP , NSE T ( 5) ) , ( NREG , NSET ( 6 ) ) , ( NRGY.NSET (7) ) , 
2(NB0UND,NSET(8)),( I SRC ,NSE Till)) , ( I NEUT, NSET ( 12 ) ) 

4 FORMAT ( 33H1 NO SOURCE. RUN TERMINATED.) 

5 FORMAT (12A6) 

6 FORMAT ( 8{ 15, 14) ) 

7 FORMAT ( 8E9.4 ) 

8 FORMAT ( 15, I4,7E9.4) 

9 FORMAT ( 6H1 12A6/13H0 C0NTR0L16I6I 

11 FORMAT ( 18H0 ZONE BN3S C0MP6(12H BNO PK )//) 

12 FORMAT (1516) 

13 FORMAT ( 27H0 BND EQ PARAMETERS...//) 

14 FORMAT (21 6, 1P7E12.4) 

15 FORMAT ( 93 HO GRP EPS/E DELTA/K Al/BO A2/B 

11 A3/B2 A4/B3 EBAR//) 

151 FORMAT ( 5H1QADX, 10X, A6, 1 6H BUILDUP IS USED) 

16 FORMAT ( 17, 5X1P9E12.4) 

17 FORMAT ( 10HL COMP/GRP 9( I 9 , 3X) ) 

18 FORM AT ( 1H0, 1 1X»9(5X, I 3,2X,A2) ) 

100 FORMAT ( 12H0 SOURCE 1P8E 12.4) 

101 FORMAT ( 12H0 R 1P8E 1 2 . 4/ ( 12 X8E1 2 . 4 ) ) 

102 FORMAT ( 12H0 Z IP 8E 12 . 4/ ( 12 X8E1 2 . 4 ) ) 

103 FORM AT ( 12 HO PH I 1P8E 1 2. 4/ ( 12 X8E1 2 . 4 ) ) 

104 FORMAT ( 12H0 F ( L ) 1P8E 12 . 4/ ( 12 X8E1 2 . 4 ) ) 

105 FORMAT ( 12H0 F ( M ) IP 8E 12. 4/ ( 12 X8E1 2 . 4 ) ) 

106 FORMAT ( 12H0 F ( N ) IP 8E 12. 4/ ( 12 X8E1 2 . 4 ) ) 

Z ER0 = 0 .0 

MON E=- 1 

READ (5 ,5) ( ID( I ) ,1 = 1 ,12) 

READ (5 ,6) INSET! I ) ,1 =1,14) 

WRITE (6, 9)( ID ( I ),l=l ,12) , (NSET (I) ,1=1 ,14) 

IF( ISRC)19, 19, 20 

19 IF( ASO )191,191,30 

191 WRITE (6»9)( ID C I ),I = 1,12) ,(NSET(I) ,1=1 ,14) 
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WRITE (6,4) 

CALL EXIT 
20 LU=L SO +1 
MU=MS0+1 
NU=N SO + 1 

READ (5 ,7)AS0, XIISO, XI 2S0,ETA1 SO, ETA2 SO 
WRITE (6, 100 )ASO, XIISO, XI2S0.ETA1 SO, ETA2S0 
READ (5 ,7) (RS0(L),L=1»LU) 

READ (5 , 7) ( ZSO(M) , M=1 * MU) 

READ (5 ,7) (PHISO(N) ,N=1 ,NU) 

WRITE (6, 101 HR SOIL ) » L = 1 ,LU) 

WRITE (6, 102 )( ZSO(M) , M=1 ,MU) 

WRITE (6, 103 )( PHISO(N) ,N = 1 ,NU) 

GO TO ( 30,22) , I SRC 
22 READ (5 , 7 ) (F L (L ) , L= 1 ,LU) 

READ (5 »7)(FM(M),M=1»MU) 

READ (5 , 7)(FN(N),N=1»NU) 

WRITE (6, 104) (FL(L) ,L=1,LU) 

WRITE (6, 105 HFM(M) ,M=1,MU) 

WRITE (6, 106)(FN(N),N=l,NU) 

30 WRITE (6, 11) 

DO 32 J = 1 » NR EG 

READ (5 ,6) IU.NCMPZN (J) , (LBD(I ,J) ,NTRYZN(I ,J) ,1=1, IU) 

NBNDZN (J ) = IU 

32 WRITE (6, 12)J»IU*NCMPZN(J) , (LBD (I , J) , NTRYZNt I , J) , 1= l, IU ) 

READ (5 , 8 ) ( K , NE Q30 ( K) , ABD ( K) ,BBD ( K) ,CBD ( K ) , XO BD ( K) , YOBD(K), 

1Z0BD1K ),DBD( K) ,I=1,NB0UND) 

WRITE (6,13) 

WRITE (6, 14) (I,NEQB3(I ) ,ABD(I ) ,BBD( I) ,CBD( I) ,XOBD( I ), YO BD( I ) , Z OB 
1D( I ) , DBD( I ), I=1,N30UND) 

READ (5,6) NBLD,(MATZ(I ) ,1=1, MAT) 

DO 45 J=1 , NCOMP 

45 READ (5 ,7) (COMP( I, J), 1=1, MAT) 

DO 41 1=1, MAT 
L=MATZ ( l ) 

NAM ( I ) =NAME ( L ) 

HYDRAT { I ) = 0.0 
IF(L.EQ.l) HYDRA T( I ) = 9.0 
41 XSECN ( I ) = XSECN T ( L ) 

READ (5,7) (EBARJI ) , I =1 ,NRGY) 

READ ( 5, 7)(GAMEN(I),I=1,NRGY) 

AL F2 = • 29 
AL F3= *83 
AL F4= « 58 
EPSLN = l % E-6 
FUDGE* l.E-3 
CALL XSECQD 

WRITE (6,151) NAM3 ( N3L0 ) 

WRITE (6, 16 )ZERO,EPSLN, FUDGE , ALF1 ,ALF2 , ALF3 , ALF4 
WRITE (6,15) 

DO 65 1=1 , NR GY 

65 WRITE (6, 16) I,GAMEN( I ) ,C ONV ( I ) ,BI LDO ( I ) ,BI LD1 ( I ) ,BILD2 ( I) , BILD31 I) 
1 , EBAR ( I ) 

DO 70 K=1 »MA T» 9 
IU=MIN0(K+8,MAT) 

WRITE (6, 17) ( I ,I=K,IU) 
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00 69 J=1 ,NCOMP 

69 WRITE 16, 16) J, ICOMPCI , J) ,I=K,1 
WRITE (6,18) <MATZ( I ) ,NAM(I ) 


00 70 J=1 ,NR. GY 

70 WRITE (6, 16) J, ( XSECG ( I , J ) , I =K ,] 
CALL WEIGHT 
RETURN 



$ I BFTC WTQAD DECK 


SUBROUTINE WEIGHT 
COMMON /QADATA/ 

1 NSET (16) , R SO (21) , ZSO (21) , PHISO (21) , 

2 NBNDZN(50) . NCMPZN(50) , LBD (6,50) , NTRYZN(6,50) , 

3 NEQBD (50) , ABO (50) , BBD (50) , CBD (50) , 

4 XOBD (50) , Y030 (50) , ZOBD (50) , DBD (50) , 

5 HYDR AT( 30 ) , XSECN (30) , XSECG (20,30), BILDO (30) , 

6 BIL01 (30) , BIL02 (30) , BILD3 (30) , GAMEN (30) , 

7 CON V (30) , COMP (20,40), EBAR (30) , MATZ (20) 

EQUIVALENCE (LS0»NSET(1)) , (MS0,NSET(2) ) , ( NSO, NSET( 3 ) ) , 

1( MAT, NSET (4) ) , ( NCOMP , NSE T ( 5 ) ) , ( NREG , NSET ( 6 ) ) , ( NRGY , NS ET ( 7 ) ) , 
2(NB0UND,NSET(8)),(ISR:,NSET(11) ),(INEUT,NSET(12)) 

DIMENSION RSQ( 6) , X(6), Y(6), Z(6), NQB(6) , CFLAG16 ) , QFLAGI6) , 

1 XFLAG(6), YFLAGI6), ZFLAS(6), V(16) , W(16) 

DATA! V( I) , 1 = 1, 16) /O. 529953E-2 ,0. 277125E-1 ,0 .67 1844E- 1 ,0 . 122298, 
XO. 191062, .270992,0.359198,0.452494,0.547506,0.640802,0.729008, 

YO. 808938, .877702,0.932816,0. 972288,0.9947/ 

DATA! W ( I ) , 1= 1,1 6) /O. 01 3 5 762 ,0.03112 68,0.0475 793,0.06 23145, 

XO. 0747 980, 0.0 845783, 0.0913 01 7, 2*0. 0947253,0.0913017,0.0845783, 

YO. 0747980, 0.0623145,0. 0475793 , 0.031 1268 , 0. 01 35762/ 

DATA PI/3.1415927/ 

500 FORMAT ( 28H1 J VOL WEIGHT) 

501 FORMAT( 14, 1P2E12.4) 

502 FORMA T ( 4H0SUM IP 2E 12.4) 

WR ITF( 6,500) 

VTOTAL =0.0 
WTOTAL =0.0 
DO 99 J=1 ,NREG 
IF( NBNDZN ( J ) ) 99,99,10 
10 IU= NBNDZN ( J ) 

VOL =0.0 
NX = 0 
NY = 0 
NZ =0 
NC = 0 
N0 = 0 

DO 20 1=1, IU 
KK=LBD( I, J ) 

K =1 ABS (KK ) 

FL AG=KK/< 

NEQ=N EQBD ( K ) 

GO T0(99, 12, 13,14,15,16) ,NEQ 

12 NQ=NQ+ 1 
NQB( NQ )=K 
QFLAGI NQ ) =FLAG 
GO TO 20 

13 NC=NC+ 1 

RSQ ( NC )=DBD( K ) 

CFL AG( NC ) =FLAG 
GO TO 20 
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14 NX=NX *■ 1 

X ( NX ) =DBD ( K ) 

XFL AG( NX ) =FLAG 
GO TO 20 

15 NY =NY + 1 

Y( N Y ) =D8D ( K ) 

YFLAG(NY)=FLAG 
GO TO 20 

16 NZ =NZ + 1 

Z ( NZ 1 = DBD ( K ) 

ZFLAGINZ )=FLAG 
20 CONTINUE 

C CALCULATE VOLUMES 

I F( N Z- 2 > 99, 26,22 

22 ZMAX=+1.0E+8 
ZM IN=-7.MA X 
00 25 1=1, NZ 

IF( ZFLAG1 I >123,99,24 

23 ZMIN=AMAXHZMIN,zm ) 

GO TO 25 

24 ZMAX=AMIN1(ZMAX,Z(I) ) 

25 CONTINUE 

OZ = ZM A X— ZM IN 
GO TO 261 

26 DZ = ZFLAGI 1 )*ZU) ♦ ZFLAG (2 ) *Z (2) 

261 I Ft (NX.NE.2) .OR. (NY.NE.2) > GO TO 27 

OX = XFLAG(l)*xm ♦ XFLAG (2 ) *X (2 ) 

OY = YFLAGI 1 I* Yf 1 ) * YF LAG ( 2 ) * Y ( 2 ) 

C RECTANGULAR PARALLELP I PE!) 

VOL =DX*DY*DZ 

27 I F( NC ) 99,32,28 

28 00 30 1=1, NC 

C CIRCULAR CYLINDER 

VV=CFLAG< 1 1*PI*RSQ< I )*DZ 
30 VOL =VQL+VV 
32 I Ft NO >99, 75,35 
35 NQQ=NQ 

DO 50 1=1, NQQ 
NQ=NQB(L> 

VV=0.0 

00 4 5 1 = 1,16 

C CALCULATE VOLUME OF FULLY ROTATED QUAORATI C SURFACE (TYPE 2>, 

C BOUN0ED BY Z(l> AND Z 1 2 > , BY 16 POINT GAUSS QUADRATURE. 

zz=z( i >m i >*iz( 2 )-z( i) > 

RR=(DBO(NQ)-CRO(NQ>*(ZZ-ZOBO(NQ) )**2> /ABD(NQ) 

C CROSS-SECTION IS ASSUMED CIRCULAR IN Z PLANE 

45 VV =VV «-W(I>*PI*RR 
50 VOL =VOL+ QFLAG(L )*VV*OZ 
75 WT=0.0 

C CALCULATE WEIGHT 

K =NCMP ZN ( J ) 

DO 80 1 = 1, MAT 
80 WT=WT+V3L*C0MP (I,K) 

WRITE! 6,5011 J,VOL,WT 
VTOTAL=VTOTAL+VOL 
WTOTAL=WTOTAL + WT 
99 CONTINUE 

WR I TE ( 6, 502 > VTOTAL , WTOTAL 

RETURN 

END 
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non 


$ I BFT C XSQAO DECK 


C 


SUBROUTINE XSECQD 
CONDON /0 ADA TA / 


1 

NSET 

( 16) , 

R SO 

(21) , 

ZSO 

(21) , 

PHIS 0 (21) 

2 

NBNDZNt 50 ) , 

n:mpzn(50) , 

LBD 

(6,50) , 

NT RYZ N( 6,50) 

3 

NEQBD 

( 50) t 

ABD 

(50) , 

BBD 

(50) , 

CBO (50) 

4 

XOBD 

(50) , 

YOBD 

(50) , 

ZOBD 

(50) , 

OBO (50) 

5 

HYDRATI30) , 

XSEC N 

(30) f 

XSECG 

(20,30) , 

8 I LOO (30) 

6 

BIL01 

(30) , 

BILD2 

(30) » 

BILD3 

(30) , 

GAMEN (30) 

7 

CONV 

(30) , 

COMP 

(20,40) , 

EBAR 

(30) , 

MATZ (20) 

8 

FL 

(21) , 

FM 

(21) , 

FN (21) , 


9 

ASO 

, XI1S0 

, XI2S0 

, ETA1SC 

1, ETA2S0, ALFl 

, ALF2 , A 

X 

AL F4 

, EPSLN 

, FUDGE 

, NBLO 




COMMON 

/STDATA/ 






I 

NAME 

(100) , 

ENRG 

(20) , 

XSECO 

(20,100) 

,XSECEA(20) 

2 

C (A, 6, 9) , 

EMIN33I8) , 

NAMB ( 8 ) , 

XSECNT(IOO) 


ALF3 


, DCONST 


DIMENSION BILD(30,4), XSEC ( L 9 , 3 0) , SLOPE ( 1 9 ,30 ) , S LOPEA ( 19) 
EQUIVALENCE < N SET( 7 > , NRG Y> , (NSEK4) * MAT! 

CALCULATE QAD BUILDUP COEFFICIENTS. 

NBL D= 1 , 2 , 3, 4, = H20, AL, FE , OR PB DOSE BUILDUP 

N8L D= 5* 6, 7, 8, = H20, AL , FE , OR PB ENERGY ABSORPTION BUILDUP 

DO 12 1=1 , NR GY 

E= EBARC I ) 

l F( NBL D. GE .8 I GO TO 5 
L IM = 5 

IF C NBLD.EQ.4) LIM = 6 

I F( E.LT.EMINBBINBLD)) E=EMI NBB < NBLO) 

I F( NBL D.N E .4 ) E = 1.0/F 
NB=N BL D 
GO TO 9 

5 IF( E-4.0) 6, 6» 7 

6 L IM =4 

IFtE.LT. 0.5) E =0.5 
E= l .0 / F 
NB = 8 
GO TO 9 


7 L IM =3 
N B = 9 

9 DO 10 11*1,4 

BILDC I, II )=0 • 

DO 10 J = 1 , L IM 

10 B IL D( I f I I ) = B ILO ( 1 1 1 1 ) + CdltJyNB) *E**IJ-1) 
BIL D0( I)= B ILD ( 1,1) 

BILDK I)= B ILD ( 1,2) 

B I L D2 ( I)= B ILD ( I , 3) 

12 BILD3( I)= B ILO (1,4) 

SELECT CROSS-SECTIONS FROM LARGE TABLE 
DO 15 M= 1 ,MA T 
L = MATZ( M ) 

DO 15 NE = 1 , 19 
XSEC(NE,M )= XSECO C NE , L ) 
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15 SLOP E ( NE, M ) = ( XSECO(NE+l»L)-XSECO(NE,L) ) / I ENRG ( NE+ 1 )-E NRG( NE )> 
DO IT 1=1,19 

CALCULATES SLOPES FOR ENERGY ABSORPTION CROSS-SECTION 
17 SL OPEA ( I ) = ( XSECEA ( I 1 )-XSECEA( I) ) / ( ENRG ( I + lt-ENRGUM 
CALCULATE CROSS-SECTIONS AT ENERGY EBAR(I) 

DO 25 1 = 1 , NR GY 
DO 20 NE= 1,19 

IF< EBARI I )-ENRG(NE + l) ) 22,22,20 
20 CONTINUE 
N E = 19 

22 CON VI I )=OCON ST* ( XSEC E A ( NE ) ♦ SLOPEA ( NE) * ( EBARII) -ENRG! NE ) ) > 

DO 25 M=1 , MA T 

XSECG(M,I) = X SEC ( NE , M ) ♦ SLOPE ( NE , M) * ( EBAR (I) -E NRG ( NE ) ) 

25 CONTINUE 
RETURN 



*IBFTC FLQAD DECK 


SUBROUTINE FINAL ( ID , NSL ICE ) 
COMMON /QADHD/ 


1 

Z (100) 

♦ 

vz 

(16) 


wz 

(16) 

♦ 

VR 

(16) 

2 

WR (16) 

? 

VP 

(16) 

* 

WP 

(16) 

» 

VOL 

(100) 

3 

DTIHE (100) 

♦ 

RANGE 

(100) 

t 

vv 

(100) 

V 



4 

SR G (16) 

t 

SRN 

(16) 

» 

SPG 

(16) 

t 

SPN 

(16) 

5 

SG (100) 

t 

$N 

(100) 

» 

ST 

(100) 

» 

QNDOT 

(100) 


6 OGDOT (100) , QTOOT (100) , TGDOT (100) , TNDOT (100) , 

7 TTDOT (ICO) , DQG (100) , DON (100) , OOT (100) ♦ 

8 DTG (100) , OTN (100) , DTT (100) , OG (100) , 

9 ON (100) , QT (100) , TG (100) , TN (100) , 

X TT ( 100) , T (100) 

DIMENSION 10(12) 

8 FORM AT ( 1H ,15, 1P7E12.4) 

9 FORMAT ( 1H ,15, 1P6E12.4) 

10 FORMAT ( 1H1 , 5X» 12A6) 

90 FORMAT ( 37H0 WITH TANK LIQUID LEVEL AT Z ( I ) ,/ 1H0 , 19X , 70HRAT E OF 

A HEAT INPUT TO TANK, C AL/SEC RATE OF TEMP RISE-DEGREES C/SEC/90 

BH I Z(I) GAMMA NEUT TOTAL GAMM 

CA NEUT TOTAL) 

91 FORMAT ( 32H0 WHEN I TH SLICE IS EMPTYI NG/1H0 ,12X .21HHEAT INPUT — 

F-CAL3R IES,15X» 29HTE MPE RA TURE RISE DEGREES CENT /78H I GAM 

GMA NEUT TOTAL GAMMA NEUT TOTAL) 

92 FORMAT (46H0 ACCUMULATIVE VALUES AT TIME T ( J ) ,S ECONDS/ 1H0 , 24X , 

J21HHEAT INPUT CALORIES ,1 5X,29HTEMPERATURE RISE DEGREES CENT /90H 

K J T ( J ) GAMMA NEUT TOTAL GAMMA 

L NEUT TOTAL) 

DATA RH0,CP/.07,2.25/ 

CR = RHO*CP 
JJ=NSL ICE +1 
QGDOT ( 2) = SG( 1 ) 

ONDOT ( 2) = SN( 1) 

QTDOT( 2)=ST( 1) 

QGDOT ( 1 )=QGDOT ( 2 ) 

QNDOT ( 1 ) =QNQOT ( 2 ) 

QTDOT( 1)=QTD0T(2) 

VV( 1 ) = VOL ( 1) 

DO 50 1 = 3, JJ 

QNDOT ( I)=QNDOT( I- IH-SNU-l ) 

QGDOT ( I) =QGDOT ( I- 1 1 +SG 1 1 - 1 ) 

QTDOT ( I)=QTDOT{ I-lJ+STCI-l) 

50 VV( I- l )=VV( I- 2 ) +VOL ( I - 1 ) 

DD 52 1=2, JJ 
V = 1 ./( CR*VV( I- 1) ) 

TGDOT ( I)=QGDOT ( I )*V 
TNDOT ( I )=QNDOT ( I )*V 
52 TTDOT ( I ) =Q TDOT ( I)*V 
T GDO T ( 1 ) = TGDOT (2) 

TNDOT ( 1 ) = TNDOT ( 2 ) 

TTDOT ( 1 ) = TTDOT (2) 
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WRITE! 6,10)( IDU 1,1=1,12) 

WR ITE C 6,90) 

WRITE! 6,8)(I,Z(I),QG30T(I) ,QND0T(I) , QTDOT ( I) , T GDOT ( I) , T NDOT (I), 
XTTDD T { I), 1=1, JJ) 

DO 55 1 = 1 ,NSL ICE 

AVQG= . 5* ( Q GDOT ( I ) +Q3D0T ! I ♦ 1) ) 

AVQN=.5*(QND0Tl I ) +QNOO T ! I + 1 > ) 

AVQT=.5*( QTDOT (I)+QT30T( 1*1) ) 

AVTG= . 5* t TGDOT! I ) +T3D0T! I ♦ 1 ) > 

AVTN=. 5*1 TNDOT! I H-TNDOT { 1 + 1 ) ) 

AVTT=. 5*! TTOOT ( I ) +TTOOT ( I + l ) ) 

DQG< I ) =AVQG*DT JME ( I ) 

OQN ( I)=AVQN*OT IME( I) 

DOT ( I)=AVQT*DTIME!I ) 

OTG! I ) =AV TG*DT IHE ! I ) 

OTN ( I ) =A VTN* DT IME ( I ) 

55 OTT( I ) =AVTT*DT IME ( I ) 

WR ITE ! 6, 1 0 ) l ID 1 1 ) , I = 1 , 12) 

WR ITE ( 6, 91 ) 

WRITE! 6,9) { I ,DQG< I ) ,3QN! I ) ,DQT « I) ,DTG( I) ,DTN! I ) , OTT ( I ) , 1= 1 , NSL ICE ) 
REVERSAL 

QG( l)=DQG(N SLICE) 

QN ( 1 ) = DQN ( NSL ICE ) 

QT! 1)=DQT!NSLICE) 

TGI 1)=DTG(NSLICE) 

TN! 1 ) =DTN ( NSL ICE ) 

TT ( 1 ) =DTT { NSL ICE ) 

T! 1 ) = DT IM E !N SL ICE ) 

00 60 1=2, NSL ICE 
N=J J— I 

Tl I) = T( I-D+DTIMEIN) 

QG( I ) =QG! I-IUDQG(N) 

QN! I ) =QN( I- 1 ) +OQN ( N ) 

QT ( I ) =QT( I-1M-0QTIN) 

TGI I ) =TG( 1-1 )*DTG(N) 

TN! I ) = TN ( 1-1 M-DTN(N) 

60 TT! I ) = TT( I- 1 ) +DTT ( N ) 

WR ITE ! 6, 1 0 ) ( 10(1 ), 1 = 1, 12) 

WR ITE! 6,92) 

WRITE! 6, 8 * ( I , T ( I ),QG( I) ,QN(I ) ,QT(I ) , TG ( I ) ,TN(I ) ,TT !I), 1 = 1, NSL ICE) 

RETURN 

END 


tIBFTC SCQAD DECK 


SUBROUTINE SOURCE 
COMMON /QADATA/ 


1 

NSET 

(16) , 

RSO 

(21) , 

zso 

(21) 

9 

PH IS 0 

(21) 

2 

NBNDZNt 50 ) , 

NC MPZN ( 50) , 

LBD 

(6,50) 

9 

NTRYZNC6, 50) 

3 

NEQBO 

(50) , 

ABO 

(50) , 

BBD 

(50) 

9 

CBD 

(50) 

4 

XOBO 

(50) , 

YOBO 

(50) , 

ZOBD 

(50) 

9 

DBD 

(50) 

5 

HYDR A T( 30 ) , 

xse:n 

(30) , 

XSECG 

(20,30) 

9 

BILDO 

(30) 

6 

BIL01 

(30) , 

BILD2 

(30) , 

BILD3 

(30) 

9 

GAMEN 

(30) 

7 

CONV 

(30) , 

COMP 

(20,40) , 

EBAR 

(30) 

9 

MATZ 

(20) 

8 

FL 

(21 ) , 

FM 

(21) , 

FN (21) 

t 



9 

A SO 

, XI1S0 

, XI2S0 

, ETA1S0 

1, ETA2S0, ALF1 

, ALF2 f A 

X 

AL F4 

, EPSLN 

, FUDGE 

, NBLD 







ALF3 


EQUIVALENCE! L SO *NSET( 1)) , ( MSO ,NSET(2 ) ) , ( NSO, NSET (3 ) I , (ISRC, 
I NSET(ll)) 


40 FORMAT (32H1 L 

41 FORMAT ( 32HO M 

42 FORMAT (32H0 N 

FL T=0 • 0 

FMT =0 • 0 


RCL) 
Z(M> 
PHI (N) 


FILI//II7 »5X1P2E12.4) ) 
F(M)//(I7 »5X1P2E12.4) ) 
F(N)//( 17 ,5X1P2E12.4) ) 


DO 15 L=1 » L SO 
GO TO (10,13), ISRC 
10 IF(XIISO) 12, 11,12 

12 FL ( L » * (CO SC XI1S0*(RS3(L*1)-XI2S0))-C0S(XI1S0*(RS0(L)-X I2S0 ) ) )/ 
11S0**2+(RS0( L* 1 )* SI N( XI ISO* (RS0(L+1 ) -XI2 SO) ) 

20(L)*SIN( XI1S0*(RS0(L)-XI2S0) ) ) /XI ISO 
GO TO 14 


, 


XI 

-RS 


11 FL(L ) = (RSO ( L +1 )**2-RS0(L ) **2) 

GO TO 14 

13 FL(L)*(FL(L+1)-FL(L) )* ( RSO ( L* 1 ) * (RSO ( L) + RS C( L+l ) )+RSO( L )**2 )/ 1 . 5 
1+(RS0( L>+RS0(L+1>)*(FL(L)*RS0(L«-1J-FL(L+1) *RSO(L) ) 

14 FLT=FLT*FL(L ) 

15 RSO(L )-IRS0(L)+RS0(L+l))/2»0 
DO 25 M= 1 » MSO 

GO TO (20,23), ISRC 

20 IF( ETA ISO >21,22,21 

21 FM(M)=(SIN(FTA1 SO* ( ZSO ( M+l ) -E TA2S0) ) 

IS IN < ETA1S0*(ZS0(M)-ETA2S0) ) 1/ETA1S0 

GO TO 24 

22 FM( M ) = ZSO (M-fl)-ZSO(M) 

GO TO 24 

23 FM(M)=(Z$0(M+1)-ZS0(M) ) * ( FM( M) *FM( M+l ) ) 

24 FMT=FM T+FM ( M ) 

25 ZS0(M)=(ZS0(M)+ZS0(M+1) )/2.0 

ASO=A S0/( FL T*F MT*IPHI SO (NS 0+1)- PH I S0(1 ) ) ) 

DO 30 N=1 ,NSO 

FN(N)=(PHIS0(N+1)-PHI SOI N) )*ASO 
30 PHIS0(N)=(PHIS0(N)+PHIS0(N+1)) /2.0 
WRITE (6,40) (L ,RSO( L) , F L (L ) ,L*1 ,LSO) 

WRITE (6,41) (M,ZS0(M) ,FM(M) ,M = 1,MS0) 

WRITE (6,42)(N,PHISO(N),FN(N),N=1,NSO) 

RETURN 

END 
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$ I BFTC LNQAD DECK 


SUBROUTINE LENGTH( SMA T ,S SO , X SO, YSO, ZSO , XRC ,YRC ,ZRC, IERR) 

COMMON /DADATA/ 

1 NSET (16) , R SO (21) , ZSO (21) , PHISO (21) , 

2 N BMDZNI 5 0 ) , NC MPZN ( 50) » LBD (6,50) , NTRYZN(6,50) , 

3 NEQBD (50) , ABD (50) , BBD (50) , CBO (50) , 

4 XOBD (50) , YOB 3 (50) , ZOBD (50) , DBO (50) , 

5 HYDRAT(BO) , XSECN (30) , XSECG (20,30), 8ILD0 (30) , 

6 BILD1 (30) , BILD2 (30) , BIL03 (30) , GAMEN (30) , 

7 CDNV (30) , COMP (20,40), EBAR (30) , MATZ (20) , 

8 FL (?1) , FM (21) , FN (21) , 

9 ASO , XI1S0 , XI2S0 , ETAISO, ETA2S0, ALF1 , ALF2 , ALF3 , 

X ALF4 , EPSLN , FUDGE , NBLO 

DIMENSION R ( 6 ) , SMAT (20) 

EQUIVALEMCE(NREG,NSET( 6) ) , ( NZSO ,NSE T( 10) ) , ( MAT , NSET (4 ) ) 

5 FORMAT (67H0 ZONE BOUNOARY DISTANCE X Y 

1 Z ) 

6 FORM AT t 2( 17, 5X ) , 1P4E 12. 4) 

9 CALL SLITET( 4,K000FX) 

GO T0(10,10), KOOOF X 

10 CALL SLITET( 1, KOOOF X) 

GO T0(101,102), KOOOF X 

101 CALL SLITE ( 1) 

WRITE (6,5) 

102 DO 11 M=1,MAT 

1 1 SMAT(M )=0.0 
X = X S3 

Y =Y S3 
Z =Z S3 ( 1) 

AL P HA * XRC - X 
BET A = YRC- Y 
GAMMA = ZRC-Z 

SRC=SORT( ALPHA**2+BETA**2+GAMMA**2 ) 

AL P HA = AL P HA / SR C 
BET A=R ETA / SRC 
GAMMA* GAMMA /SRC 
SSO=SRC 
SM IN*FUOGE 
NP=NZ SO 

100 X*X +AL PHA* SM IN 
Y = Y ♦BETA* SM IN 
Z *Z +GAMMA* SM IN 
SRC*SRC-SM IN 
20 IL =NP 
IU=NR EG 

25 DO 42 I* IL , I U 
JU*NBN DZN ( I) 

DO 41 J=l, JU 
K*LBD( J, I ) 

NEQ=N E0BD(K > 

GO T3( 31,32,33,34,35,36) ,NEQ 
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31 R ( J )=X*(ABD(K)*X +XOBO ( K) )+Y*(BBD(K> *Y+YOBD(K) ) 
l+Z*(CBDtK )*Z +Z OBD ( K ) ) -DBD IK) 

GO TO 40 

3 2 R{ J )=ABD( K )* ( X- XOBD ( K) )**2*BBD(K)*(Y-Y0BD(K) ) **2 
1+CB01K )*(Z-Z0BD(K) )* * 2-0 BD ( K) 

GO TO 40 

33 R l J ) = ( X-XOBO ( K ) >** 2* ( Y- YOBO ( K) ) **2-DBD (K) 

GO TO 40 

34 R( J )=X-OBD(K ) 

GO TO 40 

35 R( J ) = Y-DBD(K ) 

GO TO 40 

36 R( J ) =Z— DBO ( K ) 

40 IF( SIGN! 1.0,LB0{ J,I ) )*R< J) ) 41,41,42 

41 CONTINUE 
GO TO 50 

42 CONTINUE 

I FI IU-NP >45,43,43 

43 IL = 1 
IU=NP- 1 
GO TO 25 

45 I ERR = 1 

CALL SLITETl 1,K000FX) 

GO TO ( 300 , 46 ) , KOOOF X 

46 CALL SLITE ( 1) 

GO TO 9 

50 SM IN=SRC 
KK =0 

I F ( NBNDZN ( I) ) 99,99, 51 

51 DO 98 J=1,JU 
K=LBDl J, I ) 

N EQ=N EOBD ( K ) 

GO T0( 61,62,63,64,65,66) ,NEQ 

61 E=( ABDIK)*X+XOBOIK) /2.0)*ALPHA+ ( BBD ( K) *Y*YOBD ( K) /2 .0 > *BET A 

1 + (CBD(IO*Z+Z0BD<K)/2.0)*GAMMA 

GO TO 620 

62 E=ABD( K)* ALPHA* I X-X03D IK) ) +BBD « K) *B ETA* ( Y-YOBD I K ) ) 
l+CBD( K )*GAMMA* IZ-ZOBD t K 1 ) 

620 H=ABD( K)*ALPHA**2«-BBD( K ) *BE TA**2+CBD IK) *GAMMA**2 
GO TO 70 

63 E = ALPH A*( X-XOBDI K )) +BE TA* ( Y-YOBD I K) ) 

H=1.0-GA'IMA**2 

GO TO 70 

64 E=ALP HA/2 .0 
GO TO 660 

65 E=B ETA /2.0 
GO TO 660 

66 E=GAMMA/2 .0 
660 H=0 .0 

70 I F ( H) 7 1, 75, 75 

71 H=-H 
E=-E 

R( J )=-R( J ) 

K=-K 

75 I F ( H-EPSLNI76, 76,85 


76 I FI R ( J ) )80, 77, 80 

77 I ERR = 2 

CALL SLITET( 1,K000FX) 

GO TO C 300 » 7 8 ) ,KOOOFX 

78 CALL SLITE < U 
GO TO 9 

80 IFC ABSCE)-EPSLN)98,98» 81 

81 SBD=*-RCJ> /2.0/E 
I F{ SBD)98,98,95 

85 IFCRC J 1)86,77,90 

86 IFC H*$MIN**2+2.0*£*SMIN+R(J) ) 98,98,87 

87 QUAD=E**2-H*RC J) 

GO TO 92 

90 IFC 6)91,98,98 

91 QUAD=E**2-H*RC J) 

IFC QUA D-EPSLN) 98,98,92 

92 S8D-C-E+S IGN l SQRTC QUAD ) ,K) ) /H 

95 IFC SBO-S^I IN )96,98,98 

96 CALL SLITE C 4) 

SN IN = S8D 

KK=K 

NP -NTR YZN ( J , I ) 

98 CONTINUE 

99 SM IN= SM IN +FUDGE 

CALL SLITETC 1,K000FX) 

GO T0(991,992),K000FX 

991 CALL SLITE ( 1) 

WRITE (6,6)I,KK,$MIN,X,Y,Z 

992 N=NCMP ZN ( I ) 

DO 994 M= 1 ,MAT 

994 SMATCM ) = SMATCM)+SM IN* COMP ( M, N) 

CALL SLITETC 4, KOOOFX) 

GO TOC 100, 105), KOOOFX 
105 I ERR-0 
300 CALL SLITE (0) 

RETURN 

END 
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$ I BFTC KNQAD DECK 


SUBROUTINE KERNEL(SMAT,SSO,ANEUT,AGAM,UAGAM, IERRI 
COMMON /3ADATA/ 


1 

NSET 

( 16) 


R SO 

(21) 


zso 

(21) 

t 

PHIS 0 

(21 ) 

9 

2 

N BNDZ N< 50 ) 


n:mpzn 

(50) 


LBO 

(6,50 ) 

9 

NTRYZN( 6 , 50 ) 

9 

3 

NEQ8D 

< 50) 


AB3 

(50) 


BBO 

(50) 

9 

C80 

(50 ) 

9 

4 

X OBO 

( 50) 


YOBO 

(50) 


ZOBD 

(50) 

t 

OBO 

(50) 

9 

5 

HYDRAH30) 


XSECN 

(30) 


XSECG 

(20,30) 

9 

B1LD0 

(30) 

f 

6 

BIL01 

(30) 


BILD2 

(30) 


BIL03 

(30) 

9 

GAMEN 

(30 ) 

» 

7 

CON V 

(30) 


COMP 

(20,40) 


EBAR 

(30) 

» 

MATZ 

(20) 

9 

8 

FL 

(21) 


FM 

(21) 


FN (21) 

9 




9 

ASO 

, XI1S0 


9 XI2S0 

, ETA1 

SO 

, ETA2S0, ALF1 

, ALF2 , ALF3 

X 

ALF4 

t EPSLN 


♦ FUDGE 

, NBLD 









DIMENS ION SMATI 20) ,AGAM(20) ,UAGAM(2 0) 

EQU I V ALENCEI MA T,NSET( 4) ) , (NRES ,NSET (6) ) , ( NRGY * NSET { 7 ) ) , 
1( INEUT »NSET I 12) J 
AN EUT = 0.0 
I F ( INEUT) 21, 21,15 
15 SHY0=0 .0 
AN =0.0 
DO 18 1=1, MAT 
AN = AN + SMAT( I >*XSECN(I ) 

18 SHY D=S HYO+HYDR A T I I )*SMAT ( I ) 

AN EUT = ALF 1#EXP I-AN) 

I F ( SHYD) 19,19, 20 

19 I ERR = 1 
GO TO 22 

20 AN EUT =ANEUT* I SHYD ) **ALF2*£ XPI-ALF3 *( SHYD >**ALF4) 

21 I ERR = 0 

22 DO 25 J = 1 ,NR GY 
AG=0 . 0 

DO 24 1=1, MAT 

24 AG=AG+SMA T( I )* XSECG ( I , J ) 

UAGAM ( J )=EXP (- AG ) 

25 AGAMI J >=( BILDOC J) +9 1 LDlt J)*AG+BIL02 ( J) *(AG**2> 
H-BIL03(j)*(AG**2)*AG)*UASAM( j) 

RETURN 

END 


S I 8FTC X1SEC0 DECK 


BLOCK DATA 
COMMON / S T DA T A / 

1 NAME (100) , ENRS (20) , XSECO ( 20 , 1 00 ) , XS ECEA( 20 ) , 

2 C (4,6,9) , EMIN3B(8) , NA MB ( 8 ) , XSFCNT(IOO) ,DCONST 

COMMON /GAUSSQ/ V(4,16), W(4,16) 

COMMON / DOSE/ XSECE8 (20,2), A LF ( 2 ) , DC0NI2) 

DATA ALF/ 9.1429E-14, 2. 1 87E-3 /,DCON/ 2 . 6755 68E-15 , 5.76756E-5/ 

DATA (XSECFB( 1,1), 1 = 1, 20)/ 0. 0270,0.0305,0.0362,0.0406,0.0485, HYO.ENRG 
AO. 053 0, 0.0573, 0.0 58 7, 0. 0589, 0. 0588 , 0. 057 3 , 0. 05 55 ,0.0507,0.0464, ABS.XSEC 
BO. 038 9, 0.0351, 0.0316,0.02 88,0. 0249,0.0222/ NBS-H85H 

DATA (XSECFBJ 1,2) ,1 = 1,20) / 0.0409,0.0312,0.0255,0.025 2,0.0 276, MUS.ENRG 

CO. 029 7, 0.031 7, 0.0325, 0. 032 7, 0. 0326 ,0. 031 8 ,0. 0308 ,0 .028 1 ,0 .0257 , ABS.XSEC 

DO. 022 5,0. 0203, 0.0188, 0. 01 78 ,0. 0163 ,0.0154/ N3S-H8 5M 

DATA XSECNT / C.0 , 0. 14989, 0. 08763 ,0. 071 53 ,0. 05 40 1 ,0 .04063 ,0 .04283, 1-7 

1 0.0 3 728,0.04091,0.0 343 7,0. 03213 , 0. 031 2 l , 0. 02 925 ,0 .02898 , 0 .02 755, 8-15 

2 0.02706,0. 02570,0.02416,0.02443,0.02412 ,0.02270,0.02200,0.02131, 16-23 

3 0.0 2109,0.0 2 050,0.0 2136,0.01 977,0. 01940,0.01934,0.01874,0.018 13, 24-31 

4 0.01775,0.01747,0.01700,0. 016 90,0.01649,0.01632,0.01611,0.01599, 32-39 

5 0.01559,0.01563,0.0153 7,0. 01513 ,0. 01496 ,0.01483,0.0145 5,0.01447, 40-47 

6 0.01417,0.01401,0.01378,0. 01 360,0 . 01 32 7 ,0. 01 33 1 ,0 .0 1 308 , 0 .0 1 299, 48-55 

7 0.01278,0.01270,0.01264,0.01261 ,0.01246,0.01242 ,0 .01219,0.01212, 56-63 

8 0.01 191 ,0.01185,0.01171 ,0.01162,0.01154,0.01 148 ,0 .01134,0 .01 128, 64-71 

9 0.01 116,0.01108,0.01101 ,0.01092 , 0. 01 080 ,0. D 1 074 ,0 .0 1066 , 0.0 106 l, 72-79 

X 0.01C51, 0.01041, 0.01 026, 0.01 006, 0.01 026, 0.0 1024, 0.00997, 0.00994, 80-87 

Y 0.00588, 0.00986,0. 00974,0. 00 97 7,0. 00911 ,0.03964,0.00954,0.00952, 88-95 

7. 0.00 548 ,0. CO 940 ,0. 00940 , 0. 0093 0,0. 00934/ 96-100 

OATA( V ( 1, n, 1=1,2) /. 21132487, . 78867513/, (Wt 1,1) ,1=1, 2)/ 2*. 5/ 

DA T A ( V ( 2 , I ), I =1,4)/. 06943 184,. 3300094 8,. 669990 52 ,.93056816/, 

E( W ( 2, I ), 1=1, 4)/. 173 92742, 2*. 32 60725 8 ,.17392742/ 

OATA( V (3, I ), I = 1,8)/. 01 98 850 7,. 101 666 76,. 2 3 723 3 79,. 40828 268, 

P. 59171732, .762 766 21,. 8 983 3324,. 98014493/, (W (3 , I ), 1= 1 , 8 )/ .05061 4268 

Q, .1111 9052 , . 15685332, 2*. 181341 85, .1 5685332 ,.11 119052,. 0 50 614268/ 

DATA( V (4, I), 1 = 1, 16)/ 0. 52 9953E-2 ,0. 2 7712 5 E-l ,0 . 67 18 44E- 1 ,0 . 1 22298, 

XO. 1910 62, .270992,0.359198,0.4 52 494,0.547506,0.64080 2,0.7290 08, 

YO. 808938, .877702,0.932816,0. 972288, C.9947/ 

DATA! W (4, I), 1=1, 16) /O. 01 35 762 ,0. 0311268,0.0475793,0.0623145, 

XO. 0747 980, 0.084578 3, 0.091 301 7,2*0.0947253 ,0.0913017,0.0845783, 

YO. 0747980, 0.06 23145, 0.0475793 ,0.0311268,0.0135762/ 

DATA (NAMEU), 1=1,100) / 2H H ,2HHE,2HLI ,2HBE,2H B,2H C,2H N,2H 0, 1-8 

S2H F,2HNE,2HNA,2HMG,2HAL,2HSI ,2H P,2H S,2HCL ,2 HAR, 2H K,2HCA,2HSC, 9-21 
T2HT ( , 2 H V , 2HCR , 2HMN, 2HFE,2HC0 * 2HNI ,2HCU,2 HZN, 2HGA, 2HGE, 2H AS , 2HS E, 22-34 
U2HBR , 2HKR,2HRB,2HSR,2H Y,2HZR»2HNB,2HM0,2HTC »2HRU,2HRH, 2HPD, 2HAG, 35-47 
V2HCD, 2 HIN , 2HSN , 2HSB » 2HTE , 2H I ,2HXE , 2HCS , 2 HBA , 2HL A, 2HCE, 2HPR, 2HND, 48-60 
W2HPM, 2HSM , 2HEU » 2HG0 » 2H TB » 2HDY ,2HH0, 2HER ,2HTM,2HYB,2 HLU , 2HHF, 2HT A, 61-73 
X2H W, 2WE,2HOS,2HIR,2HPT, 2HAU ,2HHG , 2HTL ,2H PB ,2 HB I , 2HP0 , 2HAT , 2HRN, 74-86 
Y2HFR, 2HRA, 2HAC »2HTH,2HPA ,2H U , 2HNP , 2H PU, 2H AM, 2 HC M, 2 HB K , 2HCF, 2HES , 87-99 

Z2HFM / 100 

DATA EMINBB/ 0.255, 3*0.5, 0.255, 3*0.5/ , 

UC(I, 1,1) , 1 = 1, 4) 7 + 1.010945*0, + 1.16772 E-l , -7. 65869E-3 , + 1 .67068E-4/, 

2(C( I, 2,1) , 1=1, 4)/-6.00394E-2, + 2.321 2 5E+0 ,-l .790 23F-2 , +5 .69 2 95E-4/, 
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31 C( I, 3,1) ,I=l,4)/+7.20778E-2,-2.128 01E+0,+2.41735E-l,-7.96 332E-3/, 
4<CU,4|1I, 1*1, 4)/-3.01498E-2,+7.67783E-l,-4.34443E-2,+7.23758E-3/, 

5(C( I, 5, 1 ) , 1= 1, 4) /+3. 9 4733E-3 ,-9. 081 39E-2 ,-l. 34203E-3 »-9 .87237E-4/* 

If C< I, 1,2) ,1*1, 4)/+1.00768E+0,-1.03807E-2,+1.33705E-2,-3.29348E-4/, 

2(C( 1,2,2), I=l,4)/-4.98085E-2,+3.32216E+0,-1.58167E-l,+4.60315E-3/, 

3C CM, 3,2) ,1*1, 4) /+7.23425E-2 ,-5. 52427E+0»+6.39496E-l ,— 2.04255E-2/, 

4( Cl I, 4,2) ,1 = 1,4)/- 3. 93 841 E -2 * + 4. 167 COE+O »-5 .59836E-1 ,+2 .00554E-2/, 

5( C( I, 5, 2), 1*1,4)/+ 7. 3577 8E -3 »-l . 04638E+0 ,♦ 1.4 1 308E— 1 ,-5 .29934E-3/ , 

1(CU, 1,3) ,I=l,4)/+1.01460E+0,-4.l2104E-2,+1.8B074E-3,+l.20198E-3/, 

21 Cl 1,2,3), 1=1, 4)/- 1.8865 7E-1 , + 2.72 752E+0,+1.03217E-l,-9.83313E-3/, 

3f C I I, 3,3) , 1=1, 4) /+6. 38649E-1 »-3.76728E+0»-1.319B8E-l»+2.06002E-2/* 

41 Cl 1,4,3) , I=l»4)/-6.55159E-l*+2.42384E+0,+1.68976E-l»— 1.75251E— 2/» 

51 C ( I, 5,3) ,1*1,4 )/+1.90742E-l ,-5. 546 57E-1 ,-5.80710E-2 ,+4 .75673E-3/ 

DATA NAMB/6HH20 DS, 6H AL DS, 6H FE DS, 6H PB DS , 

X 6HH20 EA, 6H AL EA, 6H FE EA, 6H PB EA/ , 

1C Cl I, 1,4) , 1*1,4 )/+9. 59342E-1 ,+6. 782 54E-2 ,-2 .266 26E-2 ,+6 .39872E-4/, 

21 Cl I, 2,4) » I* 1»4)/+1.13722E-1 ,+4. 50412 E-l ,+7.55191E-3,+1.52094E-4/, 
3lClI»3»4)»I=l»4)/-7.39816E-2»-2.15037E-l,+5. 10254E-3,-6 .04837E-4/, 

41 Cl 1,4,4) ,1*1,4) /+1. 87 76 7E -2, + 4. 051 89E-2 ,— 1 .89332E-3, +2.42263E-4/, 

51 Cl I, 5,4) , I* 1, 4) /- 2. 04254E-3 *-3.40802E— 3*+1.93415E-4,-2 .93865E-5/, 

61 Cl I, 6, 4), 1*1, 4) /+7. 93621E-5 , + l . 06510E-4,-6. 24306E-6 , + 1.1 3914E-6/, 

11CI I, 1,5) ,I*l,4)/+1.0093 5E + 0,+ 6. 75909E-2 ,-2.88 512E-3 ,+3.541 58 E-5/, 

21 Cl I, 2, 5), 1*1, 4) /-5.051 87E-2 , + 2. 739 52E+0,-7. 1 1806E-2 ,♦ 1.98272E-3/, 

31 Cl 1,3,5) ,1*1,4) /+ 5 . 9 3 8 63E-2, -2.823 06E+0, + 3.98796E-1,-1 .19564E-2/, 

41 ClI, 4, 5), 1 = 1, 4)/- 2. 4371 7E-2, + l.ll668E+0, — 1.50948E— l,+9.52922E-3/» 

51 Cl I,5,5),I = l,4)/+3.14779E-3,-1.41350E-l, + 1.65 26 6E-2.-1 .36541E-3/, 

11 Cl I, 1,6) , 1 = 1, 4) /+9. 97130E-1 ,-3 . 852 19E-2 , +1 .278 16E-2 ,-3 .39093E-4/ , 

21 Cl 1,2,6) ,I = l,4)/ + 6. 54640E — 2 , + 3. 27946E+0 ,-l . 64393E-1 ,+5 .04304E-3/, 

31 Cl I, 3,6) ,I = 1,4)/-2.68716E-1 ,-4. 74349E+0 , + 7.45795E-l ,-2 .38346E-2/, 

41 Cl 1,4,6) ,1=1,4) /+2. 8490 8E-1.+3. 394 95E+0,-6.26720E-l ,+2 .53658E-2/, 

51 Cl I, 5, 6), 1 = 1, 4)/- 8. 331 77E-2 ,-8. 373 92E-1 , + 1.66 032 E-l ,-7 .07268E-3/ 

DATA 

11 Cl I, 1,7) ,I=l,4)/+1.00009E+0,-9.76981E— 2.+3.07 480E-3 ,+5.43036E-4/ , 

21 Cl 1,2,7) , I=l,4)/-2.69391E-2,+2.81567E+0, + 5. 4 7069E— 2 ,— 3 .058 50E-3/, 

31 Cl 1,3,7) , I=l,4)/+1. 06797E-1 ,-3. 145 90E +0,-1. 33763E— 4 , + 2 .4368 IE-3/, 

41 Cl I, 4, 7), 1*1, 4)/- 1.11 72 5E-1 , + l . 91 5 66E+0, + 1.0 3003E-1 ,-l . 570 16E-3/, 

51 Cl I, 5,7) , I = l,4)/+3.27028E-2 ,-4. 1 7820E— l ,-4.79248E-2,+5.351 1 IE— 4/, 

II Cl I, 1,8) ,I = l,4)/+9.99305£-l,-3.49764E-l,+2.99886E— 2, + l.28555E— 3/, 

21 Cl I, 2,8), I* l,4)/+1.13017E-2, + 2.79 914E+0,-7.10002E— 2 ,-4. 75607E— 3/, 

31 Cl I, 3,8) , 1= l,4)/-3.44367E-3,-2.203 99E+0,-l. 11 108E-2,+7.04843E-3/, 

41 Cl I, 4,8) ,I = l,4)/+2. 04092E-3 , + 5. 023 01 E— 1 , + 1.45 444E-2 ,— 2 .3 1052E— 3/, 

11 Cl I, l,9),I=l,4)/+1.02 626E + 0, + 4.l5814E— 1. + 5.64 179E— 2 ,-4. 83073E— 3/ , 

21 Cl I, 2,9) ,1=1, 4)/- 7. 44724E-3 ,-6.147 01 E-2 , — 1 .42 153E-2 , +1 .76 866E— 3/, 

31 Cl I, 3,9) , 1=1,4)/+ 3. 8 1090E-4 . + 3.09 491E— 3, + 7. 33689E-4, — 1 .07370E-4/ 

DAT Al ENRGI I), 1*1, 20 )/. 05, . 06, . 08, .1,. 15, .2, .3,. 4,. 5, .6, .8,1. ,1.5, 

1 2. ,3. ,4. ,5. ,6. ,8., 10. / 

DATA! XSEC01 I , 1 ) , I = 1, 20 ) /. 3 35, .326 ,. 309, . 295 , . 265 , .243 , . 2 12, . 189, HYDROGEN 
X .173, .160, .140,. 126,. 103,. 0876, .0691 ,.0579 ,.0502 ,. 0446 , .0371 ,. 032 1 Z=1 


Y / 

DATA! XSECOl I ,2) ,1 = 1,20)/. 1688,. 1640,. 1560, . 1480 ,. 1330 ,. 1220 , .1060, HE 0201 

X. 0953, .0870, .0805 , . 0707, . 063 5 ,. 0516 ,.0442, 0202 

Y. 0349, .0295, .0257, . 0231 ,. 0194 ,. 0169/ 0203 

DAT Al XSECOl I ,3) ,1 = 1, 20)/. 1469,. 1423 ,. 1 347 , . 12 5 8 , . 1 1 50 , • 1060, .092 1, L 1 0 301 

X. 0 8 2 5, .0 7 54, . 0 6 96, . 06 1 2 , . 05 5 0 ,. 044 7, .0382 , 0 302 

Y . 0304, .0257, .0225, .0203,. 0172 ,. 0152 / 0303 


DAT Al XSECOl 1 ,4 ) , I = 1 , 20 ) /. 1 51 , . 1 46 , .138 ,. 1 3 2 , . 1 19 ,.109,. 0945, .0847,8ERYL IUM 
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X. 0 773, .0715, .0628,. 0565,. 045 9,. 0394,. 0313 ,.0 266 , .0234 , .02 1 1 , .018, Z = 4 

Y. 0161 / eE 

DATA( XSECO ( I , 5 > , I = 1 , 20 ) /. 1 608,. 1 536 ,. 1463 , . 1380 , .1237 , . 11 39, .0985, B 0501 

X. 0881, .0805, .0745, .0654,. 0588, .0478, .0410, 0502 

Y. 0321 , .02 79, .0247 ,. 0224, . 01 93 ,. 0175/ 0503 

DATA! XSECO ( 1 ,6), 1 = 1, 201/. 178, .169, .157,. 149 , . 134 , . 1 22 , . 10 6, . 0953, CARBON 

X. 087, .0805, .0707, .0636,. 0518,. 0444 , . 0356 ,.0304, .027, .0245, .0213, Z = 6 

Y. 0194 / c 

DAT A( XSECO! I, 7), 1=1, 20)/. 185, .173 ,.159, .15, .134, .123, .106, .095 5, .ONITROGEN 
X869, .0805, .0 707,. 0636,. 051 7,. 0445,. 035 7 ,.0305 , .0273 , .0249, .0218, 1=7 

Y .02/ N 

DATA! XSECO! I, 8), 1 = 1, 20)/. 197, .18, .162, .151, .134, .123,. 107, .0953, OXYGEN 

X. 087, . 080 6,. 0708,. 063 6,. 0518,. 0445 ,. 0359,. 0309 ,.0276, .0254, .0224, Z = 8 

Y. 0206 / 0 

DATA! XSECO (I ,9) ,1=1, 20)/. 1987,. 1786 ,. 1 565 ,. 1 448 ,. 1275 ,. 1780, . 1016, F 0901 

X. 0905, .0826, .0762, .0670, . 0602 ,. 0491 ,.0422, 0902 

Y . 0 342, .0296, .0266, .0245, .0218, .0203/ 0903 
DATA! XSECO! I , 10) , 1 = 1, 20) /. 231 ,. 1 992 ,. 1688 , . 1 5 36 , .1 337 , . 1226, .1054, NE 1001 

X . 0947 , .0863, .0797, .0701, . 0 630 ,. 051 4, .0442 , 1002 

Y . 0359, .0312, .028 1, . 0260, . 023 3 ,. 0218/ 1003 
DATA! XSECO! I , 1 1 ) , I =1 , 20 ) /. 24 8, . 2 06 , .168 , . 1 5 1 , . 1 3 , .1 18 , . 102 , .09 12, SODIUM 

X. 0833, .077,. 06 76, .060 8,. 0496,. 0427 ,. 0348 ,. 0303 ,.0274, .0254,. 02 29, Z = ll 

Y. 0215 / na 

DATA! XSECO! I , 12 ) , I = 1 , 20) /. 293 , . 2 32 , . 1 81 , . 16 , . 1 35 , . 1 22 , . 106, .0944, MAGNSIUM 

X. 086, .079 5, .06 99,. 0627,. 0512 ,.0442 , . 036 , . 031 5 , .0286 , . 0266, .0242, Z=12 

Y. 0228 / MG 

DATA! XSECO! I , 13 ) ,1 = 1, 20) /. 326, .24 8, .186, .161, .134, .12, .103, .0922, ALUM INUM 

X.084, . 077 7, . 0683, .0614, . 05 ,. 0432,. 0353 ,. 031 , . 028 2 , .0264, .0 241 , Z = 13 

Y .0229/ AL 

DATA! XSECO! I, 14), 1 = 1, 20) /. 3 89 ,. 2 88 , .205 , . 1 72 , . 139 , . 125 , . 107, .0954, S IL ICON 

X. 0869, .0802, .0706, . 0635 ,. 051 7 ,. 0447 ,. 0367 , .3 323 , .0296, .0277, .0254, Z = 14 

Y. 0243 / SI 

DATA! XSECO! I, 15 ), I =1,20)/. 432,. 311,. 211,. 174,. 137,. 122,. 104,. 0928, PHOS- 

X. 0846, .078, .0685,. 061 7,. 0502,. 0436 ,.0358 ,.03 16, .029, .0273, .0252, PHORUS 

Y. 0242 / z=l5 

DATA! XSECO II, 16), 1 = 1, 20)/. 51 8, .36 3, .2 34, .189, .144, .127, .108, .0958, SUL PHER 

X. 0874, .0806, .0707, . 0635, . 051 9 ,. 044 8 ,. 0371 , .0 328 , .3 302 , .0284, .0266, Z=16 

Y. 0255 / s 

DATA! XSECO! I, 17), 1 = 1, 20)/. 5 84, .3957, .2451,. 190 2,. 1403,. 1223, .1033, CL 1701 

X. 0922, .0839, .0775, .0680, .0609,. 0498, .0432, 1702 

Y . 0358, .0318, .0295, .0 280, .0262 ,. 0252 / 1703 

DATA! XSECO! I ,18), 1 = 1, 20) /. 629,. 420, .24 9, . 188 ,. 135 ,. 117 , .0977, , 0867 ARGON 

X, . 079, .07 3, .0638, .0573,. 0468,. 04 07 ,.0338 ,.0301 ,.0279, .0266, .0248, Z = 18 

Y. 0241 / A 

DATA! XSECO! I, 19), 1 = 1, 20)/. 782 ,. 514 , .2 93 , . 2 15 , . 149 ,. 12 7, . 106 , .0936, POTAS IUM 

X. 08 52, .0786, .0689,. 061 8,. 0 505,. 0438 ,. 0365 ,. 3 327 , .3305 , .0289, .0 274, Z = 19 

Y . 0267 / K 

DATA! XSECO! I , 20 ) , I = 1 , 20 ) /. 92 9, . 597 , . 33 ,. 238 , . 158 , . 132 , . 109, . 0965, CALCIUM 
X.08 76, .0809, .0708,. 0634, .051 8, .0451 ,.0376,. 033 8 ,.0316, .0302,. 0285, Z=20 

Y • 0 2 8 / 

DATA! XSECO! I ,21), 1=1, 20 )/. 999, . 631 7 ,. 3387 , .2377 ,. 1507 ,. 1247, . 1025, SC 2101 

X. 0902, .08 19,. 0 755,. 0662,. 0593,. 0484,. 042 1 , 2102 

Y. 0 353, .0318, .0299, .028 7, . 0273 ,. 0267/ 2103 

OAT A! XSECO! I , 2 2 ) , I = 1 , 20) /I . 1 25 , . 7062 ,. 36 77 , . 25 17 , . 1 53 1, . 1 253 , . 101, T I 2201 

X. 0892,. 0809, .0745,. 0653, .0586, .04 78, .04 16, 2202 

Y . 0350 , .03 16, .0298, .028 7,. 02 74,. 0269/ 2203 
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DATA( XSECO ( l,23),I=l,20)/1.250,. 785 0,. 3999 ,.2661 , .1565 ,.1249, .100, V 2301 

X. 0880, .0797, .0734, . 0642 ,. 0576 ,. 0470 ,. 0410 , 2302 

Y. 0345, .0313, .0296, .0285, . 02 74 ,. 0270/ 2303 

DATA! XSECO ( I ,241, 1 = 1, 20 1/1. 443,. 900 5,. 45 16 , .29 30 ,. 1665, .1302, . 103, CR 2401 

X. 0904, .0817, .0 751,. 0657,. 0589,. 0482 ,.0419, 2402 

Y . 0 3 54, .03 2 2, .0 3 05 , .0 2 9 5, . 0 2 8 5 ,. 0282 / 240 3 
DATA( XSECO ( I ,25 ),I =1 ,201 /l. 599, . 9942 ,. 4925 ,. 3 108 ,. 1709 , . 13 18, . 102, MN 2 501 

X . 0895, .0807, .0742, . 0648 ,. 0582 ,. 0474 ,.0414, 2 502 

Y. 0351, .0320, .0304, .0295 ,. 0286 ,. 0283 / 2503 

DATA! XSECO ( 1 , 26 ) , I = 1, 20) / 1 . 83 ,1 .13 , .555 , . 344 , . 18 3 , . 138 , . 106, .0919, IRON 

X. 0828, .0762, .0664, . 0595, . 0485 ,. 0424 ,. 036 , . 03 3 , .031 3 , .0304, . 0295, Z= 26 

Y . 0294 / FE 

END 
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t I RFTC X2SEC0 DECK 


BLOCK DATA 

COMMON /STDATA/ * 

1 NAME (100) , ENR 3 (20) , XSECO < 20 , 1 00 ) , XS ECE A ( 20 ) , 

2 C ( A, 6 » 9 ) , E M I NO 3 ( 8 ) , NA MB ( 8 ) , XSECNT(IOO) ,DCONST 

DATA( XSECOC I ,27) ,1 = 1, 20) /2. 01 * 1 .235 , . 5999 , . 3669 , . 18 80 , . 1 38 7, . 10 52, CO 2701 
X. 0910, .0817, .0750, .0655, .0588,. 0479,. 0419, 2702 

V.0356, .03 27, .0312, .0304, .02 97, . 02 96/ 2703 

DATA( XSECO (I, 28), I =1,20) /2. 32, 1.427,. 6832,. 4127 ,.2550, .1488, .1108, NI 2801 

X. 0951, .0854, . 0783,. 0682,. 0611 ,. 0499,. 043 7, 2802 

Y. 0372, .0343, .0329, .0320, .0314, .0314/ 2803 

DATA( XSECO ( 1,29) ,1=1, 20) /2. 45,1. 51 ,. 713..427 , .206, . 147,. 108, .0916, COPPER 

X . 082, *0751, % C654, .058 5,. 04 76,. 041 8 , . 0357 , . 03 3 , .0 316 ,.0 309, .0303, Z = 29 

Y. 0305 / CU 

DAT A( XSECO( I ,30) ,1 = 1, 20) /2.71 ,1.661 ,.7773 ,. 4630,. 2182, *151 9, *10 95, ZN 3001 

X. 092 7, .08 27, . 0 756 ,. 0659 ,. 0590 , . 0481 ,.0422, 3002 

Y. 036 2, .0335, .0322,. 03 16,. 031 2,. 0313 / 3003 

DATA( XSECO ( I ,3 1 )»I = 1, 201/2. 88,1. 768,. 8150, .434 1,.2228 ,. 1527, .1074, GA 3101 

X . 0904, .0806, .0734, .0639,. 0571 ,. 0467, .0410, 310? 

Y. 0353, .03 27, .0 31 5,. 0309,. 03 07,. 0309/ 310 3 

DATA! XSECO ( 1 ,3 2 ) ,1 = 1, 20) /3. 12 ,1.905 ,. 8876, .5194 , . 2 3 23 , . 15 62 , . 1 0 8 1, GE 3 201 

X. 0904, .0802, .0731, .0634, .0568, .0463, .0407, 3202 

Y. 0351, .0327, . 0 31 5 , . 0 3 1 0 , . 03 08 ,. 0311 / 3203 

DATA! X SECO ( 1 , 3 3 ) , I = 1 , 20 ) / 3 . 3 9 , 2 . C8 , . 96 06 , . 55 9 2 , . 2440 , . 1 6 17 , , 1099, AS 3301 

X . 0910 , .08 06, .0 733,. 0 63 6,. 056 8,. 0464,. 0407, 3 302 

Y. O 35 3, .0329, .031 8,. 03 13,. 031 3,. 0316/ 3303 

DATA! X SEC 0(1, 34), 1 = 1, 20)/3.62 , 2. 20,1. 029, .5586, .2515, .1634, ,1087, SE 3401 

X. 0897, .0792, . 0 720 ,. 0623 ,. 0556 , . 0454 ,. 0399 , 3402 

Y. O 346, .03 24, .0315,. 0310,. 0310,. 0314/ 3403 

DATA ( XSECO! I ,35 ) ,1 =1, 20/3. 99,2.43 ,1.135 ,.6436 ,.2691, ,1727, ,1124, BR 3501 

X .0922 , .0810, .0 735,. 063 5,. 0567,. 0462 ,.0406 , 3502 

Y.0354, .03 3 2, .0 32 3,. 03 18,. 0320,. 032 5/ 3 50 3 

DATA( XSECO! I, 36), (=1,20/ 4. 24, 2. 56, 1.2 03, ,6791, .278 5, ,175 3, .1122, KR 3601 

X. 0913, .08 00, .0 727,. 0 625,. 0 557,. 0454,. 0400, 3602 

Y . 0349 , .0329, .0320, .0317,. 0318,. 0324/ 3603 
DATA! XSECO ( 1 , 37 ) , I = 1 , 20 ) / 4. 60 , 2 . 79 , 1 . 309 , , 733 5 , , 29 6 3 , , 1 826 , , 1 1 49, RB 3701 

X. 0 92 8, .0811, .0732,. 0630,. 0562 ,. 045 8,. 0404, 3702 

Y. C353, .0333, .0325, .0322,. 0325 ,. 0331 / 3703 

DATA! XSECO! I, 38), 1=1, 20)/ 4. 96, 3. 03, 1.406, ,7882, .3122, ,1900, ,1176, SR 3801 

X. 0941, .0817, .0736, .0633, .0564, .0459, .0405 , 3802 

Y. O 355, .0336, .0329,. 03 26,. 0330,. 033 7/ 3803 

DATA! XSECO! I, 39), I =1,20)/ 5. 40, 3.29, 1.542,. 85 55 , . 33 28 , . 200 2, , 1 2 1 2 , Y 3901 

X . 0960 , .0832, . 0743, . 0642 ,. 0 572 ,. 0465 ,.0411, 3902 

Y. O 361, .034 3, .0336,. 03 34,. 033 8,. 0346/ 390 3 

DATA! XSECO! I ,40) ,1 = 1, 20) / 5. 78,3. 52, 1.658,. 91 00 ,.35 17,. 2085, ,1239, ZR 4001 

X. f'978, .0839, .0 752, . 0644,. 05 73 ,. 0466 ,.0412 , 400 2 

Y . 0363 , .03 45, .0339 ,. 03 3 7, . 0 343 ,. 0351 / 4003 
DATA! XSECO! I, 41), I =1,20/ 6. 2 3, 3. 80, 1.764,. 98 20 ,. 37 2 4,. 2 178,. 12 75, NB 4101 

X . 0 996 , .08 50, .0 766,. 0650,. 0578,. 0470 ,.0416 , 4102 

Y . 0367 , .0350,. 0344,. 0343,. 0349,. 0359/ 4103 

DATA! XSECO! I , 4 2 ) , I =1 , 20 ) /6. 62 , 4. 04 , 1. 86 , 1. 03 , . 3 89 ,, 22 5 ,, 1 3 , . 0998, MOLYB- 
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X. 0851,. 0761, .0 648,. 057 5,. 0467,. 041 4 ,. 0365 ,. 0349 ,.0 344 , .0344, .0 349, DENUM 

Y . 0359 / z = 42 

OAT At XSECOt 1, 43), 1=1, 201/7.02, 4.28, 1.958, 1.089 , .4065, *2328, *1310, TC 4301 

' X. 1003, .0851, .0757, .0646, .0573, .0465, .0412, 4302 

Y.0365, .0349, . 0 345, . 0344 ,. 0353, .0363/ 4303 

DAT At X SECO (I, 44), 1=1, 20) /7. 43 ,4. 54 , 2. 075 ,1 .156 , .4274, . 2411 , *1327, RU 4401 
- X. 1011, .0855, .0759, .0645, .0571 ,.0464, .0411, 4402 

Y.0 365, .03 51, .0347,. 034 7,. 0356,. 0367/ 4403 

DATA! XSECOt 1,45 ),I =1, 20) /7. 99, 4. 89 , 2. 220 ,1 . 243 , .4548 , . 2 543 , . 1 373, RH 4501 

X. 1039, .0872, .0773, .0655,. 0579,. 0470, .0417, 4502 

Y. 0371 , .0357, .0354, .0354,. 0364 ,. 0376 / 4503 

DATA! XSECOC 1,46) , I *1, 20) / 8. 39 , 5. 1 2 ,2. 330 , 1 .3 0 1 , .47 10 , . 260 8 , . 1391 , PO 4601 

X. 1048, .0869, .0767, .0648,. 0573 ,. 0464, .0413, 4602 

Y. 0368, .0354, .0352,. 0353,. 0363,. 0375 / 4603 

DATA( X SECO t I ,47 ) , I = 1 , 20) / 8. 97 ,5 .49,2. 50 , 1 .388 , .50 1 4, ,2 759, .1447, AG 4701 

X . 10 76 , .0886, .0781, .0659, .0582, .0470, .0418, 4702 

Y. 0374, .0361, .0359,. 0360,. 0371 ,. 0384/ 4703 

DAT At XSECOt I ,48 ) , I =1 , 20 ) / 9. 31 , 5 . 71 ,2. 59, l . 43 7 , . 5168 , . 28 16 , . 1455, CD 4801 

X. 1065, .0878, .0771, .0648,. 0571,. 0462 ,.0410, 4802 

Y. 0368, .0356, . 0355, . 03 57,. 0368 ,. 0382 / 4803 

DATAt XSECOt I ,49 ) , I =1 , 20) /9. 85 ,6. 02 , 2. 75 , 1 . 515 , .5 4 3 6 , .2 9 40 , .149 2 , IN 4901 

X. 1083, .0887, .0778, .0652, .0574, .0464, .0412, 4902 

Y. 0370, .0359, . 0358, . 0360, . 03 72 ,. 0387/ 4903 

DATAt XSECOt 1 ,50), I = 1,20) /i 0.2, 6. 28, 2. 87, 1.58,. 563,. 303,. 153, .109, TIN 

X. 0886, .0776, . 0647, . 0568 ,. 0459,. 0408 ,. 0367 ,. 3 35 5 , .0355 ,. 0358, .0368, Z=50 

Y. 0383 / SN 

DATAt XSECOt I , 5 1 ) , I = l , 20) /1 0. 74,6. 57 ,3 . 01 ,1 . 650 , .5980 , . 3150 , . 1557, SB 5101 

X. 1101, .0892, .0780, .0647, .0568, .0458, .0407, 5102 

Y. 0367 , .0357, . 0357 ,. 0361 ,. 0374, . 0389/ 5103 

DATAt XSECOt 1 , 52 ) , I = 1, 20) /1 1 . 03 ,6. 75 ,3 . 1 0 , 1 . 700 , .6040 , . 319 1 , . 1557, TE 5201 

X. 1090, .0879, .0767, .0633, .0555, .0447, .0397, 5202 

Y . 0359, .0350, .0351, . 0354,. 0367,. 0383/ 5203 

DATAt XSECOt 1, 5 3), I =1,20) /1 1.9, 7. 26, 3. 34, 1.83,. 648,. 339,. 165,. 114, IOOINE 

X.0913, .0792, .0653,. 0571,. 046,. 04 09, . 037,. 036 ,. 0 361 , .0365 , .0377, Z=53 

Y • 0394 / I 

DATAt XSECOt I, 54), I =1,20)/ 12. 11, 7. 47 ,3. 46, 1.905,. 6640,. 3471,. 1679, XE 5401 

X. 1146, .0910, .0787, .0647, .0565, .0453 ,.0404, 5402 

Y. 0366, .0358, .0360, .0365,. 0379,. 0396/ 5403 

DATAt XSECOt 1,55 ), I *1,20) /1 2. 60, 7. 88 ,3. 67, 2. 01 2,. 7000,. 3642, .1745, CS 5501 

X. 1174, .0929, .0800,. 0655,. 0572 ,. 045 8,. 0408, 5502 

Y. 0371, .0363, .0366, . 0371 ,. 0385 ,. 0403 / 5503 

DATAt XSECOt I ,56), 1=1, 20) /I 3. 1 0 , 8. 1 6 ,3 . 79,2 .079 , . 7 180 , . 3719, . 1763, BA 5601 

X. 1174, .0928, .0798,. 0650,. 0566,. 0453 ,.0404, 5602 

Y. 0367, .0360, .0363,. 0369,. 0384,. 0402/ 5603 

DATAt XSECOt I ,57), 1=1, 20) /13. 69, 8. 55 ,4. 00 ,2. 1 95 , .7560, . 389 1, . 1839, LA 5701 

X. 1211, .0 9 47, .08 1 2,. 0 6 5 9,. 0 5 73,. 04 5 8,. 0408 , 5 7 0 2 

Y. 0372, .0365, .0369, .0375,. 0390,. 0409/ 5703 

DATAt XSECOt I ,5 8 ) ,1 =1 , 20) / 14. 2 8 , 8.99 ,4. 22 , 2. 3 1 2 , .7990 , . 408 1, , 1869, CE 5801 

X. 1249, .0973, .0828,. 0670,. 0582,. 0464,. 0414, 5802 

Y. 0377 , .0370, .0375, .0382, .0398, .0417/ 5803 

DATAt XSECOt I , 5 9 ) , 1=1 , 20) / 1 5. 06 ,9. 45 ,4. 49,2 .460 , . 8440 , . 428 1 , # 1989, PR 5901 

X . 1286, .0991, .0847, .0683, .0591, .0470, .0420, 5902 

Y. 03 83, .03 7 7, .0 3 8 3,. 03 9 0,. 04 07 ,. 042 7 / 5 9 0 3 

END 
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tIBFTC X3SEC0 DECK 


BLOCK DATA 

COMMON /STDATA/ ' 

1 NAME (100) , ENRG (20) , XSECO ( 20 , 1 00 ) , XS ECEA ( 20 ) , 

2 C (4,6,9) , EMI N3 3 ( 8) , NAMB18) , XSECNT(IOO) ,DCONST 

DATA( XSECO ( 1,60 ),I= 1,20) /I 5. 56, 9. 75,4. 64, 2.560,. 8720,. 4425,. 2037, NO 6001 

X. 1304, .1000, .0852, .0684,. 0591, .0469, .0419, 6002 * 

Y. 0383, .03 78, .0384,. 0391,. 04 08,. 0429/ 6003 

DATA! XSECO ( I , 6 1 ) , I =1 , 20) /I 6. 1 5 , 1 0. 1 4 ,4 .84 , 2 . 66 , .9 10 , .4578 , . 2093, PM 6101 

X. 1322, .1019, .0 861,. 068 8,. 0594,. 0470,. 0420, 6102 

Y. 0 3 84, .0380, . 0386, . 0394, . 041 1 ,. 0432 / 610 3 

DAT A( XSECO ( I ,62 ) ,1 = 1, 20) / 1 6. 54 , 1 0. 44 ,5 . 02 , 2 .76 , .940 , .4723 , . 2 140, SM 6201 

X. 1351, .10 2 8, .0867,. 0690,. 0594,. 0469,. 041 9, 6202 

Y. 0384, .0380, . 0387 ,. 0395 ,. 0412 ,. 0434 / 6203 

DATA( XSECO ( I ,63 ) , I = l, 20) /I 7. 23 ,1 0. 83 , 5. 24 , 2 .90 , .98 6 , .4924, . 221 6, EK 6301 

X. 1388, .1047, .0882, .0700, .0601, .04 74, .0424, 6302 

Y. 0389, .0385, .0393, . 0401 , . 041 9 ,. 0441 / 6303 

DATA( XSECO ( 1,6 4 ),I=1, 20) /2. 56 , 1 0. 93 , 5. 35 ,2 .96 , 1 . 0 1 5 , .50 3 1 , . 2245, GO 6401 

X. 139 7, .10 56, .0 880,. 069 5,. 05 96,. 0469,. 041 9, 6402 

Y. 03 85, .03 8 2, . 0390 , . 03 9 8, . 041 6 ,. 0439/ 6403 

DATA! XSECO! I ,65 ) ,1=1, 20 )/ 2. 64, 11. 32 ,5.57,3.09, 1.054,. 5 233, .2311, T8 6501 

X.1455, .1075, .0893, .0702, .0600, .0472 ,.0421, 6502 

Y .0388 , .03 84, .0393, . 0402, . 0420 ,. 0443 / 6503 

OATA( XSECO ( I, 66), 1=1, 20)/2.75, 11. 52 ,5. 76, 3. 20 , 1 .09 3 , . 5407 , .2368, OY 6601 

X. 1454, .1084, .0 900,. 0706,. 0602 ,. 0471 ,.042 2 , 660? 

Y. 0 389, .03 85, . 0395, . 040 4, . 0422 ,. 0446 / 6603 

DATA! X SEC0( 1, 67), 1 = 1, 20) /2. 91 , 11. 80, 5. 98, 3. 32, 1.131 ,.5 599, .2434, HO 6701 

X. 1492, .1103, .091 4,. 0 71 3,. 0607,. 04 74,. 0424, 6702 

Y. 0391, .0388, . 0398 , . 04 08 , . 042 6 , . 0450 / 6703 

DATA( XSECO! I, 68), 1=1, 20 J/3. 05, 12. 10, 6. 21, 3. 45, 1.180, .5 791,. 2510, ER 6801 

X. 1530, .1121, .0 928,. 07 20,. 061 2,. 0477,. 042 7, 680? 

Y. 0394, .0392, .0402 , . 041 2 , . 043 1 ,. 0455/ 6803 

DATA! XSECO( I ,69) ,1 = 1, 20) /3. 19,12.39,6.42 ,3 .58 ,1.228 ,.6012, .2586, TM 6901 

X. l 568, .1150, .0943, .0729,. 061 8,. 0481 ,.043 0, 6902 

Y. 0397, .03 95, . 0406 , . 04 1 6 , . 0435 ,. 0460/ 690 3 

DATA< XSECO ( I ,70) ,1 = 1,20) /3. 34,2. 01 ,6.60,3.70 ,1.257 , .6177, .263, YB 7001 

X. l 596, .11 59, .0951,. 0732,. 062 0,. 0481 ,.042 8, 7002 

Y. 0398, .0396, . 0406, . 04 1 7 , . 0436 , . 0462 / 7003 

0ATA( XSECO ( I ,71 ) ,1= l, 20) /3. 47,2. 13 ,6. 88,3.83, 1.315 , .6430, .2719, LU 7101 

X • 1634, .1188, .0 966,. 0 74 3,. 062 7,. 0496 ,.043 3, 7102 

Y .0402, .04 00, .04 11,. 042 1,-0441 ,. 0467/ 7103 

DATA! X SEC 0( 1 ,72), I = 1,20) /3.61 ,2. 19, 7. 04, 3. 92, 1.3 54, .6620, .2777, HE 7201 

X .166 2, .12 06, .0975,. 0 74 7,. 062 8,. 0485 ,.043 3 , 7 202 

Y.0402, .0401, .041 1, .042 2,. 0442 ,. 0469/ 7203 

DATA! XSECO! I ,73 ),I =1,20) /3. 76,2.26,7.2 9,4.08 , 1.40 2 ,.6 850, .2853, TA 7301 

X. 1700, .1225, .0994, .0757, . 0635, .0489, .0436, 730? 

Y. 0404, .0405, . 04 1 5, . 042 6 , . 044 7 , . 0474 / 7303 

DATA! XSECO! I ,74) ,1=1, 20) /3. 94,2. 34 ,7.49,4.21,1.44, .708 ,.293, .174, TUNGSTEN 

X. 125, . 101 , .0 76 3, .064,. 0492 ,. 043 7 ,. 0405 ,.04 02 ,.0409 , .0418, .0438, Z= 74 

Y .0465/ W 

DATA! XSECO! I, 75), 1 = 1, 20 1/4. 12, 2. 43, 7. 72, 4. 35 , 1 . 500 , . 729 , . 300 5, RE 7501 
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’ . 0484/ 

DAT A( XSECCM I ,82) ,1 = 1,201/5. 19,3.15,1.41,5.29 
X. 145, . 114, .0836,. 0684,. 051 2,. 0457, .0421,. 042 
Y .0489/ 

DAT A( XSECO ( I, 83), 1 = 1, 201/5.38 

X. 71 40. -148?. .1 165. .0850.. 0694 


DATAC XSECO I I, 85), 1=1, 20)/ 5. 85 

X. 2284, .1577, .1232, .0891 ,. 0722 

Y. 0442 , .0446, .0456, . 0464, . 0492 
DATA! XSEC0<I,86) ,1=1, 20)/5. 80 

X. 2246, .1549, . 1 203, . 0868, . 0701 

Y. 0428, .0432, .0441, .0449,. 0476 
DATAC XSECO ( I, 87), 1 = 1, 20 1/6. 05 

X. 2323, .1587, .1241, .0887, .0716 

Y. 0434, .0438, .0448, .0455,. 0483 
DATAC XSECOC I, 88), 1=1, 20) /6. 26 

X. 2371, .1626, .1261,. 0900,. 0723 

Y. 0436, .0440, .0448,. 0458,. 0485 
DATAC XSECOC I , 89 ) ,1 =1 , 20) /6. 52 

X. 2448, .1674, .1299, .0921, .0738 

Y. 0443 , .0446, .0456, .0456, .0492 
DATAC XSECOC I, 90), I = 1,20) /6. 70 

X. 2478, .1693, .1309,. 0925,. 0739 
✓-0447. .0444. .0457.-0464.-04^0 


3.29,1.46,5.46 


DATAC XSECOC I, 8 3), I = 1,20)/ 5. 38, 3 

X.2140, .1482, .1165, .0850, .0694, .0518, .0461, 
C.0429, .0432, .0442, . 0450, . 0477 ,. 0507 / 


.0523/ 
3.57,1.62,5.79 
. 0519, .0460, 
.0507/ 

3. 65. 1.68. . 920 
.052 8,. 0467, 
.0514/ 

3. 79. 1.72. . 930 
.0531 ,.0470, 
.0517/ 

3. 90. 1.80. . 970 
.0540,. 0477, 

. 0524/ 

4.06,1.82,1.00 
. 0539,. 0476, 
.0522/ 


,1.539, .750, .3063, 

7502 

7503 
OS 7601 

,1.587, .771, .3140, 

7602 

7603 
IR 7701 

,1.64, .795, .324, .191, 

7702 

7703 
PLATINUM 

,0418, .0427, .0448, 

7 = 78 

,1.694, .825, .3331, 

GO 7901 

,1.742, .848, .3398, 

7902 

7903 
HG 800 1 

,1.8, .866, .346, .204, 

8002 

8003 

THALLIUM 

,.0423, .0433, .0454, 

Z = 8 1 

,1.84,. 896, .356,. 208, 

THALLIUM 

LEAD 

,.0426, .0436, .0459, 

Z=82 

,1.908, .924, .3667, 

PB 

B I 8301 

,1.975, .962, .3772, 

8302 

8303 
PO 8401 

,2. 053, 1.001, .3916, 

8402 

8403 
AT 8501 

,2.025, .982, .3861, 

8502 

8503 
RN 8601 

,2. 103, 1.02, .3986, 

8602 
8603 
FR 8701 

,2. 161, 1.039, .4103, 

8702 

8703 
RA 8801 

,2.239 ,1.078, .4227, 

8802 
8803 
AC 8901 

,2. 269, 1.097, .4307, 

8902 

8903 
TH 9001 

,2. 376,1.146, .4489, 

9002 

9003 
PA 9101 

,2. 42, 1.17, .452, .259, 

9102 

9103 
URANIUM 

,.0446, .0455,. 0479, 

Z = 92 


U 
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QADH 

HEATING IN AXIAL TANK 


CONTROL A 

6 

2 

3 

A 

13 

13 

12 

-0 

1 1 

0 

-0 

-0 

SOURCE 3.3000E 

19 

-0. 


-0. 


-0. 


-0. 






R 0. 


1. 5000E 

01 

3. OOOOE 

01 

A. OOOOE 

01 

5. OOOOE 

01 





Z 0. 


3. OOOOE 

01 

6. OOOOE 

01 

9. OOOOE 

01 

1.0500E 

02 1.2000E 

02 

1.3000E 

02 


PHI 0. 


1.5700E 

00 

3.1A00E 

00 








ZONE 

BN OS COMP 

BND 

PK 1 

BND 

PK 

BND 

PK 

BND 

PK 

BND PK 

BND 

PK 


1 

3 1 

-1 

13 

7 

3 

2 

2 







2 

3 2 

- 2 

1 

7 

3 

3 

A 







3 

A A 

-1 

13 

10 

12 

3 

A 

-7 

1 





A 

3 A 

*3 

2 

10 

12 

A 

5 







5 

A 3 

-A 

A 

11 

6 

5 

8 

6 

11 





6 

5 A 

-A 

A 

10 

12 

5 

9 

-11 

5 

-12 7 




7 

3 3 

-A 

A 

12 

6 

5 

10 







8 

3 3 

- 5 

5 

8 

9 

6 

11 







9 

5 A 

-5 

6 

10 

12 

6 

11 

-8 

8 

-9 10 




10 

3 3 

-5 

7 

9 

9 

6 

11 







1 1 

-1 A 

-6 

8 











12 

-3 A 

-10 

3 

-1 

13 

6 

11 







13 

-1 A 

1 

1 











BN D 

EQ PARAMETERS. 

. . 











1 

6 -0. 


-0. 


-0. 


-0. 


-0. 

-0. 


0. 


2 

6 -0. 


-0. 


-0. 


-0. 


-0. 

-0. 


1.3000E 

02 

3 

6 -0. 


-0. 


-0. 


-0. 


-0. 

-0. 


1.6000E 

02 

A 

6 -0. 


-0. 


-0. 


-0. 


-0. 

-0. 


5. OOOOE 

02 

5 

6 -0. 


-0. 


-0. 


-0. 


-0. 

-0. 


7.5000E 

02 

6 

6 -o. 


-0. 


-0. 


—0 • 


-0. 

-0. 


1.6000E 

03 

7 

3 -0. 


-0. 


-0. 


0. 


0. 

-0. 


2.5000E 

03 

8 

3 -0. 


-0. 


-0. 


0. 


0. 

-0. 


3.6000E 

05 

9 

3 -0. 


-0. 


-0. 


1.2500E 

03 

0 . 

-0. 


3.6000E 

05 

10 

3 -0. 


-0. 


-0. 


0 . 


0 . 

-0. 


A. OOOOE 

06 

11 

2 1.000 OE 

00 

1. OOOOE 

00 

-1. OOOOE 

00 

0. 


0 . 

1.5000E 

02 

0 . 


12 

2 1 • OOOOE 

00 

1. OOOOE 

00 

-1. OOOOE 

00 

1.2500E 

03 

0 . 

1.5000E 

02 

0 . 
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QADH 


QADX 

PB EA BUILDUP IS USED 



0 

1.0000E-06 

1.0000E-C3 

9. 1429E- 

-14 

GRP 

EPS/E 

3ELTA/K 

A1 /BO 


1 

1.2470E 00 

1 • 5331E- 16 

1.0245E 

00 

2 

3. 2360E 00 

1 • 5672E-1 6 

1. 0167E 

00 

3 

2.4210E 00 

1.4593E-16 

1. 0083E 

00 

4 

1.7610E 00 

1.3450E-16 

1.0057E 

00 

5 

1. 282 0E 00 

1 • 243 8 E — 1 6 

1.0044E 

00 

6 

1. 0490E 00 

1.1652E-16 

1.0036E 

00 

7 

7. 75006- 01 

1.0910E-16 

1. 0031 E 

00 

8 

6.5600E-01 

1. 01546-16 

1.0025E 

00 

9 

5 • 0700E- 01 

9.64546-17 

1 • 0021 E 

00 

10 

3. 2600E-01 

9.1852E-17 

1.0016E 

00 

11 

1 • 8200E-0 1 

8 . 735 7E “17 

9. 9968E- 

-01 

12 

1.13006-01 

8.2 675E-1 7 

9. 9767E- 

-01 

13 

1. 2000E-01 

7.4865E-1 7 

9.9416E- 

-01 


COMP/GRP 

1 

2 

3 

1 

0 . 

1.5000E 00 

0 . 

2 

1. 1000E-01 

0 . 

5.4900E 00 

3 

7.0000E-C2 

0 . 

0 . 

4 

0 . 

0 . 

0 . 


2.9000E-01 

8.3000E-01 

5.8000E-01 



A2/B1 

A3/B2 

A4/B3 

EBAR 


4.5096E-O1 

—4.1 100E-02 

1 .48306-03 

3.0000E- 

•01 

5.4913E-01 

— 5.3537E-02 

2.2546E-03 

6.3000E- 

■01 

7.5082E-01 

-3.3812E-02 

1 .0 51 IE- 03 

l.iOOOE 

00 

6.7385E-01 

— 1 .7 537E-02 

5 .30456-04 

1.5500E 

00 

5 • 6403 E— 0 1 

-7.6500E-03 

3 .82 24E-04 

1.9900E 

00 

4.7451E-01 

-1.7261E-03 

3.601 5E-04 

2.3800E 

00 

4.008 2E-01 

2. 40056-03 

3.7699E-04 

2.7500E 

00 

3. 1748E— 0 1 

6. 51 42E-03 

4.2214E-04 

3.2500F 

00 

2.4947E-01 

9.5409E-03 
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I GAMMA NEUT TOTAL GAMMA NEUT TOTAL 

1 5.2868E 05 0. 5.2868E 05 1.5173E-01 0. 1.5173E-01 

2 9.4299E 05 0. 9.4299E 05 1.7326E-01 0. 1.7326E-01 

3 1.7602E 06 0. I.7602E 06 1.6902E-01 0. 1.6902E-01 

A 6.3I20E 06 0. 6.3120E 06 3.2597E-01 0. 3.2597E-01 

5 1.3868E 07 0. 1.3868E 07 3.5965E-01 0. 3.5965E-01 

6 2.0317E 07 0. 2.0317E 07 2.9389E-01 0. 2.9389E-01 

7 2.0893E 07 0. 2.0893E 07 I.9732E-0L 0. 1.9732E-01 

8 3.1554E 07 0. 3.1554E 07 2.1231E-01 0. 2.1231E-01 


HEATING IN AXIAL TANK 

ACCUMULATIVE VALUES AT TIME TIJ), SECONDS 

HEAT INPUT CALORIES TEMPERATURE RISE DEGREES CENT 

J TIJ I GAMMA NEUT TOTAL GAMMA NEUT TOTAL 

1 3.5633E 02 3.1554E 07 0. 3.I554E 07 2.I231E-01 0. 2.1231E-01 

2 5.9388E 02 5.2448E 07 0. 5.2448E 07 4.0963E-01 0. 4.0963E-0I 

3 8.3144E 02 7.2765E 07 0. 7.2765E 07 7.0352E-01 0. 7.0352E-01 

4 1.0096E 03 8.6633E 07 0. 8.6633E 07 1.0632E 00 0. 1.0632E 00 

5 1.1097E 03 9.2945E 07 0. 9.2945E 07 1.3891E 00 0. 1.3891E 00 

6 1.1469E 03 9.4705E 07 0. 9.4705E 07 1.5582E 00 0. 1.5582E 00 

7 1.1768E 03 9.5648E 07 0. 9.5648E 07 1.7314E 00 0. 1.7314E 00 

8 1.2000E C3 9,.6177E 07 0. 9.6177E 07 1.8831E 00 0. 1.8831E 00 
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Cleveland, Ohio, March 24, 1967, 
121-30-02-02-22. 
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(a) Conical-head cylindrical tank. 



(b) Elliptic-head conical tank. 
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(c) Elliptic-head, no tank, or no-head elliptic tank. 
Figure 3. - QADH and QADD tank geometries. 


n Q AD boundary number 
Q QAD region number 


Region 

Description 

1 

Core 

2 

Complete shadow shield 

5,7 

Liquid hydrogen 

10 

Spacecraft walls 

All others 

Void 



Figure 4. - QADD Core-shield propellant tank of spacecraft. 
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Boundary number 
Region number 


Region 

Description 

Tank emptying 
parameters 

Value 

1 

Bottom support 

JHEAD 

15 

2 

Core 

JTANK 

16 

3,4 

Reflector 

JWHED 

17 

5 

Pressure vessel 

JWTNK 

14 

7 

Top support 

ZPRIME 

530 

8,9 

Shield 

BURN 

1185 sec 

11, 12 

Liquid hydrogen 

Z-Slice 

531, 550, 575, 600, 700, 

17, 18, 19 

Tank wall 

coordinates 

800, 894, 1000, 1100, 1200, 

6, 10, 13 
14, 15, 16 

Void 

Outside regions 
(also void) 

(NSLICE = 20) 

1300, 1400, 1500, 1600, 1800, 
2000, 2200, 2400, 2600, 2800, 
3544 


KNUCKL 7 
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